Tổng quan nghiên cứu

Năm 2005 được ghi nhận là một mốc thời gian lịch sử về thảm họa thiên tai tại Việt Nam khi 8 cơn bão nhiệt đới liên tiếp đổ bộ vào các dải ven biển. Đặc biệt, cơn bão số 7 (bão Damrey) với sức gió mạnh cấp 12 (118 đến 133 km/h) kết hợp cùng triều cường dâng cao đã làm vỡ hơn 28 km đê biển tại các tỉnh Nam Định và Thanh Hóa, gây ngập lụt trên 130.000 ha đất canh tác lúa và làm sập đổ 1.194 ngôi nhà. Thiệt hại kinh tế ước tính trực tiếp thời điểm đó vượt mức 50 triệu USD, gióng lên hồi chuông cảnh báo cấp thiết về mức độ an toàn của hệ thống công trình phòng chống lũ lụt ven bờ.

Thực trạng quản lý công trình đê biển tại Việt Nam trong nhiều thập kỷ qua chủ yếu dựa vào phương pháp tiền định truyền thống với chu kỳ lặp thiết kế cố định chỉ từ 1/20 năm hoặc ngắn hơn. Cách tiếp cận này bộc lộ những lỗ hổng nghiêm trọng do đánh giá thấp tải trọng thủy lực cực đoan và chưa tính toán đầy đủ các yếu tố bất định của môi trường biển. Mục tiêu cốt lõi của nghiên cứu là đánh giá toàn diện mức độ an toàn và độ tin cậy thực tế của hệ thống đê biển hiện hữu, phát triển khung lý thuyết thiết kế theo xác suất và quản lý rủi ro ngập lụt ven biển, từ đó đề xuất các tiêu chuẩn an toàn mới cùng giải pháp cải tạo hệ thống đê điều một cách bền vững.

Phạm vi nghiên cứu tập trung vào vùng duyên hải Đồng bằng sông Hồng, lấy hệ thống đê biển tỉnh Nam Định với chiều dài gần 70 km làm khu vực điển hình. Không gian nghiên cứu được giới hạn tự nhiên bởi 4 cửa sông lớn gồm Ba Lạt, Hà Lạn, Lạch Giang và Đáy, dựa trên chuỗi số liệu quan trắc thủy văn và biến động đường bờ từ năm 1970 đến năm 2005. Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc chuyển đổi mô hình thiết kế từ tiền định sang xác suất tin cậy, giúp nâng chuẩn bảo vệ an toàn kinh tế từ mức 1/20 năm hiện hành lên 1/50 năm và 1/100 năm trong tương lai, đồng thời thiết lập ngưỡng an toàn rủi ro xã hội ở mức 1/150 năm nhằm bảo vệ tính mạng cho hơn 1 triệu người dân vùng duyên hải.

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 sự kết hợp chặt chẽ giữa các nền tảng lý thuyết kỹ thuật công trình thủy và kinh tế học rủi ro:

Thứ nhất, lý thuyết độ tin cậy kết cấu công trình (Structural Reliability Theory) được áp dụng để mô hình hóa trạng thái an toàn của từng phân đoạn đê. Cốt lõi của lý thuyết này là hàm trạng thái giới hạn (Limit State Function) biểu diễn mối quan hệ giữa sức kháng của công trình và tải trọng thủy lực môi trường tác động. Khác với hệ số an toàn tĩnh, lý thuyết này xem xét toàn bộ các tham số đầu vào như chiều cao sóng, mực nước triều dâng, chất lượng khối phủ kè và độ bền của nền đất dưới dạng các biến ngẫu nhiên mang phân phối xác suất.

Thứ hai, khung đánh giá rủi ro ngập lụt tích hợp (Integrated Flood Risk Framework) kế thừa và phát triển từ mô hình của Vrijling, kết hợp với mô hình tối ưu hóa chi phí - lợi ích kinh tế mở rộng từ nguyên lý van Dantzig. Khung lý thuyết này định nghĩa rủi ro là tích số giữa xác suất xảy ra biến cố ngập lụt và tổng thiệt hại phát sinh. Trong đó, tổn thất nhân mạng được lượng hóa qua hai chỉ số chính: rủi ro cá nhân (xác suất tử vong của một người tại vị trí nhất định trong một năm) và rủi ro xã hội (xác suất xảy ra thảm họa với số người thiệt mạng vượt quá mức quy định, biểu diễn qua đường cong FN).

