Tổng quan nghiên cứu

Hệ thống đê điều tại Việt Nam trải dài hơn 13.200 km, bao gồm khoảng 10.600 km đê sông và hơn 2.600 km đê biển, giữ vai trò đặc biệt quan trọng trong chiến lược bảo vệ tính mạng nhân dân và cơ sở hạ tầng kinh tế xã hội. Riêng hệ thống đê sông chính từ cấp III trở lên đã vượt quá 2.500 km. Tuy nhiên, theo số liệu từ các cơ quan khí tượng thủy văn, thiên tai bão lũ trong giai đoạn 10 năm từ 1997 đến 2006 đã làm chết và mất tích gần 7.500 người, gây thiệt hại vật chất ước tính chiếm khoảng 1,5% GDP mỗi năm. Thực tế cho thấy phương pháp thiết kế truyền thống theo hệ số an toàn tất định bộc lộ nhiều bất cập khi đối mặt với tác động bất thường của biến đổi khí hậu và hiện tượng nước biển dâng.

Điển hình tại huyện Giao Thủy, tỉnh Nam Định, hệ thống phòng lũ bao gồm 31,161 km đê biển trải dài từ K0+000 đến K31+161 và 11,702 km đê hữu sông Hồng từ K208+000 đến K219+702. Toàn huyện có diện tích tự nhiên 23.799,64 ha với quy mô dân số 205.075 người (theo thống kê năm 2010), nằm kẹp giữa cửa Ba Lạt của sông Hồng và cửa Hà Lạn của sông Sò. Tuyến đê biển tại đây có lịch sử hình thành hơn 250 năm nhưng chỉ được thiết kế để chống chịu bão cấp 9 và triều trung bình tần suất 5%. Khi cơn bão số 7 năm 2005 đổ bộ với cường độ cấp 10 đến 12 kết hợp triều cường, đê biển đã bị vỡ và hư hại nghiêm trọng.

Nghiên cứu này được thực hiện nhằm xây dựng phương pháp đánh giá mức độ an toàn và tối ưu hóa tiêu chuẩn thiết kế công trình phòng chống lũ dựa trên lý thuyết độ tin cậy và phân tích rủi ro. Thông qua việc lượng hóa xác suất xảy ra sự cố và đường cong thiệt hại, đề tài cung cấp cơ sở khoa học để xác lập cao trình đỉnh đê, giảm thiểu rủi ro kinh tế và bảo vệ an toàn tính mạng cho cộng đồng dân cư ven biển.

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 độ tin cậy kết cấu hiện đại kết hợp với khung lý thuyết quản lý rủi ro thiên tai tích hợp. Trạng thái an toàn của công trình được mô tả thông qua hàm độ tin cậy tổng quát Z = R - S, trong đó R đại diện cho sức kháng phá hoại của kết cấu và S là tải trọng tác động từ môi trường tự nhiên. Sự cố phá hoại xảy ra khi giá trị Z nhỏ hơn hoặc bằng 0, tương ứng với xác suất xảy ra sự cố Pf. Độ tin cậy của công trình được xác định bằng xác suất để Z lớn hơn 0.

Lý thuyết phân tích rủi ro định nghĩa rủi ro tổng thể bằng tích số giữa xác suất xảy ra ngập lụt và hậu quả thiệt hại kinh tế xã hội. Hệ thống kiểm soát rủi ro kết hợp hai nhóm giải pháp: nhóm công trình nhằm giảm thiểu xác suất vỡ đê và nhóm phi công trình nhằm giảm nhẹ hậu quả ngập lụt.

Sáu cơ chế phá hoại cơ bản được tích hợp vào mô hình cây sự cố bao gồm:

  • Cơ chế sóng tràn và chảy tràn qua đỉnh đê khi mực nước vượt quá cao trình đỉnh đê Zc.
  • Cơ chế mất ổn định kết cấu bảo vệ mái kè dưới áp lực động của sóng biển theo công thức không thứ nguyên của Van der Meer.
  • Cơ chế xói lở chân đê dựa trên mô hình động lực học của Sumer và Fredsoe.
  • Cơ chế xói ngầm và đẩy trồi nền đất cát pha dưới thân đê theo tiêu chuẩn Blight và Sellmeijer.
  • Cơ chế mất ổn định trượt mái dốc tổng thể của đê phía biển và đê phía đồng theo phương pháp mặt trượt trụ tròn Bishop.

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

