Tổng quan nghiên cứu

Tỉnh Bình Định sở hữu 165 hồ chứa nước thủy lợi với tổng dung tích thiết kế đạt 583 triệu mét khối, phục vụ tưới tiêu, cắt lũ và cấp nước sinh hoạt cho hơn 1,6 triệu người dân trên diện tích tự nhiên 6.026 kilomet vuông. Trong tổng thể mạng lưới này, có 133 hồ đập đất có chiều cao dưới 15 mét (chiếm tỷ lệ 80,6%) và phần lớn được xây dựng bằng phương pháp thủ công trong giai đoạn từ năm 1975 đến 1990. Trải qua hơn 30 năm khai thác dưới tác động khắc nghiệt của thời tiết và ảnh hưởng từ biến đổi khí hậu, có tới 89/133 hồ (tỷ lệ 66,9%) bị hư hỏng, xuống cấp nghiêm trọng. Đáng báo động, hiện tượng thấm qua nền và thân đập xuất hiện tại 91/133 hồ (chiếm 68,4%), trở thành nguyên nhân hàng đầu gây sạt trượt mái, xói ngầm và đe dọa an toàn hồ chứa.

Mục tiêu cụ thể của nghiên cứu là đánh giá thực trạng hư hỏng, phân tích bản chất cơ học đất và cơ chế hình thành dòng thấm tại các công trình đập đất vừa và nhỏ trên địa bàn tỉnh Bình Định. Từ đó, tác giả nghiên cứu đề xuất các giải pháp kỹ thuật chống thấm tối ưu về mặt kinh tế - kỹ thuật và ứng dụng tính toán chi tiết cho công trình hồ chứa nước Mỹ Thuận.

Phạm vi nghiên cứu bao quát 165 công trình hồ chứa toàn tỉnh với trọng tâm thực địa tại hồ chứa nước Mỹ Thuận từ năm 1975 đến năm 2016. Nghiên cứu mang ý nghĩa thực tiễn cấp bách, cung cấp cơ sở khoa học giúp khống chế gradient thấm, kiểm soát đường bão hòa và duy trì hệ số an toàn ổn định công trình đạt mức tối thiểu từ 1,30 trở lên, trực tiếp phục vụ chương trình sửa chữa và nâng cấp 38 hồ chứa có nguy cơ mất an toàn cao do Ngân hàng Thế giới tài trợ.

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 ứng dụng hệ thống lý thuyết thủy lực công trình ngầm và cơ học đất chuyên sâu để mô tả chuyển động của chất lỏng trong môi trường lỗ rỗng:

  • Đị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 dính, đất cát mịn với vận tốc thấm tỷ lệ bậc nhất với gradient thủy lực theo phương trình vi phân cơ bản.
  • Các định luật thấm phi tuyến: Áp dụng lý thuyết thấm chảy rối của Xêzi-Kraxnoponxki và phương trình vi phân Prôni để mô phỏng sự vận động của nước ngầm qua các tầng cuội sỏi, đá nứt nẻ và vật liệu hạt thô khi gradient thủy lực tăng cao.
  • Lý thuyết thấm trong đất bão hòa và không bão hòa nhiều pha: Phân tích sự hình thành cột nước thủy lực tổng, ảnh hưởng của áp lực nước lỗ rỗng và lực hút dính lên độ bền chống cắt của đất.
  • Khái niệm then chốt: Nghiên cứu làm rõ 4 khái niệm nền tảng gồm hệ số thấm, đường bão hòa thấm, chiều cao thoát nước giới hạn và gradient thấm tới hạn chống xói ngầm.

Mô hình nghiên cứu được xây dựng trên nguyên lý phần tử hữu hạn hai chiều (FEM) kết hợp phương pháp cân bằng giới hạn lát cắt của Bishop và Morgenstern-Price để phân tích đồng thời trường áp lực thấm và ổn định trượt mái dốc.

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

Nghiên cứu sử dụng nguồn dữ liệu đa tầng gồm:

  • Toàn bộ hồ sơ khảo sát địa chất, khí tượng thủy văn của 133 hồ đập đất quy mô vừa và nhỏ do Sở Nông nghiệp và Phát triển nông thôn Bình Định quản lý.
  • Bộ chỉ tiêu cơ lý đất đắp và nền móng tại hồ chứa nước Mỹ Thuận (dung trọng khô, lực dính, góc ma sát trong, hệ số thấm k).