Thứ ba, mô hình thủy động lực học và hình thái học bờ biển hai chiều (UNIBEST) tích hợp các công thức vận chuyển bùn cát dọc bờ của Bijker và Van Rijn nhằm đánh giá quá trình tương tác giữa sóng biển, dòng chảy triều và sự suy giảm bãi trước đê.

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

Nguồn dữ liệu sử dụng trong nghiên cứu bao gồm chuỗi số liệu quan trắc mực nước biển và sóng ven bờ liên tục trong 35 năm (1970 - 2005) tại các trạm hải văn ven bờ phía Bắc, dữ liệu đo đạc mặt cắt bãi biển định kỳ 1 đến 2 năm/lần từ năm 1972 đến năm 2003 tại bờ biển Hải Hậu, cùng cơ sở dữ liệu thống kê thiệt hại do bão lũ trong hơn 30 năm do Cục Quản lý Đê điều và Phòng chống lụt bão thu thập.

Phương pháp chọn mẫu cực trị vượt ngưỡng (Peak-Over-Threshold) kết hợp với phương pháp lấy mẫu giá trị cực đại hàng năm (Annual Maxima) được lựa chọn để trích xuất các biến cố bão lũ lớn nhất. Nghiên cứu áp dụng phương pháp Phân tích tần suất vùng (Regional Frequency Analysis - RFA) và Khớp hàm phân phối hỗn hợp (Mixed Distribution Fitting) để ngoại suy các mức tải trọng thủy lực cực đoan ứng với chu kỳ lặp từ 100 đến 10.000 năm. Lý do lựa chọn RFA là nhằm khắc phục hạn chế về số lượng trạm đo và độ dài chuỗi dữ liệu tại từng vị trí đơn lẻ; việc tích hợp dữ liệu từ các trạm lân cận có cùng đặc trưng thủy văn giúp tăng bậc tự do thống kê và giảm thiểu độ bất định của các giá trị dự báo.

Quy trình phân tích độ tin cậy được thực hiện thông qua mô hình cây sự cố (Fault Tree Analysis) và phân tích độ nhạy của các biến ngẫu nhiên. Toàn bộ các mô hình phân tích và thuật toán tối ưu hóa được hoàn thiện và kiểm chứng trong giai đoạn nghiên cứu từ năm 2006 đến năm 2010.

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

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

Quá trình phân tích độ tin cậy và rủi ro đối với hệ thống đê biển Nam Định đã mang lại các phát hiện cốt lõi:

Thứ nhất, mức độ an toàn thực tế của tuyến đê biển thấp hơn rất nhiều so với tiêu chuẩn thiết kế quy định. Mặc dù tiêu chuẩn hiện hành yêu cầu tuyến đê phải chịu được chu kỳ bão 1/20 năm (xác suất hư hỏng 0,05/năm), kết quả tính toán độ tin cậy hệ thống cho thấy xác suất phá hủy thực tế lên tới 0,15/năm. Điều này chỉ ra rằng công trình có nguy cơ bị phá vỡ trung bình 7 năm một lần, trùng khớp với tần suất hư hỏng thực tế được quan sát trong giai đoạn 1992 - 2006.

Thứ hai, cơ cấu đóng góp của các cơ chế phá hủy đê được định lượng rõ ràng: Hiện tượng sóng tràn qua đỉnh đê gây xói lở mái trong là cơ chế phá hủy áp đảo nhất, đóng góp xấp xỉ 44% vào tổng xác suất hư hỏng của toàn hệ thống; mất ổn định kết cấu các khối phủ mái ngoài (revetment blocks) đóng góp 25%; mất ổn định địa kỹ thuật trượt mái trong đóng góp 20%; hiện tượng nước dâng chảy tràn qua đỉnh đê đóng góp 10%. Các cơ chế hư hỏng khác như thấm ngầm (piping) và mất ổn định chân đê chỉ đóng góp dưới 1%.