Nghiên cứu áp dụng phương pháp phân tích ngẫu nhiên cấp độ II (Level II Reliability Analysis). Lý do lựa chọn phương pháp này là vì dữ liệu khí tượng thủy văn vùng ven biển có tính bất định cao, trong khi phương pháp cấp độ III (như tích phân đa lớp Monte Carlo) đòi hỏi nguồn tài nguyên tính toán rất lớn và chuỗi số liệu ngẫu nhiên đầy đủ. Bằng cách sử dụng phương pháp khai triển chuỗi Taylor quanh điểm thiết kế kết hợp thuật toán tối ưu Hasofer-Lind và kỹ thuật chuyển đổi biến phi chuẩn Rackwitz-Fiessler, chỉ số độ tin cậy beta và xác suất sự cố được xác định nhanh chóng với độ chính xác cao.

Cỡ mẫu nghiên cứu bao gồm chuỗi số liệu quan trắc mực nước, chiều cao sóng, áp lực gió và các trận bão lịch sử được thống kê từ năm 1930 đến năm 2012. Phương pháp chọn mẫu toàn diện được áp dụng đối với chuỗi số liệu thủy văn tại trạm Văn Lý và Giao Xuân, xử lý phân phối xác suất thông qua phần mềm BestFit.

Mô hình thủy lực hai chiều MIKE 21 với lưới tính phi cấu trúc được thiết lập để mô phỏng trường dòng chảy và lan truyền ngập lụt trên toàn bộ diện tích 23.799,64 ha của huyện Giao Thủy. Mô hình được kiểm định bằng chuỗi số liệu thực đo liên tục trong 10 ngày từ ngày 18 đến ngày 28 tháng 7 năm 2009, bảo đảm độ tin cậy của biên thủy động lực học với sai số tính toán dưới 5%.

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

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

Kết quả phân tích độ tin cậy cho thấy hệ thống đê biển Giao Thủy hiện tại tiềm ẩn nguy cơ mất an toàn nghiêm trọng trước các tác động tổ hợp của bão và triều cường. Tuyến đê biển dài 31,161 km có nhiều đoạn xung yếu; cụ thể đoạn từ K0 đến K12+600 chỉ đạt cao trình +3,20 m, thiếu hụt từ 0,80 m đến 1,10 m so với yêu cầu phòng chống sóng tràn an toàn trong điều kiện bão cấp 9.

Trong 6 cơ chế phá hoại được xem xét, cơ chế sóng tràn qua đỉnh và cơ chế bóc dỡ tấm lát bảo vệ mái đê có xác suất xảy ra sự cố cao nhất. Khi gặp sóng biển dâng cao kết hợp gió giật, áp lực sóng âm khi rút nước làm dão và bật các khối đá lát khan 0,35 m. Phân tích ổn định trượt mái dốc xác định hệ số an toàn nhỏ nhất của mái đê phía biển là 1,501 và mái đê phía đồng là 1,805.

Mô phỏng ngập lụt bằng MIKE 21 trên địa bàn huyện Giao Thủy đưa ra các thông số cụ thể:

  • Kịch bản vỡ đê hữu sông Hồng: Khi mực nước sông đạt đỉnh lũ thiết kế, vết vỡ đê sau 10 giờ làm ngập sâu 6 xã ven đê với độ sâu ngập phổ biến từ 1,2 m đến 2,8 m, tốc độ dòng chảy lan truyền đạt từ 1,5 m/s đến 2,2 m/s.
  • Kịch bản vỡ đê biển: Nước biển tràn sâu vào nội đồng, làm ngập hơn 60% diện tích đất nông nghiệp và đầm nuôi trồng thủy sản, đe dọa trực tiếp đến tính mạng của 205.075 cư dân địa phương.

Khi nâng cao trình đỉnh đê từ +3,20 m lên mức +4,50 m đến +5,20 m, xác suất vỡ đê tổng thể giảm từ mức 0,085 xuống dưới 0,001 mỗi năm. Phân tích kinh tế cho thấy chi phí đầu tư nâng cấp công trình tăng thêm 35% đến 45% nhưng giúp cắt giảm hơn 70% rủi ro thiệt hại kinh tế kỳ vọng quy đổi.

Thảo luận kết quả

Nguyên nhân cốt lõi khiến đê biển Giao Thủy dễ bị phá hủy là do cấu tạo địa chất thân đê và nền đê chủ yếu là cát mịn bọc đất thịt hình thành trên vùng bồi tích trẻ. Dưới chênh lệch cột nước áp lực giữa phía biển và phía đồng, hiện tượng xói ngầm và đẩy trồi diễn ra nhanh chóng, làm rỗng ruột thân đê trước khi đỉnh đê bị tràn.

