Tổng quan nghiên cứu

Lưu vực sông Đáy giữ vị trí chiến lược trọng yếu đối với vùng đồng bằng sông Hồng với diện tích lưu vực rộng hơn 7.500 km² và chiều dài dòng chính khoảng 247 km. Trên địa bàn thành phố Hà Nội, tuyến đê tả Đáy đoạn qua huyện Ứng Hòa có chiều dài 36,322 km (từ K43+700 giáp Thanh Oai đến K80+022 giáp Kim Bảng, Hà Nam), đóng vai trò then chốt trong việc bảo vệ an toàn dân sinh và cơ sở hạ tầng nông nghiệp cho hơn 179.000 người dân thuộc 29 đơn vị hành chính. Tuy nhiên, sự biến đổi bất thường của khí hậu cùng chế độ thủy văn cực đoan đã làm gia tăng độ chênh lệch mực nước giữa mùa lũ và mùa kiệt, với biên độ đỉnh lũ lịch sử đo được đạt 7,68 m và mực nước kiệt thiết kế tần suất P=95% hạ thấp xuống mức +0,74 m.

Tình trạng sạt lở bờ tả sông Đáy diễn ra đặc biệt nghiêm trọng sau các cơn bão số 10, số 11 và đợt mưa lũ lịch sử vào tháng 10 năm 2017. Trọng điểm sạt trượt xuất hiện tại đoạn qua xã Vạn Thái (từ K67+000 đến K68+000) với chiều dài cung sạt liên tục kéo dài 500 m đến 600 m. Mép sạt lở ăn sâu vào bãi sông, cuốn trôi diện tích canh tác hoa màu và chỉ còn cách chân đê tả Đáy 38 m, đồng thời đe dọa trực tiếp sự an toàn của tuyến Quốc lộ 21B cùng các khu dân cư lân cận với khoảng cách nguy hiểm từ 5 m đến 10 m.

Luận văn thạc sĩ chuyên ngành Kỹ thuật xây dựng công trình thủy tập trung nghiên cứu, đánh giá toàn diện cơ chế xói lở và đề xuất giải pháp công trình kè lát mái kết hợp chân khay đá hộc kiên cố. Đề tài hướng tới mục tiêu ổn định bền vững tuyến đê cấp IV, duy trì hệ số an toàn ổn định công trình đạt mức quy chuẩn từ 1,25 đến 1,30, góp phần ngăn chặn triệt để nguy cơ vỡ đê và thúc đẩy tăng trưởng kinh tế địa phương đạt trên 11,2% 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 lý thuyết tương tác thủy động lực học - địa kỹ thuật công trình, mô phỏng quá trình làm việc của kết cấu bảo vệ bờ qua hệ tương tác 3 môi trường: Nước - Đất - Công trình. Mô hình tính toán phân chia cơ chế tác động thành ba giai đoạn chính: tác động thủy lực biên bên ngoài Pn(y,t) do dòng chảy và sóng, quá trình truyền dẫn áp lực thấm Pt(y,t) qua bề mặt tiếp xúc đất - thân đê, và sự suy giảm ứng suất cắt nội bộ của khối đất khi mực nước hạ đột ngột.

Khung lý thuyết tính toán ổn định bờ đê dựa trên nguyên lý cân bằng giới hạn cơ học đất, sử dụng các phương pháp phân tích trượt cung tròn cải tiến như Bishop, Janbu, Spencer và Morgenstern-Price. Hệ thống lý luận của đề tài bao quát 4 khái niệm cốt lõi:

  • Tuyến chỉnh trị sông và lưu lượng tạo lòng ứng với tần suất thiết kế từ 5% đến 10%.
  • Cơ chế xói sâu lòng dẫn và sự dịch chuyển của lạch sâu áp sát bờ lõm.
  • Hiện tượng xói ngầm, xói trôi hạt mịn qua tầng lọc ngược khi xuất hiện gradien thấm dòng triệt tiêu.
  • Ổn định đẩy nổi và ổn định trượt phẳng của lớp gia cố mái dốc dưới tác động tổng hợp của áp lực thủy tĩnh và tải trọng sóng gió.

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

Nguồn dữ liệu đầu vào được thu thập đồng bộ từ Trạm thủy văn Ba Thá (K46 tả Đáy) với chuỗi số liệu quan trắc liên tục trên 50 năm (từ năm 1965 đến nay). Dữ liệu địa chất công trình được trích xuất từ 5 hố khoan khảo sát hiện trường với chiều sâu từ 2 m đến 16 m, kết hợp dữ liệu đo đạc địa hình đáy sông chi tiết trên 28 mặt cắt ngang đại diện (điển hình như mặt cắt C19 và C28) trải dài dọc đoạn sông nghiên cứu.