Thứ ba, kích thước hình học hiện hữu của đê biển hoàn toàn không đáp ứng yêu cầu kỹ thuật. Để đảm bảo đúng tiêu chuẩn an toàn 1/20 năm theo phương pháp thiết kế tin cậy, cao trình đỉnh đê bắt buộc phải đạt mức 6,80 mét, trong khi cao trình thực tế của tuyến đê Nam Định chỉ đạt 5,50 mét (thiếu hụt 1,30 mét).

Thứ tư, nghiên cứu phát hiện tốc độ dâng của mực nước biển tại dải ven biển Nam Định đạt khoảng 0,68 cm/năm, đi kèm 3 chu kỳ dao động thủy văn đặc thù: chu kỳ 6 tháng (tác động của hai mùa gió mùa), 12 tháng (mùa bão hàng năm) và 224 tháng (xấp xỉ 18,6 năm, tương ứng với một chu kỳ triều thiên văn). Cùng với đó, tốc độ xói lở bờ biển tại khu vực Hải Hậu diễn biến nghiêm trọng với mức thoái lui trung bình từ 10 đến 20 mét/năm, cá biệt tại xã Hải Lý tốc độ xói lở lên tới 22 đến 30 mét/năm.

Thảo luận kết quả

Kết quả nghiên cứu làm sáng tỏ nguyên nhân sâu xa của tình trạng vỡ đê liên tục trong các mùa mưa bão. Phương pháp thiết kế tiền định cũ chỉ sử dụng các giá trị tải trọng tĩnh và bỏ qua sự xói sâu của bãi trước đê (trung bình từ 0,2 đến 0,6 mét/năm). Khi bãi nông phía trước bị hạ thấp, các đợt sóng lớn ngoài khơi với chiều cao từ 3,0 đến 6,0 mét dễ dàng truyền thẳng vào chân đê mà không bị tiêu hao năng lượng, gây nên lưu lượng sóng tràn thực tế lên tới 190 lít/giây/mét, vượt xa ngưỡng cho phép thông thường từ 1 đến 10 lít/giây/mét.

Trong phân tích rủi ro, dữ liệu tử vong do lũ lụt có thể được trực quan hóa thông qua biểu đồ đường cong FN. Đường cong này phản ánh rõ nét mức độ rủi ro xã hội do ngập lụt tại Việt Nam đang cao gấp 5 đến 10 lần so với các quốc gia phát triển như Hà Lan. Bên cạnh đó, việc biểu diễn phân bố tổn thất kinh tế qua hàm phân phối log-normal cho thấy giá trị thiệt hại kỳ vọng hàng năm do ngập lụt tại Việt Nam vào khoảng 589 triệu USD, tương đương mức tổn thất trung bình từ 1,0% đến 1,5% GDP quốc gia mỗi năm trong suốt một thập kỷ.

Khi so sánh với mô hình kinh tế van Dantzig truyền thống vốn chỉ áp dụng cho địa hình lòng chảo bằng phẳng, việc mở rộng công thức cho các vùng đất dốc ven biển đã chứng minh tính ưu việt. Mô hình kinh tế chỉ ra rằng mức an toàn tối ưu cho đê biển Nam Định trong điều kiện hiện tại phải là 1/50 năm và cần nâng lên mức 1/100 năm khi quy mô tích lũy tài sản kinh tế gia tăng trong tương lai.

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