Cỡ mẫu nghiên cứu bao gồm 133 hồ chứa nhỏ trên toàn tỉnh, áp dụng phương pháp chọn mẫu phân tầng kết hợp chọn mẫu điển hình tập trung vào 38 hồ chứa hư hỏng nặng xuất hiện hiện tượng thấm qua thân đập.

Phương pháp phân tích được thực hiện trên tổ hợp phần mềm GeoStudio với mô-đun SEEP/W (tính toán thấm) và SLOPE/W (kiểm tra ổn định). Lý do lựa chọn giải pháp phần mềm này là khả năng mô phỏng chính xác trường thấm hai chiều phi tuyến qua môi trường vật liệu phân tầng không đồng nhất và giải quyết hiệu quả các điều kiện biên thủy lực phức tạp. Timeline nghiên cứu và phân tích số liệu được hoàn tất vào tháng 5 năm 2016.

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

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

Quá trình điều tra thực tế và phân tích mô phỏng đã đưa ra 4 phát hiện quan trọng:

  • Tỷ lệ hư hỏng do thấm chiếm ưu thế tuyệt đối: Trong 133 hồ đập nhỏ tại Bình Định, có 91 hồ bị thấm qua thân và nền (chiếm 68,4%), 31 hồ rò rỉ khớp nối cống lấy nước (chiếm 23,3%) và 82 hồ hỏng tràn xả lũ (chiếm 60%), trong đó 39 hồ bị hư hỏng rất nặng (chiếm 30%).
  • Lỗ hổng nghiêm trọng trong quản lý vận hành: Có tới 91,7% số hồ (122/133 hồ) không có hồ sơ thiết kế lưu trữ; chỉ 15% hồ (20/133 hồ) có phương án phòng chống lụt bão và chỉ 42,9% (57/133 hồ) thực hiện đăng ký an toàn đập theo Nghị định số 72/2007/NĐ-CP.
  • Thực trạng mất an toàn tại hồ Mỹ Thuận: Khảo sát mặt cắt D31 cho thấy mái thượng lưu chưa được gia cố, đường bão hòa dâng cao xuất lộ trực tiếp trên mái hạ lưu, làm giảm khả năng chống cắt của đất đắp và đẩy hệ số ổn định trượt xuống dưới ngưỡng an toàn thiết kế.
  • Hiệu quả của giải pháp tường nghiêng sân phủ: Kết quả mô phỏng trên mô-đun SEEP/W cho thấy giải pháp tường nghiêng kết hợp sân phủ bằng đất sét đầm chặt có hệ số thấm nhỏ $k = 7,5 \times 10^{-6}$ cm/s đã hạ thấp đường bão hòa về đáy đập, triệt tiêu nguy cơ xói ngầm và nâng hệ số an toàn ổn định trượt $K_{min}$ lên trên 1,30.

Thảo luận kết quả

Nguyên nhân cốt lõi gây ra các sự cố thấm đập tại Bình Định xuất phát từ quá trình thi công đắp đất thủ công trước năm 1990 không kiểm soát được dung trọng khô thiết kế. Tình trạng đất đắp không đồng nhất, độ chặt không đều (dung trọng khô thực tế chỉ đạt 1,4 đến 1,5 T/m3 thay vì 1,84 T/m3) dẫn đến lún sụt không đều khi ngập nước. Khi kết hợp với loại đất có tính tan rã và trương nở mạnh của miền Trung, dòng thấm nhanh chóng tạo thành các rãnh xói ngầm xuyên thân đập.

Hiện tượng này hoàn toàn tương đồng với các bài học sự cố vỡ đập Suối Trầu (4 lần gặp sự cố trong giai đoạn 1977-1983) và đập Suối Hành (vỡ năm 1986 tại Khánh Hòa làm hư hỏng trên 100 ha đất canh tác) do bỏ sót thí nghiệm tan rã và xác định sai dung trọng đất đắp.