Nghiên cứu áp dụng phương pháp chọn mẫu có chủ đích tại 3 vị trí sạt lở nguy hiểm nhất:

  • Vị trí 01 tại K67+700 (chiều dài cung sạt 15 m, ăn sâu 3 m vào bờ, đỉnh sạt +2,79 m).
  • Vị trí 02 tại K67+870 (chiều dài cung sạt 80 m, lấn sâu 3,5 m đến 4,0 m, đỉnh sạt +3,30 m).
  • Vị trí 03 tại K67+950 (chiều dài cung sạt 70 m, lấn sâu 3,0 m, đỉnh sạt +4,11 m).

Về phương pháp phân tích, đề tài sử dụng bộ phần mềm chuyên ngành GEO-SLOPE (mô-đun SEEP/W và SLOPE/W) để mô phỏng dòng thấm bão hòa - không bão hòa và tính toán hệ số ổn định trượt mái dốc. Lý do lựa chọn GEO-SLOPE là nhờ khả năng xử lý bài toán điều kiện biên phức tạp khi mực nước sông rút nhanh từ cao trình lũ +5,00 m về mực nước kiệt +1,94 m và +0,70 m, phản ánh sát thực tế tương tác giữa dòng ngấm và trọng lượng riêng bão hòa của các lớp đất yếu.

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

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

Phân tích bình đồ và trắc ngang địa hình đáy sông cho thấy lòng dẫn đoạn qua xã Vạn Thái bị xâm thực dữ dội, cao độ đáy sông trung bình từ -4,0 m bị đào sâu cục bộ xuống mức -6,0 m, hình thành vách bờ dốc đứng với hệ số mái m < 1,0. Sự biến đổi này tập trung tại bờ lõm của khúc sông cong, nơi dòng chủ lưu áp sát bờ tạo ra lưu tốc dòng chảy mùa lũ vượt ngưỡng 1,5 m/s.

Thống kê diễn biến bùn cát khẳng định tác động rõ rệt từ quy trình vận hành cắt lũ, xả kiệt của Nhà máy thủy điện Hòa Bình từ sau năm 1989. Tổng lượng bùn cát lơ lửng qua sông Hồng hiện giảm 42% (chỉ còn 58% so với thời kỳ 1961-1970 với lưu lượng phù sa 2.487 kg/s), độ đục trung bình giảm còn 60% đến 65% (đạt 261 g/m³). Thành phần bùn cát có xu hướng thô hóa dần với đường kính hạt trung bình dbq = 0,063 mm và hạt lớn nhất dmax = 1,06 mm, biến dòng chảy thành dòng nước đói phù sa làm gia tăng khả năng cuốn trôi lòng sông.

Kết quả kiểm toán ổn định trên mô hình số khi chưa có công trình kè cho thấy sự suy giảm hệ số an toàn nghiêm trọng. Tại mặt cắt C19 và C28, khi mực nước sông duy trì ở mức bình thường +1,94 m, hệ số ổn định chỉ dao động trong khoảng 1,05 đến 1,12. Đặc biệt, trong kịch bản nước rút đột ngột từ cao trình đỉnh lũ +5,00 m về +1,94 m, hệ số an toàn Fs giảm sâu xuống mức 0,82 đến 0,94 (đều nhỏ hơn 1,0), giải thích nguyên nhân sụp đổ hàng loạt các vách đất bãi sông trong các năm 2015-2017.

Thảo luận kết quả

Nguyên nhân gốc rễ gây mất ổn định bờ đê xuất phát từ cấu trúc địa tầng bất lợi gồm lớp đất sét pha lẫn cát dày 2 m đến 16 m nằm trên lớp bùn hữu cơ yếu dày 1,3 m đến 6,0 m. Khi mực nước sông rút nhanh, nước ngầm trong thân đê không kịp thoát ra ngoài, sinh ra gradien thấm vượt quá giới hạn cho phép và làm xuất hiện áp lực kẽ rỗng đẩy nổi các khối đất chân bờ.