Nhằm giải quyết triệt để các hạn chế của hệ thống phòng hộ ven biển, nghiên cứu đưa ra 4 nhóm giải pháp chiến lược:

  1. Nâng cấp cao trình và gia cố kết cấu bề mặt đê biển: Bộ Nông nghiệp và Phát triển Nông thôn cùng Sở Nông nghiệp các tỉnh ven biển cần thực hiện điều chỉnh cao trình đỉnh đê từ mức 5,50 mét hiện hữu lên tối thiểu 6,80 mét tại các vị trí xung yếu. Đồng thời, cần tăng cường khả năng chống xói của đỉnh đê và mái trong bằng các khối phủ liên kết bền vững, chịu được lưu lượng sóng tràn trên 10 lít/giây/mét, hoàn thành việc nâng cấp trong giai đoạn 2026 - 2030.

  2. Triển khai giải pháp bảo vệ bãi biển tích hợp: Viện Khoa học Thủy lợi Việt Nam cùng chính quyền địa phương cần kết hợp đê biển với các công trình mỏ hàn (groins), đê ngầm giảm sóng ngoài khơi và công nghệ nuôi bãi cát nhân tạo (sand nourishment). Mục tiêu là duy trì độ nông của bãi trước đê, giảm tốc độ xói lở đường bờ từ mức 20 mét/năm xuống dưới 3 mét/năm, thực hiện thí điểm trong vòng 3 năm tới trước khi nhân rộng toàn dải Đồng bằng sông Hồng.

  3. Hoàn thiện quy chuẩn kỹ thuật thiết kế đê biển theo xác suất: Cục Quản lý Đê điều và Phòng chống Thiên tai cần chủ trì chuyển đổi khung tiêu chuẩn thiết kế quốc gia từ phương pháp tiền định sang phương pháp thiết kế theo độ tin cậy. Mục tiêu là nâng chuẩn thiết kế kinh tế cho các vùng đê bao ven biển lên 1/50 năm ở giai đoạn trước mắt và 1/100 năm vào năm 2035, đồng thời đảm bảo chuẩn an toàn rủi ro xã hội đạt mức 1/150 năm.

  4. Hiện đại hóa hệ thống quan trắc và cảnh báo sớm thủy văn hải văn: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần chuẩn hóa việc ứng dụng phương pháp phân tích tần suất vùng RFA và mô hình mô phỏng bão số trị. Mục tiêu là nâng cao độ chính xác dự báo nước dâng do bão và sóng cực trị trước năm 2028, cung cấp chuỗi số liệu đầu vào chuẩn xác cho công tác vận hành công trình.

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

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

  1. Cơ quan quản lý nhà nước và cơ quan lập chính sách (Bộ Nông nghiệp và Phát triển Nông thôn, Ban Chỉ đạo Quốc gia về Phòng chống thiên tai): Tài liệu cung cấp cơ sở khoa học định lượng về rủi ro xã hội và kinh tế, hỗ trợ quá trình phân bổ ngân sách đầu tư công trình bảo vệ bờ biển một cách tối ưu, tránh lãng phí nguồn lực quốc gia.

  2. Kỹ sư tư vấn thiết kế và các viện nghiên cứu chuyên ngành công trình thủy lợi, công trình biển: Nắm bắt phương pháp tính toán độ tin cậy, ứng dụng các hàm trạng thái giới hạn để tính toán hình học mặt cắt đê, kích thước khối kè và kiểm soát rủi ro từ các cơ chế phá hủy phức tạp.

  3. Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Công trình Thủy, Quản lý Bờ biển và Biến đổi khí hậu: Sử dụng tài liệu như một giáo trình tham khảo chuyên sâu về phân tích thống kê cực trị vùng (RFA), mô hình hóa toán học các biến ngẫu nhiên và lý thuyết thiết kế theo xác suất tiên tiến.

  4. Các tổ chức quốc tế và các nhà quản lý dự án hỗ trợ kỹ thuật (UNDP, Ngân hàng Thế giới, JICA): Tham khảo mô hình phân tích kinh tế rủi ro thiên tai để lập báo cáo khả thi cho các dự án tăng cường năng lực thích ứng với biến đổi khí hậu tại các vùng đồng bằng ngập lũ ở các nước đang phát triển.

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

  1. Phương pháp thiết kế đê biển theo xác suất có ưu điểm gì vượt trội so với phương pháp tiền định truyền thống? Phương pháp thiết kế theo xác suất lượng hóa trực tiếp các yếu tố ngẫu nhiên như sóng, nước dâng và sức kháng vật liệu dưới dạng các biến phân phối thống kê. Thay vì chỉ dựa vào hệ số an toàn tĩnh mang tính cảm tính, phương pháp này cho phép xác định chính xác xác suất phá hủy thực tế của toàn bộ tuyến đê theo chiều dài công trình, giúp tối ưu hóa chi phí đầu tư và nâng cao độ bền vững.

  2. Tại sao cơ chế tràn sóng lại là nguyên nhân phá hủy đê biển phổ biến nhất tại vùng duyên hải Đồng bằng sông Hồng? Thực tế các trận bão lớn như bão Damrey năm 2005 cho thấy bão kết hợp triều cường tạo nên sóng lớn với lưu lượng tràn lên tới 190 lít/giây/mét. Do đỉnh đê hiện trạng chỉ cao 5,50 mét và mái trong chưa được gia cố chống xói, dòng chảy tràn với động năng mạnh nhanh chóng làm xói rỗng thân đê từ phía sau, chiếm tới 44% tổng xác suất hư hỏng công trình.

  3. Quá trình xói lở bờ biển tại khu vực Hải Hậu (Nam Định) đe dọa độ an toàn của đê biển như thế nào? Bờ biển Hải Hậu bị thoái lui từ 10 đến 20 mét mỗi năm và xói sâu bãi trước đê từ 0,2 đến 0,6 mét/năm do thiếu hụt bùn cát từ hệ thống sông Hồng. Bãi trước bị hạ thấp làm mất tác dụng giảm sóng tự nhiên, khiến sóng lớn ngoài khơi tác động trực tiếp vào chân đê, gây mất ổn định lớp kè ngoài và tăng gấp đôi nguy cơ vỡ đê.

  4. Chuẩn an toàn rủi ro xã hội 1/150 năm cho đê biển Nam Định được xác lập dựa trên nguyên lý nào? Chuẩn này được xác lập dựa trên khung đánh giá rủi ro của Vrijling, kết hợp giữa tỷ lệ tử vong nền do các nguyên nhân tự nhiên với mức độ chấp nhận rủi ro của cộng đồng. Khi phân bổ rủi ro quốc gia cho từng vùng theo đường cong FN, vùng trũng thấp Nam Định bắt buộc phải đạt mức bảo vệ 1/150 năm để đảm bảo an toàn tính mạng cho người dân.

  5. Mô hình kinh tế van Dantzig mở rộng mang lại lợi ích gì trong việc lựa chọn chiều cao đê tối ưu? Mô hình mở rộng giúp cân bằng giữa tổng kinh phí đầu tư xây dựng đê với giá trị thiệt hại kinh tế kỳ vọng trong suốt vòng đời công trình tại các vùng đồng bằng dốc. Tại Nam Định, mô hình chỉ ra rằng việc nâng chuẩn an toàn lên 1/50 năm và tiến tới 1/100 năm sẽ giảm thiểu mức thiệt hại hàng năm (vốn đang chiếm từ 1,0% đến 1,5% GDP quốc gia).