Các kết quả tính toán có thể được trực quan hóa thông qua bảng so sánh kinh tế - kỹ thuật giữa phương án tường nghiêng sân phủ và tường hào bentonite, kết hợp biểu đồ tương quan giữa mực nước thượng lưu với độ cao thoát nước giới hạn ($MNTL \sim a_{gh}$) tại các mặt cắt D14, D31, D38. Biểu đồ cho thấy khi mực nước thượng lưu dâng cao, điểm xuất lộ của dòng thấm được khống chế an toàn trong phạm vi lăng trụ thoát nước hạ lưu.

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

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

  • Thi công giải pháp tường nghiêng kết hợp sân phủ hoặc chân răng thượng lưu: Áp dụng đất sét đầm chặt đạt hệ số thấm $k \le 10^{-5}$ cm/s cho 38 hồ chứa bị thấm xung yếu tại Bình Định trong giai đoạn 2016-2020. Giải pháp này giúp tận dụng vật liệu địa phương, giảm 20% đến 30% chi phí xây lắp so với phương án khoan phụt cao áp. Chủ thể thực hiện: Ban Quản lý Dự án Nông nghiệp và Phát triển nông thôn tỉnh Bình Định.
  • Tiêu chuẩn hóa quy trình tính toán kiểm tra thấm bằng phần mềm GeoStudio SEEP/W: Bắt buộc các đơn vị tư vấn áp dụng mô hình phần tử hữu hạn để kiểm toán 100% các mặt cắt xung yếu, đảm bảo hệ số ổn định trượt $K_{min} \ge 1,30$ đối với mọi trường hợp mực nước dâng bình thường và mực nước rút đột ngột. Timeline hoàn thành trước khi phê duyệt hồ sơ thiết kế bản vẽ thi công.
  • Nâng cấp và số hóa hồ sơ an toàn đập: Hoàn thành việc đo đạc, lập lại hồ sơ kỹ thuật cho 122 hồ chứa còn thiếu; lắp đặt thiết bị quan trắc đường bão hòa và đo áp lực thấm định kỳ hàng tháng cho 44 hồ chứa đã sửa chữa. Chủ thể thực hiện: Chi cục Thủy lợi phối hợp UBND các huyện, thị xã.
  • Ứng dụng công nghệ tường hào bentonite và khoan phụt vữa sét - xi măng: Áp dụng đối với các tuyến đập có tầng thấm ngầm dày trên 5 mét hoặc nền đá nứt nẻ phức tạp, đồng thời tổ chức tập huấn kỹ thuật quản lý vận hành đập định kỳ hàng quý cho 100% cán bộ thủy nông cấp cơ sở.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  • Các kỹ sư thiết kế và đơn vị tư vấn thủy lợi: Khai thác quy trình tính toán trường thấm hai chiều bằng phần mềm SEEP/W và các thông số cơ lý đất đắp điển hình duyên hải Nam Trung Bộ để tối ưu hóa mặt cắt đập cho 165 hồ chứa vừa và nhỏ.
  • Ban quản lý dự án và các nhà thầu xây dựng: Ứng dụng quy trình kỹ thuật thi công đắp đất chống thấm, xử lý nối tiếp cổ cống lấy nước và thi công chân khay chống thấm trên nền địa chất phức tạp nhằm đảm bảo hệ số an toàn $K \ge 1,30$.
  • Cơ quan quản lý nhà nước về thủy lợi (Sở Nông nghiệp & PTNT, Chi cục Thủy lợi): Sử dụng các biểu đồ quan hệ mực nước và độ cao thoát nước giới hạn ($MNTL \sim a_{gh}$) làm tiêu chuẩn kiểm tra an toàn đập định kỳ trước mùa mưa bão hàng năm, phục vụ chương trình hiện đại hóa 38 hồ đập xung yếu.
  • Giảng viên, học viên cao học và sinh viên chuyên ngành Công trình Thủy: Tiếp cận case study thực tế về mô hình hóa dòng thấm bão hòa và không bão hòa kết hợp phân tích các bài học sự cố vỡ đập điển hình tại Việt Nam.

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

