Tổng quan nghiên cứu
Tỉnh An Giang nằm ở khu vực đầu nguồn Đồng bằng sông Cửu Long, chịu tác động trực tiếp từ hai nhánh sông lớn là sông Tiền và sông Hậu với lưu lượng dòng chảy bình quân đạt 13.800 m³/s trong năm và lên đến 24.000 m³/s vào mùa lũ. Nhằm phục vụ sản xuất nông nghiệp trên tổng diện tích 280.000 ha, địa phương đã xây dựng hệ thống công trình gồm 1.069 tuyến đê bao với tổng chiều dài 5.384 km, cao trình đỉnh đê dao động từ +2,8 m đến +6,0 m. Hệ thống này đảm nhận nhiệm vụ kiểm soát lũ cho 624 tiểu vùng với quy mô 247.183 ha, đồng thời kết hợp làm đường giao thông nông thôn và cụm tuyến dân cư.
Tuy nhiên, theo Báo cáo số 196/BC-SNN&PTNT của Sở Nông nghiệp và Phát triển nông thôn tỉnh An Giang, chỉ riêng trong năm 2017 đã ghi nhận 66 tuyến đê bao bị sạt lở với tổng chiều dài 36.200 m, đe dọa ngập úng nghiêm trọng cho 18.400 ha đất canh tác lúa vụ ba. Vấn đề nghiên cứu trọng tâm xuất phát từ thực trạng thi công đắp đê mang tính tự phát: người dân và các đơn vị quản lý thường tận dụng bùn sét nạo vét từ đáy kênh đổ trực tiếp lên đỉnh đê mà không qua đầm nén kỹ thuật.
Mục tiêu cụ thể của luận văn là làm rõ cơ chế phá hoại, xác định nguyên nhân gây sạt trượt sâu và xây dựng mô hình dự báo ổn định mái dốc dưới tác động đồng thời của tải trọng đắp bùn, chiều sâu nạo vét và biến động thủy triều. Nghiên cứu giới hạn không gian tại các tuyến đê bao xung yếu thuộc huyện Phú Tân với cao trình đỉnh đê trên +4,0 m trên nền đất sét yếu bão hòa, khai thác chuỗi dữ liệu thủy văn thực đo 20 năm từ năm 1998 đến năm 2017 tại trạm Tân Châu. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc xác lập hệ số an toàn ổn định công trình đạt ngưỡng quy chuẩn thiết kế tối thiểu 1,25 theo QCVN 04-05:2012/BNNPTNT, giúp giảm thiểu rủi ro vỡ đê và thiệt hại kinh tế hàng tỷ đồng mỗi năm.
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 phương pháp cân bằng giới hạn (Limit Equilibrium Method - LEM) làm nền tảng toán học cốt lõi để phân tích độ ổn định mái dốc. Trong các mô hình phân mảnh cổ điển, thuật toán Bishop đơn giản hóa (1955) được lựa chọn nhờ khả năng thiết lập cân bằng lực tổng thể và mô men một cách chuẩn xác, mô tả mặt trượt nguy hiểm dạng cung tròn đi sâu vào tầng đất yếu bên dưới thân đê. Độ ổn định của công trình được định lượng thông qua hệ số an toàn Fs, biểu thị tỷ số giữa tổng mô men kháng trượt do sức bền chống cắt của đất tạo ra và tổng mô men gây trượt do trọng lượng khối đất sinh ra.
Về mô hình ứng xử vật liệu, luận văn áp dụng định luật Mohr-Coulomb kết hợp chặt chẽ với lý thuyết sức kháng cắt không thoát nước (Undrained Shear Strength - Su) dành riêng cho đất sét yếu bão hòa theo Tiêu chuẩn ngành 22 TCN 262-2000. Ba khái niệm chuyên ngành then chốt được tích hợp đồng bộ gồm:
- Sức kháng cắt nguyên trạng Su xác định từ thí nghiệm cắt cánh hiện trường (VST).
- Hệ số hiệu chỉnh sức kháng cắt mu theo lý thuyết của Bjerrum (1973), biến thiên từ 0,70 đến 1,09 tương ứng với chỉ số dẻo Ip từ 10 đến 70 nhằm bù trừ hiệu ứng bất đẳng hướng và tốc độ trượt cắt trong đất.
- Công thức tương quan ứng suất hiệu dụng của Ladd và Foote (1974) với lực dính tính toán Cu = 0,25 nhân với ứng suất pháp hữu hiệu theo phương thẳng đứng đối với các lớp đất bề mặt bị khô nứt.
Phương pháp nghiên cứu
Nguồn dữ liệu của luận văn được tổng hợp từ hồ sơ khảo sát địa chất công trình và chuỗi quan trắc mực nước tự động tại trạm Tân Châu trong giai đoạn 1998 - 2017. Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu đã tiến hành thu thập và phân tích 6 mẫu đất nguyên dạng trên mỗi hố khoan khảo sát sâu 12 m, với khoảng cách bố trí hố khoan từ 250 m đến 500 m dọc theo các tuyến kênh sạt lở thực tế như Kênh Km 5 và Kênh 26/3 theo đúng Tiêu chuẩn quốc gia TCVN 10404:2015 và TCVN 9153:2012.
Lý do lựa chọn phần mềm chuyên dụng GeoStudio với mô đun SLOPE/W phiên bản 2012 là vì công cụ này giải quyết xuất sắc bài toán lặp đa biến tìm mặt trượt nguy hiểm nhất trên nền địa chất phân tầng phức tạp. Khác với phương pháp phần tử hữu hạn (như Plaxis 2D) đòi hỏi các thông số thí nghiệm nén ba trục phức tạp và tốn kém vốn hiếm khi được thực hiện ở các dự án thủy lợi cấp huyện, SLOPE/W cho phép tận dụng tối đa chỉ tiêu cắt cánh hiện trường Su và độ ẩm tự nhiên, mang lại kết quả hội tụ nhanh chóng, chuẩn xác và phù hợp tuyệt đối với điều kiện công trình nông thôn. Toàn bộ tiến trình mô phỏng, kiểm định thông số đầu vào và tính toán lặp được thực hiện trong thời gian 8 tháng từ tháng 9 năm 2017 đến tháng 4 năm 2018.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Phân tích mô hình số trên các mặt cắt điển hình đã chỉ ra bốn phát hiện khoa học mang tính quyết định:
- Thứ nhất, hoạt động nạo vét mở rộng và hạ sâu lòng kênh đơn thuần chỉ làm suy giảm nhẹ mức độ an toàn của mái đê. Khi hạ cao độ đáy kênh xuống 2,3 m, hệ số an toàn Fs chỉ suy giảm trong khoảng từ 0,05 đến 0,10 (tương đương mức giảm biên độ dưới 8%), chưa đủ để kích hoạt hiện tượng trượt cung tròn sâu nếu thân đê không bị gia tải.
- Thứ hai, việc đắp bùn sét nạo vét lên đỉnh đê mà không lu lèn là tác nhân phá hoại nghiêm trọng nhất. Việc nâng chiều cao đỉnh đê thêm 1,0 m bằng bùn ướt có dung trọng 17,3 kN/m³ khiến hệ số an toàn Fs sụt giảm đột ngột đến 0,35 đơn vị (giảm hơn 26% độ ổn định), kéo hệ số an toàn từ mức đạt chuẩn 1,32 xuống dưới ngưỡng giới hạn phá hoại 1,00.
- Thứ ba, tồn tại một ngưỡng hình học tới hạn: độ chênh cao 7,0 m giữa đỉnh đê bao và đáy lòng kênh là điểm giới hạn gây ra hiện tượng sạt trượt sâu trên hầu hết các loại nền sét yếu bão hòa tại khu vực huyện Phú Tân.
- Thứ tư, chế độ thủy văn và biên độ dao động triều trong ngày tác động trực tiếp đến áp lực thủy tĩnh giữ mái dốc. Khi mực nước rút nhanh từ đỉnh lũ cao trình +4,5 m xuống mức kiệt âm dưới đáy sông với biên độ triều đạt 1,41 m, hệ số ổn định của bờ đê giảm mạnh từ 20% đến 35%, biến các đoạn đê đắp bùn thành những khối trượt tiềm năng.
Thảo luận kết quả
Cơ chế gây mất ổn định trượt sâu được làm sáng tỏ dựa trên sự thay đổi trạng thái ứng suất nội tại. Bùn nạo vét từ đáy sông có độ bão hòa nước cao và độ sệt lớn, khi đổ lên đỉnh đê không chỉ làm gia tăng tải trọng gây trượt thẳng đứng mà còn thẩm thấu nước vào thân đê cũ, triệt tiêu lực dính đơn vị Cu và giảm góc ma sát trong của đất nền. Kết quả này hoàn toàn tương đồng với các công bố quốc tế của VandenBerge và cộng sự (2013) về sự suy giảm 70% sức kháng cắt không thoát nước khi gia tải đột ngột trên mái dốc, đồng thời giải thích rõ diễn tiến sạt lở hình thành hàm ếch và trượt khối theo mô hình của Thorne và Lewin (1979).
Về mặt trình bày dữ liệu, các mối quan hệ tương quan phi tuyến giữa hệ số an toàn Fs và chiều cao đắp đê được trực quan hóa tối ưu thông qua hệ thống đồ thị đường cong suy giảm dốc đứng khi chiều cao đắp vượt quá 0,5 m. Bên cạnh đó, các bảng tra cứu mặt cắt hình học cho thấy sự vượt trội của mái dốc thiết kế m = 1,5 so với mái dốc truyền thống m = 1,0: hệ số mái dốc m = 1,5 giúp mở rộng bán kính cung trượt nguy hiểm, nâng hệ số an toàn trung bình lên thêm 18% trong cùng điều kiện địa chất và tải trọng bùn thải.
Đề xuất và khuyến nghị
Dựa trên kết quả mô phỏng định lượng, luận văn đưa ra bốn nhóm giải pháp hành động cụ thể nhằm bảo vệ an toàn công trình đê bao:
- Chấm dứt triệt để tập quán đổ bùn nạo vét trực tiếp lên đỉnh đê hiện hữu: Ban Quản lý dự án nông nghiệp và các nhà thầu thi công phải khống chế chiều cao đắp đất bùn không quá 0,5 m cho mỗi chu kỳ nạo vét, bắt buộc triển khai giải pháp phơi khô và đầm nén cơ giới từng lớp đạt hệ số đầm chặt K tối thiểu 0,85 trong vòng 30 ngày kể từ khi thi công nhằm ngăn chặn sự gia tăng tải trọng trượt.
- Chuẩn hóa thông số hình học mặt cắt đê thiết kế: Các đơn vị tư vấn thiết kế công trình thủy lợi cần chuyển đổi toàn bộ hệ số mái dốc taluy từ m = 1,0 sang m = 1,5 hoặc m = 2,0 đối với các tuyến kênh có cao trình đỉnh đê trên +4,0 m, đồng thời khống chế khoảng cách chênh cao an toàn giữa đỉnh đê và đáy kênh luôn nhỏ hơn 6,0 m, hoàn thành cập nhật hồ sơ thiết kế trước quý 2 năm 2019.
- Quy hoạch lại quy trình nạo vét thông luồng: Cơ quan quản lý đường thủy và thủy lợi cần giới hạn chiều sâu đào vét mỗi đợt không vượt quá 1,5 m so với đáy tự nhiên, đồng thời thiết lập hành lang lưu không cách chân mái đê tối thiểu 5,0 m để bố trí bãi tập kết bùn thải riêng biệt có bờ bao ngăn bùn chảy tràn.
- Xây dựng hệ thống cảnh báo và quan trắc sạt lở tự động: Chi cục Thủy lợi tỉnh An Giang phối hợp cùng chính quyền địa phương lắp đặt mốc quan trắc chuyển vị tại 66 vị trí xung yếu có nguy cơ sạt trượt cao, tổ chức cắm biển cảnh báo và điều tiết giao thông tải trọng nặng qua thân đê trong các khung giờ triều rút kiệt có biên độ dao động trên 1,41 m trước mùa mưa lũ năm 2020.
Đối tượng nên tham khảo luận văn
Công trình nghiên cứu mang giá trị thực tiễn và học thuật sâu sắc, đặc biệt hữu ích cho bốn nhóm đối tượng chuyên môn:
- Các cơ quan quản lý nhà nước về nông nghiệp và đê điều: Điển hình là Sở Nông nghiệp và Phát triển nông thôn, Chi cục Thủy lợi tỉnh An Giang và các tỉnh vùng bán đảo Cà Mau, Đồng Tháp Mười. Dữ liệu luận văn cung cấp căn cứ khoa học xác đáng để quy hoạch bảo vệ 624 tiểu vùng kiểm soát lũ và ban hành quy định kỹ thuật về nạo vét kênh mương thủy lợi nội đồng.
- Các công ty tư vấn thiết kế và thẩm tra công trình xây dựng thủy lợi: Tài liệu cung cấp trọn bộ quy trình tính toán ổn định trượt sâu bằng phương pháp Bishop kết hợp hiệu chỉnh sức kháng cắt Su theo chỉ số dẻo Ip, giúp kỹ sư thiết kế lựa chọn mặt cắt hình học đê tối ưu với hệ số mái dốc m = 1,5 thay vì áp dụng máy móc các quy chuẩn cũ.
- Các đơn vị nhà thầu thi công nạo vét và cải tạo đê bao: Luận văn là cẩm nang hướng dẫn thực tế về giới hạn tải trọng đất đắp, giúp chỉ huy trưởng công trường xác định đúng vị trí đổ đất thải và chiều dày lớp đắp an toàn nhằm tránh gây ra các sự cố trượt mái đê trong quá trình nạo vét.
- Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Xây dựng công trình: Đây là nguồn tài liệu tham khảo mẫu mực về phương pháp kết hợp giữa lý thuyết cơ học đất hiện đại, phần mềm mô phỏng địa kỹ thuật SLOPE/W và số liệu kiểm chứng hiện trường trong điều kiện đất yếu Đồng bằng sông Cửu Long.
Câu hỏi thường gặp
Tại sao việc đắp đất bùn nạo vét lên đỉnh đê lại gây nguy cơ sạt trượt cao hơn việc nạo vét sâu lòng kênh? Bởi vì việc nâng cao đỉnh đê thêm 1,0 m bằng bùn sét bão hòa làm gia tăng trực tiếp mô men gây trượt lên đỉnh cung trượt nguy hiểm, khiến hệ số an toàn Fs giảm mạnh tới 0,35 đơn vị. Trong khi đó, việc nạo vét hạ sâu đáy kênh thêm 2,3 m chỉ làm giảm Fs khoảng 0,08 đơn vị do tác động cục bộ ở chân khay mái dốc.
Ngưỡng chênh cao giới hạn giữa đỉnh đê bao và đáy kênh để đảm bảo an toàn chống trượt sâu là bao nhiêu? Mô hình tính toán xác định độ chênh cao 7,0 m là ngưỡng phá hoại tới hạn khiến hệ số an toàn Fs tụt xuống dưới 1,00. Do đó, trong thiết kế và quản lý vận hành, độ chênh cao hình học giữa đỉnh đê và cao độ đáy kênh bắt buộc phải được khống chế dưới 6,0 m nhằm duy trì hệ số an toàn Fs đạt mức yêu cầu từ 1,25 trở lên.
Ưu điểm nổi bật của việc áp dụng mô hình Bishop và phần mềm SLOPE/W trong nghiên cứu này là gì? Phương pháp Bishop trong SLOPE/W giải quyết tối ưu bài toán trượt tròn trên nền đất sét yếu bão hòa nước theo chuẩn Tiêu chuẩn ngành 22 TCN 262-2000. Phương pháp này chỉ yêu cầu các chỉ tiêu thí nghiệm cắt cánh Su và chỉ số dẻo Ip dễ thu thập, khắc phục triệt để sự thiếu hụt thông số thí nghiệm nén ba trục đắt đỏ của phương pháp phần tử hữu hạn.
Biến động thủy triều ảnh hưởng như thế nào đến sự ổn định của mái đê bao? Khi thủy triều rút kiệt với biên độ dao động lên đến 1,41 m trong ngày, áp lực nước thủy tĩnh tác dụng giữ mái ngoài của đê bị triệt tiêu đột ngột trong khi nước trong thân đê chưa kịp thoát ra ngoài. Áp lực nước lỗ rỗng dư này làm giảm đáng kể sức kháng cắt của khối đất, khiến hệ số an toàn Fs suy giảm từ 20% đến 35%.
Kết quả nghiên cứu này có thể áp dụng cho các khu vực khác ngoài huyện Phú Tân không? Hoàn toàn có thể áp dụng cho toàn bộ các huyện thuộc tỉnh An Giang và các tỉnh lân cận tại Đồng bằng sông Cửu Long có điều kiện địa chất trầm tích sét yếu bão hòa tương đồng và đê bao có cao trình trên +4,0 m, với điều kiện kỹ sư cần cập nhật chính xác chỉ số dẻo Ip và số liệu thủy văn cục bộ của từng tiểu vùng.
Kết luận
- Luận văn đã xác định chính xác nguyên nhân cốt lõi gây ra các sự cố sạt trượt sâu đê bao tại An Giang là do hoạt động đắp bùn sét nạo vét không qua đầm nén lên đỉnh đê, làm giảm hệ số an toàn Fs tới 0,35 đơn vị.
- Nghiên cứu đã tìm ra ngưỡng hình học tới hạn với độ chênh cao 7,0 m giữa đỉnh đê và đáy kênh, đồng thời chỉ rõ sự suy giảm độ ổn định từ 20% đến 35% khi mực nước rút kiệt.
- Công trình đã kiểm chứng thành công tính thực tiễn của phương pháp cân bằng giới hạn Bishop kết hợp công thức hiệu chỉnh sức kháng cắt Su theo Bjerrum và Ladd & Foote trên phần mềm chuyên dụng SLOPE/W.
- Đề tài đã đề xuất thành công bộ giải pháp kỹ thuật toàn diện, trọng tâm là chuyển đổi mặt cắt mái dốc từ m = 1,0 sang m = 1,5 và kiểm soát tải trọng đắp thải nhằm bảo vệ 1.069 tuyến đê toàn tỉnh.
- Về lộ trình tiếp theo, các cơ quan chức năng cần hoàn thiện việc cập nhật quy chuẩn thiết kế đê bao trước quý 2 năm 2019 và hoàn thành lắp đặt hệ thống cảnh báo sạt lở giai đoạn 2019 - 2021.
Hãy áp dụng ngay các thông số hình học và quy trình kiểm soát tải trọng đắp bùn từ luận văn này vào công tác lập dự án, thiết kế và thi công đê bao thủy lợi để bảo vệ vững chắc hệ thống công trình nông nghiệp trước diễn biến biến đổi khí hậu ngày càng phức tạp.