Dữ liệu tính toán từ GEO-SLOPE có thể được trực quan hóa thông qua biểu đồ biến thiên hệ số an toàn theo các pha mực nước rút và biểu đồ phân bố đường trượt nguy hiểm nhất đi qua đáy lăng thể đá. So sánh với các giải pháp dân gian truyền thống (như trồng tre chắn sóng hay thả rọ đá thô sơ), các công trình mềm chỉ phát huy hiệu quả khi vận tốc dòng chảy dưới 1,0 m/s. Sau 3 đến 5 năm khai thác, rọ đá thường bị đứt dây thép do va đập bùn cát đáy, còn cọc tre bị mục nát do mực nước dao động liên tục, dẫn đến hiện tượng xói ngầm làm rỗng đất sau lưng kè như từng ghi nhận tại bờ tả sông Đáy đoạn qua thị trấn Vân Đình (K65-K67) và xã Đồng Tiến (K64-K65).

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

  • Xây dựng hoàn chỉnh tuyến kè lát mái kết hợp chân khay đá hộc đổ tự lựa từ K67+000 đến K68+000. Chân kè sử dụng lăng thể đá hộc có đường kính hạt D từ 0,20 m đến 0,35 m nhằm chống xói sâu đáy sông ở cao độ -6,0 m; thân kè lát tấm bê tông đúc sẵn M200 trong khung dầm bê tông cốt thép M250, lót dưới bằng lớp vải địa kỹ thuật lọc ngược chuyên dụng dày 1,5 mm nhằm triệt tiêu áp lực thấm. Mục tiêu nâng hệ số an toàn Fs lên trên 1,25 trong mọi kịch bản lũ rút. Thời gian hoàn thành trong vòng 12 tháng do Ban Quản lý Dự án Đầu tư Xây dựng Công trình Nông nghiệp và Phát triển Nông thôn Hà Nội làm chủ đầu tư.
  • Lắp đặt hệ thống 04 mốc quan trắc chuyển vị ngang và 02 giếng quan trắc áp lực nước ngầm tự động tại bờ tả xã Vạn Thái. Mục tiêu thiết lập cảnh báo sớm dịch chuyển khối trượt khi độ lún vượt quá 5 mm/tháng. Thời gian triển khai trong Quý II/2020 do Chi cục Đê điều và Phòng chống lụt bão Hà Nội chủ trì thực hiện.
  • Quy hoạch chỉnh trị cục bộ lòng dẫn sông Đáy thông qua việc giải tỏa vật cản và nạo vét bãi bồi phía bờ hữu đối diện đoạn K67-K68 nhằm mở rộng mặt cắt ướt thêm 15% đến 20%, kéo dòng chủ lưu dịch chuyển ra xa bờ tả ít nhất 25 m. Thời gian thực hiện từ năm 2020 đến 2022 do Sở Nông nghiệp và Phát triển Nông thôn Hà Nội phối hợp cùng UBND huyện Ứng Hòa điều hành.
  • Ban hành quy chế kiểm soát tải trọng giao thông trên mặt đê tả Đáy (tải trọng tối đa 10 tấn) và cấm tuyệt đối các hoạt động khai thác cát trong phạm vi 1.500 m thượng và hạ lưu tuyến kè. Mục tiêu giảm 100% các xung lực cơ học bất lợi tác động lên mái bờ sông. Thời gian thực hiện liên tục từ năm 2019 do UBND xã Vạn Thái và Công an huyện Ứng Hòa kiểm soát.

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

  • Các cơ quan quản lý nhà nước và Ban Quản lý Dự án (Sở Nông nghiệp và PTNT Hà Nội, Ban QLDA Nông nghiệp Hà Nội): Sử dụng toàn bộ số liệu khảo sát, trắc dọc, trắc ngang và kết quả mô phỏng thủy lực làm luận cứ khoa học để lập báo cáo nghiên cứu khả thi, phê duyệt dự án đầu tư kiên cố hóa đê điều cấp bách trên địa bàn Thủ đô.
  • Đơn vị tư vấn thiết kế và nhà thầu xây dựng công trình thủy lợi: Ứng dụng quy trình tính toán kích thước đá đổ rơi tự do (đường kính D từ 0,15 m đến 0,35 m), cấu tạo tầng lọc ngược bằng vải địa kỹ thuật kết hợp dăm đệm và biện pháp thi công chân kè dưới nước cho các dự án kè sông tương tự.
  • Viện nghiên cứu, giảng viên và học viên cao học chuyên ngành Kỹ thuật công trình thủy, Địa kỹ thuật: Sử dụng làm tài liệu tham khảo điển hình về phương pháp ứng dụng phần mềm GEO-SLOPE trong bài toán tương tác đất - nước - công trình và đánh giá tác động của dòng chảy thô hóa hạ lưu các hồ chứa thủy điện lớn.
  • Chính quyền địa phương xã Vạn Thái và huyện Ứng Hòa: Tham khảo số liệu dự báo xu thế sạt lở tại 3 cung trượt trọng điểm để xây dựng kịch bản ứng phó thiên tai, di dời dân cư khẩn cấp trong mùa mưa bão và quản lý hành lang bảo vệ đê điều an toàn cho 179.000 nhân khẩu.

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

  • Tại sao sạt lở bờ tả sông Đáy đoạn qua xã Vạn Thái diễn biến đặc biệt nghiêm trọng sau năm 2017? Sau các trận bão số 10, 11 và mưa lũ tháng 10 năm 2017, mực nước sông Đáy dâng cao đến +5,00 m rồi rút nhanh xuống mực nước kiệt +0,74 m. Chênh lệch áp lực nước lỗ rỗng lớn trong nền đất yếu (bùn sét dày 1,3 m đến 6,0 m) kết hợp với dòng chủ lưu áp sát bờ lõm đã làm sập đổ 3 cung trượt liên tiếp dài trên 165 m, khoét sâu tới 4,0 m vào bãi đê.

  • Giải pháp chân kè bằng lăng thể đá hộc đổ tự lựa có ưu điểm gì vượt trội so với rọ đá và cọc tre? Lăng thể đá hộc với viên đá đường kính D từ 0,20 m đến 0,35 m có khả năng tự chìm lấp vào các hố xói sâu từ -4,0 m đến -6,0 m, tạo thành chân khay vững chắc chịu được vận tốc dòng chảy mùa lũ trên 1,5 m/s mà không bị mục nát như cọc tre hay đứt gãy lưới thép như rọ đá truyền thống.

  • Phần mềm GEO-SLOPE đóng vai trò như thế nào trong quá trình kiểm toán ổn định công trình? GEO-SLOPE cho phép phân tích 28 mặt cắt ngang điển hình, liên kết giữa mô hình thấm SEEP/W và mô hình trượt SLOPE/W. Phần mềm giúp xác định chính xác mặt trượt cung tròn nguy hiểm nhất và tính toán hệ số an toàn Fs trong điều kiện lũ rút đột ngột từ +5,00 m về +1,94 m, làm cơ sở lựa chọn kích thước hình học thân kè.

  • Vì sao công trình kè lát mái bắt buộc phải bố trí lớp vải địa kỹ thuật kết hợp tầng lọc ngược? Vải địa kỹ thuật dày 1,5 mm có chức năng giữ lại các hạt sét pha và cát mịn có đường kính dưới 0,063 mm không bị dòng thấm cuốn trôi ra sông khi mực nước hạ nhanh. Điều này giúp ngăn ngừa hoàn toàn hiện tượng tạo khoang rỗng dưới mái dốc, tránh gây lún sụt gãy vỡ các tấm bê tông mặt kè.

  • Hoạt động điều tiết của Nhà máy thủy điện Hòa Bình ảnh hưởng như thế nào đến sự ổn định bờ sông Đáy? Thủy điện Hòa Bình giữ lại phần lớn bùn cát, khiến tổng lượng phù sa về hạ lưu giảm 42% và độ đục giảm xuống 261 g/m³. Hiện tượng dòng nước đói bùn cát kết hợp với việc xả - đóng nước phát điện hàng ngày làm mực nước dao động bất thường, gia tăng xói sâu lòng dẫn và thúc đẩy sạt lở chân bờ.