Nguyên nhân chính dẫn đến hiện tượng thấm nghiêm trọng tại các đập đất nhỏ ở Bình Định là gì?
Nguyên nhân chủ yếu do 80,6% đập đất được đắp thủ công từ những năm 1975-1990, đất đắp không đạt dung trọng khô thiết kế (chỉ đạt 1,4-1,5 T/m3) và thiếu thí nghiệm tan rã, trương nở của đất. Khi tích nước, đất bị lún không đều tạo ra các rãnh xói ngầm làm 68,4% số hồ xuất hiện thấm qua thân và nền.

Tại sao phương pháp phần tử hữu hạn bằng SEEP/W lại ưu việt hơn các công thức giải tích cổ điển?
Mô hình SEEP/W cho phép phân tích dòng thấm hai chiều phi tuyến trong môi trường đất phân tầng không đồng nhất, mô phỏng chính xác đường bão hòa trong cả trạng thái bão hòa và không bão hòa. Phương pháp này giúp xác định chính xác vị trí xuất lộ dòng thấm tại mái hạ lưu với sai số thấp hơn nhiều so với giải tích phẳng.

Trong trường hợp nào nên lựa chọn giải pháp tường nghiêng sân phủ thay vì tường hào bentonite?
Giải pháp tường nghiêng kết hợp sân phủ bằng đất sét đầm chặt thích hợp cho các công trình sửa chữa có chiều dày tầng thấm đáy $T > 5,0$ mét và dồi dào mỏ vật liệu đất dính gần công trình. Giải pháp này giúp tiết kiệm 20-30% chi phí xây dựng so với tường hào bentonite mà vẫn đảm bảo hệ số an toàn $K \ge 1,30$.

Tiêu chí độ cao thoát nước giới hạn ($a_{gh}$) có ý nghĩa như thế nào trong giám sát an toàn đập?
Độ cao thoát nước giới hạn ($a_{gh}$) là điểm cao nhất mà đường bão hòa thấm được phép xuất lộ trên thiết bị thoát nước hạ lưu. Nếu đường bão hòa thực tế vượt quá ngưỡng $a_{gh}$, dòng thấm sẽ tràn ra mái đất hạ lưu gây trôi hạt mịn, xói lở và dẫn đến trượt mái đập nguy hiểm.

Bài học đắt giá nhất từ sự cố vỡ đập Suối Trầu và Suối Hành là gì?
Bài học cốt lõi là không được xem nhẹ công tác khảo sát cơ lý đất đắp. Sự cố Suối Trầu (vỡ 4 lần) do chọn sai dung trọng thiết kế 1,5 T/m3, còn đập Suối Hành vỡ năm 1986 do bỏ sót chỉ tiêu tan rã và lún ướt của đất, dẫn đến phá hủy thân đập chỉ sau vài giờ tích nước.

Kết luận

  • Luận văn đã phân tích toàn diện hiện trạng 165 hồ chứa tại Bình Định, xác định 68,4% số hồ bị thấm và nhận diện 38 hồ chứa xung yếu cần xử lý khẩn cấp.
  • Thiết lập cơ sở lý thuyết vững chắc về dòng thấm hai chiều bão hòa và không bão hòa, ứng dụng thành công phần mềm SEEP/W trong mô phỏng trường thấm đập đất.
  • Khẳng định giải pháp tường nghiêng kết hợp sân phủ/chân răng bằng đất sét là giải pháp kỹ thuật tối ưu, giúp nâng hệ số ổn định trượt tại hồ Mỹ Thuận đạt $K_{min} \ge 1,30$.
  • Xây dựng bộ biểu đồ tiêu chí thấm ($MNTL \sim a_{gh}$) phục vụ công tác quan trắc và cảnh báo an toàn đập trong suốt quá trình vận hành khai thác.
  • Đề xuất lộ trình nâng cấp an toàn đập giai đoạn 2016-2025 gắn liền với nguồn vốn tài trợ của Ngân hàng Thế giới và ngân sách địa phương.

Các chủ hồ, đơn vị tư vấn và cơ quan quản lý thủy lợi hãy áp dụng ngay bộ tiêu chí giám sát thấm và giải pháp mô phỏng SEEP/W từ công trình hồ Mỹ Thuận để đảm bảo an toàn tuyệt đối cho hệ thống hồ đập trước mùa mưa bão.