Kết luận

  • Khẳng định tính tất yếu của việc chuyển đổi phương pháp thiết kế công trình bảo vệ bờ bờ biển từ tiền định sang mô hình độ tin cậy và xác suất.
  • Xác định rõ mức độ an toàn thực tế của đê biển Nam Định chỉ đạt chu kỳ 7 năm/lần (xác suất 0,15/năm) và cơ chế sóng tràn là nguyên nhân hư hỏng chính chiếm 44%.
  • Xây dựng thành công bộ công cụ phân tích cực trị thủy lực tích hợp phương pháp phân tích tần suất vùng RFA cùng các chu kỳ dao động tự nhiên 6, 12 và 224 tháng của mực nước biển.
  • Đề xuất khung tiêu chuẩn thiết kế kinh tế mới ở mức 1/50 năm và 1/100 năm, đi kèm chuẩn an toàn rủi ro xã hội tối thiểu 1/150 năm.
  • Mở rộng thành công mô hình kinh tế van Dantzig ứng dụng cho vùng đồng bằng dốc ven biển tại các quốc gia đang phát triển.

Nghiên cứu là một bước tiến quan trọng trong việc chuyển giao tri thức quốc tế vào thực tiễn đê điều Việt Nam. Trong lộ trình từ năm 2026 đến năm 2030, các cơ quan quản lý và giới nghiên cứu công trình thủy cần khẩn trương áp dụng khung phương pháp luận này vào quy chuẩn quốc gia, nhằm kiến tạo hệ thống phòng hộ ven biển hiện đại, vững chắc và chủ động ứng phó với biến đổi khí hậu toàn cầu.