Kết luận

  • Luận văn đã làm rõ bức tranh sạt lở nghiêm trọng tại 3 cung trượt bờ tả sông Đáy (K67+700 đến K67+950) với chiều dài ăn sâu vào bãi sông từ 3,0 m đến 4,0 m, đe dọa an toàn tuyến đê tả Đáy chỉ cách mép lở 38 m.
  • Xác định nguyên nhân gốc rễ gây mất ổn định là do dòng chủ lưu áp sát bờ lõm kết hợp hiện tượng dòng nước đói phù sa (lượng bùn cát giảm 42%) và địa tầng nền chứa lớp bùn sét yếu dày 1,3 m đến 6,0 m.
  • Đề xuất giải pháp kỹ thuật tối ưu gồm kết cấu chân kè lăng thể đá hộc đổ tự lựa (D=0,20 - 0,35 m) kết hợp thân kè lát tấm bê tông đúc sẵn M200 trong khung dầm BTCT và tầng lọc ngược vải địa kỹ thuật, nâng hệ số an toàn Fs từ 0,82 lên trên 1,25.
  • Đóng góp bộ cơ sở dữ liệu thủy văn, địa chất và quy trình ứng dụng GEO-SLOPE hoàn chỉnh, phục vụ trực tiếp cho công tác tư vấn thiết kế và quản lý đê điều tại các khu vực sông cong đồng bằng Bắc Bộ.
  • Lộ trình giai đoạn 2020-2025: Triển khai thi công hoàn chỉnh 1.000 m kè Vạn Thái, thiết lập hệ thống quan trắc chuyển vị tự động và mở rộng nạo vét chỉnh trị dòng chảy, đảm bảo an toàn tuyệt đối cho hành lang đê điều và phát triển kinh tế bền vững cho huyện Ứng Hòa.