So sánh với phương pháp tính toán hệ số an toàn truyền thống vốn coi tải trọng là đại lượng cố định, phương pháp ngẫu nhiên đã làm rõ tác động phân tán của từng biến số. Chiều cao sóng và mực nước triều đóng góp tới hơn 58% vào tổng xác suất phá hỏng kết cấu mái kè. Kết quả này tương đồng với các công bố của dự án FLOODsite tại châu Âu và chương trình đồng bằng Deltaplan của Hà Lan, khẳng định tính tất yếu của việc chuyển dịch sang tiêu chuẩn thiết kế theo độ tin cậy.

Dữ liệu nghiên cứu được biểu diễn trực quan qua hệ thống bản đồ phân vùng ngập lụt hai chiều và đồ thị đường cong quan hệ giữa tần suất phòng lũ với tổng chi phí hệ thống. Sơ đồ cây sự cố giúp các nhà kỹ thuật dễ dàng nhận diện khâu xung yếu nhất của tuyến đê để ưu tiên phân bổ ngân sách gia cố.

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

  1. Nâng cấp toàn diện trắc dọc và cao trình đê biển: Tôn cao 31,161 km tuyến đê biển Giao Thủy từ cao độ hiện hữu +3,20 m lên mức thiết kế tối ưu từ +4,50 m đến +5,20 m, mở rộng mặt đê lên 4,0 m đến 5,0 m kết hợp làm đường cứu hộ giao thông. Mục tiêu hạ xác suất sự cố tổng thể xuống dưới 0,001/năm. Thời gian hoàn thành trong giai đoạn 2014-2020 do Ủy ban nhân dân tỉnh Nam Định và Bộ Nông nghiệp và Phát triển nông thôn phối hợp chỉ đạo.
  2. Cải tiến kết cấu áo kè và gia cố chân khay chống xói: Thay thế toàn bộ lớp đá lát khan truyền thống bằng cấu kiện bê tông đúc sẵn tự chèn dày 0,35 m trên lớp vải địa kỹ thuật lọc cát. Xây dựng hệ thống chân khay bằng ống buy bê tông cốt thép mác 200, đường kính 100 cm, ngập sâu 1,5 m chứa đá hộc để chống xói chân đê dưới tác động sóng bão cấp 11 đến 12. Tiến độ triển khai từ năm 2014 đến 2018 do Chi cục Đê điều và Phòng chống lụt bão Nam Định thực hiện.
  3. Thiết lập hệ thống cảnh báo sớm và phương án sơ tán dân số: Tích hợp bản đồ ngập lụt số hóa từ mô hình MIKE 21 vào hệ thống quản lý rủi ro thiên tai huyện Giao Thủy, xây dựng kịch bản sơ tán chi tiết cho 205.075 người dân tại 22 xã, thị trấn khi xảy ra sự cố vỡ đê. Mục tiêu giảm thiểu 95% nguy cơ thương vong về người. Thời gian triển khai giai đoạn 2015-2017 do Ban Chỉ huy Phòng chống thiên tai địa phương chủ trì.
  4. Ban hành quy chuẩn thiết kế công trình phòng lũ theo độ tin cậy: Xây dựng tiêu chuẩn kỹ thuật quốc gia về đánh giá an toàn đê điều theo chỉ số độ tin cậy beta thay thế dần phương pháp hệ số an toàn đơn nhất cho hơn 2.500 km đê từ cấp III trở lên. Lộ trình thực hiện 2015-2025 do Bộ Khoa học và Công nghệ cùng Trường Đại học Thủy lợi phối hợp thực hiện.

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

  1. Kỹ sư thiết kế và tư vấn xây dựng công trình thủy lợi: Luận văn cung cấp công thức chi tiết, quy trình giải thuật phân tích độ tin cậy Cấp độ II và phương pháp thiết lập trạng thái giới hạn cho 6 cơ chế phá hoại đê sông, đê biển.
  2. Cán bộ quản lý quy hoạch và đê điều tại các tỉnh duyên hải: Cung cấp cơ sở khoa học để đánh giá mức độ rủi ro, phân bổ vốn đầu tư nâng cấp cho 31,161 km đê biển Nam Định cũng như các hệ thống đê lân cận tại Thái Bình, Hải Phòng.
  3. Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Thủy công và Kỹ thuật Biển: Là tài liệu học thuật giá trị về mô hình hóa ngập lụt bằng MIKE 21, phân tích thống kê BestFit và ứng dụng lý thuyết ngẫu nhiên trong công trình thủy.
  4. Các tổ chức quốc tế và cơ quan tài trợ dự án thích ứng biến đổi khí hậu: Sử dụng phương pháp luận tối ưu hóa kinh tế và rủi ro cá nhân để thẩm định tính khả thi và hiệu quả của các dự án ODA nâng cấp đê điều phòng chống thiên tai tại Việt Nam.

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

  1. Điểm khác biệt lớn nhất giữa thiết kế đê theo độ tin cậy và thiết kế truyền thống là gì? Phương pháp truyền thống dựa trên các đại lượng tải trọng và độ bền cố định cùng một hệ số an toàn đơn nhất, chưa phản ánh được sự biến thiên của tự nhiên. Phương pháp theo độ tin cậy xem xét toàn bộ các yếu tố sóng, gió, địa chất dưới dạng các biến ngẫu nhiên có quy luật phân bố xác suất, giúp xác định chính xác xác suất xảy ra sự cố của từng bộ phận kết cấu.

  2. Tại sao hệ thống đê biển Giao Thủy dễ bị tổn thương trước các cơn bão mạnh? Tuyến đê biển Giao Thủy dài 31,161 km chủ yếu đắp bằng đất cát pha trên nền trầm tích bở rời hơn 250 năm. Cao trình đê tại nhiều vị trí chỉ đạt +3,20 m và mái đê gia cố bằng đá lát khan thô sơ. Khi bão vượt cấp 9 kết hợp triều cường, sóng tràn qua đỉnh gây xói lở mặt trong và áp lực sóng bóc tung mái kè.

  3. Mô hình MIKE 21 đóng góp vai trò gì trong việc đánh giá an toàn đê điều? Mô hình MIKE 21 tính toán dòng chảy hai chiều trên lưới phi cấu trúc phủ khắp 23.799,64 ha của huyện Giao Thủy. Mô hình giúp mô phỏng chi tiết vận tốc dòng nước từ 1,5 m/s đến 2,2 m/s và độ sâu ngập lụt từ 1,2 m đến 2,8 m theo thời gian thực khi xảy ra sự cố vỡ đê, hỗ trợ lập bản đồ rủi ro.

  4. Tiêu chuẩn an toàn tối ưu theo quan điểm kinh tế được tính toán như thế nào? Tiêu chuẩn tối ưu đạt được tại điểm mà tổng chi phí toàn vòng đời công trình là nhỏ nhất. Tổng chi phí này bao gồm chi phí đầu tư xây dựng nâng cấp mặt cắt đê, chi phí duy tu bảo dưỡng hàng năm quy về hiện tại, và giá trị thiệt hại kinh tế do lũ lụt gây ra nhân với xác suất xảy ra sự cố ngập lụt tương ứng.

  5. Luận văn này có thể nhân rộng cho các hệ thống đê điều khác tại Việt Nam không? Hoàn toàn có thể mở rộng. Bằng cách thu thập chuỗi số liệu khí tượng thủy văn địa phương, phân tích thống kê bằng BestFit và thiết lập mô hình thủy lực 2 chiều, quy trình phân tích độ tin cậy Cấp độ II hoàn toàn áp dụng được cho hơn 13.200 km đê sông và đê biển trên phạm vi cả nước.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp luận đánh giá an toàn hệ thống đê vùng cửa sông ven biển dựa trên lý thuyết độ tin cậy kết cấu Cấp độ II và phân tích rủi ro.
  • Lượng hóa chính xác xác suất hư hỏng của 6 cơ chế phá hoại đối với 31,161 km đê biển và 11,702 km đê hữu Hồng thuộc huyện Giao Thủy, tỉnh Nam Định.
  • Thiết lập mô hình ngập lụt MIKE 21 trên diện tích 23.799,64 ha, cung cấp dữ liệu thủy lực chi tiết phục vụ bảo vệ an toàn cho 205.075 người dân.
  • Xác định điểm cân bằng tối ưu giữa chi phí đầu tư nâng cấp công trình và mức độ giảm thiểu rủi ro thiệt hại kinh tế xã hội.
  • Đề xuất 4 nhóm giải pháp kỹ thuật công trình và quản lý phi công trình đồng bộ, khả thi cho khu vực ven biển Đồng bằng Bắc Bộ.

Đóng góp quan trọng nhất của công trình là chuyển đổi phương thức đánh giá an toàn công trình phòng chống lũ tại Việt Nam từ kinh nghiệm tất định sang định lượng ngẫu nhiên khoa học. Trong giai đoạn 2015-2025, hướng nghiên cứu cần tiếp tục mở rộng mô phỏng không gian ba chiều và thuật toán Monte Carlo cho các hệ thống đê phức tạp. Các cơ quan quản lý và viện nghiên cứu chuyên ngành cần sớm tham khảo và áp dụng kết quả luận văn vào việc cập nhật quy chuẩn thiết kế đê điều quốc gia.