Tổng quan nghiên cứu

Vùng hạ lưu hệ thống sông Đồng Nai – Sài Gòn với diện tích tự nhiên khoảng 452.000 ha là trung tâm kinh tế, công nghiệp và đô thị năng động bậc nhất cả nước. Nhằm bảo đảm an ninh năng lượng và điều tiết nguồn nước, công trình thủy điện Trị An được xây dựng từ năm 1988 trên lưu vực rộng 14.800 km2, sở hữu dung tích toàn bộ đạt 2,567 tỷ m3 cùng dung tích phòng lũ 635 triệu m3. Công trình đóng vai trò chiến lược trong việc phát điện với công suất 400 MW, đẩy mặn mùa khô và cắt giảm lũ cho vùng hạ du. Tuy nhiên, sau hơn 25 năm vận hành, dưới tác động của biến đổi khí hậu toàn cầu và sự gia tăng các hình thái thời tiết cực đoan, nguy cơ mất an toàn đập trở thành mối hiểm họa tiềm tàng to lớn.

Thực tế lịch sử ngành thủy lợi thế giới ghi nhận hơn 300 sự cố đập lớn trong thế kỷ 20, gây tổn thất nghiêm trọng về người và tài sản, điển hình như thảm họa vỡ đập Vaiont năm 1963 cướp đi sinh mạng của 2.600 người. Nhận thức rõ mức độ nguy hiểm đối với khu vực hạ lưu trải dài 150 km từ chân đập Trị An đến cửa Soài Rạp, luận văn thạc sĩ chuyên ngành Xây dựng công trình thủy của tác giả Lê Ngọc Anh (Trường Đại học Bách Khoa – ĐHQG TP.HCM) đã thực hiện đề tài nghiên cứu ứng dụng mô hình toán số TELEMAC-2D nhằm đánh giá toàn diện ảnh hưởng của hiện tượng vỡ đập Trị An. Mục tiêu cụ thể là mô phỏng quá trình truyền sóng lũ thủy lực hai chiều, định lượng chính xác mực nước lớn nhất, độ sâu ngập, vận tốc dòng chảy và thời gian truyền sóng lũ trong chu kỳ 96 giờ. Kết quả nghiên cứu cung cấp cơ sở khoa học thiết thực để bảo vệ tính mạng cho hơn 10 triệu người dân và bảo vệ các trung tâm kinh tế trọng điểm phía Nam.

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 trên nền tảng cơ học thủy khí hiện đại kết hợp lý thuyết thủy động lực học dòng chảy hở hai chiều không ổn định. Cơ chế thủy lực hạ lưu chịu sự chi phối phức tạp của chế độ bán nhật triều không đều tại Biển Đông với biên độ dao động triều lớn từ 3,0 m đến 4,0 m trong các kỳ triều cường.

Khung lý thuyết cốt lõi ứng dụng hệ phương trình Saint-Venant hai chiều trung bình theo chiều sâu (2DH), bao gồm phương trình liên tục và hệ hai phương trình bảo toàn động lượng theo phương nằm ngang (trục x và y). Mô hình tích hợp đầy đủ các đại lượng vật lý gồm cao trình mặt nước Z, gia tốc trọng trường g = 9,81 m/s2, hệ số nhớt rối động lượng, lực ma sát đáy và thành phần lực quán tính Coriolis.

Bên cạnh đó, đề tài ứng dụng bộ mô hình toán thủy lực mã nguồn mở TELEMAC-2D do Phòng Thí nghiệm Thủy lực và Môi trường Quốc gia (EDF - Pháp) phát triển. TELEMAC-2D tích hợp phương pháp phần tử hữu hạn (Finite Element Method) giải sơ đồ ẩn có tính ổn định cao và phương pháp thể tích hữu hạn phi cấu trúc (Finite Volume Method) áp dụng phép xấp xỉ Godunov giải bài toán Riemann. Cấu trúc thuật toán này cho phép mô phỏng chính xác sóng xung kích gián đoạn sinh ra khi vỡ đập đột ngột trên nền địa hình lòng dẫn phức tạp.

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

Nghiên cứu kết hợp phương pháp mô hình hóa toán số với phân tích không gian GIS (MapInfo, ArcGIS) và phần mềm tạo lưới Blue Kenue để thiết lập sơ đồ tính toán thủy lực:

  • Cỡ mẫu và cấu trúc lưới tính: Miền tính toán hạ lưu Đồng Nai – Sài Gòn được rời rạc hóa thành 171.919 phần tử tam giác phi cấu trúc. Kích thước diện tích phần tử được làm mịn đạt 36 m2 tại các khúc sông hẹp, chân công trình và mở rộng tối đa khoảng 1.500 m2 tại các vùng bãi tràn ngập lũ. Phương pháp chọn lưới phi cấu trúc giúp mô tả chuẩn xác các bờ uốn khúc phức tạp và tối ưu hóa thời gian tính toán so với lưới trực giao truyền thống.
  • Nguồn dữ liệu đầu vào: Địa hình mặt đất được trích xuất từ mô hình số độ cao DEM kích thước 25x25 m, kết hợp dữ liệu bình đồ lòng dẫn sông 2D và mô hình độ sâu Biển Đông lưới 1.800x1.800 m. Chuỗi số liệu thủy văn biên thu thập từ 11 trạm quan trắc thực tế, gồm Trị An, Dầu Tiếng, Phước Hòa, Biên Hòa, Nhà Bè, Phú An, Bến Lức, Tân An, Bến Đá, Mộc Hóa và Vàm Kênh.
  • Hiệu chỉnh và kiểm định: Mô hình được hiệu chỉnh qua trận lũ lịch sử tháng 10/2000 kéo dài 10 ngày (240 giờ) với tổng lưu lượng xả lũ toàn vùng đạt khoảng 5.000 m3/s, đồng thời được kiểm định độc lập với chuỗi số liệu trận lũ tháng 10/2007. TELEMAC-2D được lựa chọn nhờ khả năng xử lý bài toán ngập khô linh hoạt và tốc độ tính toán song song vượt trội so với các phần mềm thương mại đắt đỏ như MIKE 21 hay HEC-RAS.

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

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

Mô phỏng thủy lực sự cố vỡ đập Trị An trong thời gian 96 giờ theo hai kịch bản tổ hợp lũ và triều đã làm sáng tỏ các thông số động lực học nguy hiểm:

  • Phát hiện 1 (Kịch bản 1 - Vỡ đập đơn lẻ): Khi đập vỡ đột ngột trong điều kiện lũ thượng nguồn sông Bé và sông Sài Gòn chưa kịp đổ về, sóng lũ hình thành ngay sau đập tạo lưu tốc dòng chảy cực đại vượt mức 8,0 m/s. Sóng lũ dồn dập làm mực nước tại trạm Biên Hòa tăng vọt vượt đỉnh lũ lịch sử năm 2000 (2,19 m) tới hơn 150%, gây ngập sâu trung bình từ 3,0 m đến 5,5 m tại các vùng trũng ven sông chỉ sau 6 đến 12 giờ truyền lũ.
  • Phát hiện 2 (Kịch bản 2 - Vỡ đập tổ hợp lũ lớn và triều cường): Khi đập vỡ sau 24 giờ kết hợp cùng lúc đỉnh lũ sông Bé đạt trên 1.400 m3/s, hồ Dầu Tiếng xả 600 m3/s và triều cường trạm Vũng Tàu đạt đỉnh 1,45 m, quy mô ngập lụt gia tăng đột biến. Mực nước đỉnh lũ tại trạm Biên Hòa và trạm Nhà Bè dâng cao thêm từ 0,6 m đến 1,2 m so với Kịch bản 1, tương đương mức tăng áp lực ngập từ 25% đến 35%.
  • Phát hiện 3 (Vận tốc và phạm vi lan truyền): Vận tốc truyền sóng lũ dọc dòng chính sông Đồng Nai duy trì ở mức cao từ 2,5 m/s đến 3,8 m/s tại các phân lưu sông Nhà Bè và sông Lòng Tàu. Thời gian sóng lũ di chuyển toàn tuyến dài 150 km ra đến cửa Soài Rạp mất khoảng 28 đến 36 giờ, làm ngập úng nghiêm trọng trên 80% diện tích các khu vực trũng thấp thuộc quận Thủ Đức, quận 2, quận 7, huyện Nhà Bè, Bình Chánh và Cần Giờ.

Thảo luận kết quả

Diễn biến thủy lực khẳng định hiện tượng vỡ đập Trị An sẽ tạo nên một thảm họa thiên tai có sức tàn phá nghiêm trọng đối với toàn vùng hạ lưu. Nguyên nhân chính khiến mức độ ngập lụt gia tăng khốc liệt trong Kịch bản 2 là do sự cộng hưởng thủy lực hai chiều: dòng lũ khổng lồ với dung tích hơn 2,5 tỷ m3 từ hồ Trị An tràn về bị chặn đứng bởi sóng triều cường Biển Đông dâng cao với biên độ triều 3,0 m đến 4,0 m, làm tê liệt khả năng tự tiêu thoát của cửa sông Soài Rạp và vịnh Gành Rái.

So sánh với các trận lũ cực hạn trong quá khứ như trận lũ ngày 20/10/1952 (mực nước trạm Biên Hòa đạt đỉnh 4,73 m) hay trận lũ năm 2000 (tổng lượng lũ 30 tỷ m3), cường suất gia tăng mực nước do vỡ đập lớn gấp 4 đến 6 lần, vượt xa năng lực bảo vệ của các tuyến đê bao đô thị hiện hữu.

Về phương thức trực quan hóa, toàn bộ dữ liệu thủy lực được biểu diễn sinh động qua hệ thống đồ thị đường quá trình mực nước Z(t) tại 12 vị trí mặt cắt đặc trưng, kết hợp bản đồ số 2D thể hiện trường phân bố độ sâu ngập và trường vectơ vận tốc dòng chảy. Bảng so sánh chỉ số ngập chi tiết giữa hai kịch bản giúp các nhà quản lý dễ dàng nhận diện ranh giới phân vùng nguy hiểm trực tiếp và vùng giáp nước giao thoa lũ - triều.

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

Dựa trên kết quả tính toán thủy lực, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược nhằm giảm thiểu tối đa rủi ro thiệt hại:

  1. Thiết lập hệ thống quan trắc và cảnh báo sớm tự động: Lắp đặt mạng lưới cảm biến đo áp lực thấm, độ lún thân đập và đo mưa tự động với tần suất truyền tin 5 phút/lần tại hồ Trị An; bảo đảm phát tín hiệu cảnh báo thảm họa trước tối thiểu 12 đến 24 giờ cho hạ du. Kế hoạch thực hiện trong 18 tháng do Công ty Thủy điện Trị An chủ trì phối hợp cùng Ban Chỉ huy Phòng chống thiên tai tỉnh Đồng Nai.
  2. Quy hoạch không gian thoát lũ và gia cố công trình bờ bao: Nạo vét trục thoát lũ dòng chính sông Đồng Nai – Nhà Bè – Soài Rạp bảo đảm khả năng thông qua lưu lượng tối thiểu 6.000 m3/s; nâng cấp cao trình đê bao xung yếu ven đô thị lên mức an toàn từ 2,5 m đến 3,0 m nhằm chống chịu dòng chảy lưu tốc trên 3,0 m/s. Timeline triển khai giai đoạn 2015–2020 do Ủy ban nhân dân TP.HCM và tỉnh Đồng Nai phụ trách.
  3. Số hóa bản đồ ngập lụt và lập phương án sơ tán dân cư khẩn cấp: Tích hợp bản đồ ngập lụt 2D từ mô hình TELEMAC-2D vào hệ thống thông tin địa lý để tổ chức di dời hơn 200.000 hộ dân tại các vùng ngập sâu trên 2,0 m; khoanh định 15 vị trí gò đồi tự nhiên cao từ 5 m đến 15 m tại Củ Chi, Thủ Đức, Long Bình làm nơi trú ẩn an toàn tuyệt đối. Kế hoạch hoàn thiện trong 12 tháng do cơ quan phòng chống thiên tai địa phương chủ trì.
  4. Tích hợp mô hình số vào quy trình vận hành liên hồ chứa theo thời gian thực: Định kỳ 6 tháng một lần thực hiện cập nhật dữ liệu lòng dẫn và vận hành mô phỏng liên hồ chứa Trị An, Thác Mơ, Dầu Tiếng nhằm chủ động điều tiết xả lũ sớm trước khi lũ đạt đỉnh. Đơn vị thực hiện là Viện Quy hoạch Thủy lợi Miền Nam và các cơ quan nghiên cứu chuyên ngành.

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

Công trình nghiên cứu mang giá trị học thuật và tính ứng dụng thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

  1. Cơ quan Quản lý Nhà nước và Ban Chỉ đạo Phòng chống thiên tai: Sử dụng các số liệu định lượng về độ sâu ngập, lưu tốc và thời gian truyền sóng lũ để xây dựng phương án ứng phó khẩn cấp và bảo vệ an toàn cho 452.000 ha vùng hạ lưu kinh tế trọng điểm.
  2. Các đơn vị quản lý, vận hành hồ chứa thủy điện và thủy lợi: Cung cấp tài liệu chuyên sâu cho đội ngũ kỹ sư tại Thủy điện Trị An, Dầu Tiếng, Thác Mơ nhằm đánh giá rủi ro an toàn đập dung tích trên 2,5 tỷ m3 và tối ưu hóa quy trình xả lũ liên hồ.
  3. Đơn vị tư vấn quy hoạch đô thị và hạ tầng giao thông thủy: Tham khảo bản đồ phân vùng vận tốc dòng chảy cực đại từ 2,5 m/s đến 8,0 m/s để tính toán cốt nền xây dựng đô thị, thiết kế đê kè cảng biển và hạ tầng kỹ thuật tại TP.HCM, Bình Dương và Đồng Nai.
  4. Giảng viên, học viên cao học và chuyên gia ngành Kỹ thuật tài nguyên nước: Khai thác quy trình ứng dụng phần mềm mã nguồn mở TELEMAC-2D, kỹ thuật chia lưới phi cấu trúc 171.919 phần tử và lập trình Fortran chuyên sâu trong giải bài toán thủy lực sóng vỡ đập 2D.

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

Mô hình TELEMAC-2D có ưu điểm gì nổi bật so với các mô hình thủy lực khác trong bài toán vỡ đập?
TELEMAC-2D là mô hình mã nguồn mở miễn phí, hỗ trợ tính toán song song hiệu năng cao và giải bài toán sóng gián đoạn vỡ đập chính xác bằng phương pháp thể tích hữu hạn Godunov trên lưới phi cấu trúc 171.919 phần tử. Nền tảng cho phép lập trình tùy biến linh hoạt bằng ngôn ngữ Fortran, vượt trội về mặt chi phí so với các gói phần mềm thương mại.

Thời gian truyền sóng lũ từ đập Trị An đến trung tâm TP.HCM và cửa biển mất bao lâu?
Kết quả tính toán chỉ ra sóng lũ di chuyển với vận tốc cực đại trên 8,0 m/s ở thượng lưu và giảm xuống 2,5 đến 3,8 m/s ở hạ lưu. Sóng lũ bắt đầu tác động trực tiếp đến TP. Biên Hòa sau 2 đến 4 giờ, lan truyền đến trung tâm TP.HCM sau 10 đến 16 giờ và ra đến cửa Soài Rạp sau 28 đến 36 giờ.

Yếu tố nào làm gia tăng nghiêm trọng mức độ ngập lụt hạ du khi xảy ra sự cố vỡ đập?
Tổ hợp bất lợi nguy hiểm nhất xảy ra khi vỡ đập trùng thời điểm đỉnh lũ sông Bé đạt 1.400 m3/s, hồ Dầu Tiếng xả 600 m3/s và gặp kỳ triều cường Biển Đông với biên độ dao động từ 3,0 m đến 4,0 m. Triều dâng hoạt động như một đập chắn thủy lực, làm mực nước hạ lưu dâng cao thêm từ 0,6 m đến 1,2 m.

Bản đồ ngập lụt trong luận văn hỗ trợ công tác cứu hộ và di dời dân cư như thế nào?
Bản đồ mô phỏng 2D chi tiết tại 12 vị trí mặt cắt giúp lực lượng cứu nạn khoanh vùng chính xác các điểm ngập sâu trên 3,0 m để sơ tán khẩn cấp. Đồng thời, đề tài xác định rõ các khu vực gò đồi cao từ 5 m đến 15 m tại Củ Chi, Thủ Đức, Long Bình làm vị trí lập trạm sơ tán an toàn tuyệt đối.

Làm cách nào để nâng cao độ chính xác mô phỏng trong các nghiên cứu thủy lực tiếp theo?
Cần định kỳ cập nhật mô hình số độ cao DEM với độ phân giải chi tiết hơn 25x25 m và đo đạc lại địa hình đáy sông Đồng Nai – Nhà Bè – Soài Rạp. Sự biến đổi bồi xói đáy sông tự nhiên qua thời gian là nguyên nhân chính dẫn đến sai lệch đỉnh lũ giữa mô phỏng và thực đo tại các trạm kiểm định.

Kết luận

  • Luận văn đã ứng dụng thành công mô hình toán số mã nguồn mở TELEMAC-2D để mô phỏng toàn diện bài toán thủy lực hai chiều khi xảy ra hiện tượng vỡ đập Trị An trên phạm vi lưu vực hạ lưu 452.000 ha.
  • Thiết lập hoàn chỉnh sơ đồ lưới phi cấu trúc gồm 171.919 phần tử, được hiệu chỉnh qua trận lũ tháng 10/2000 và kiểm định chuẩn xác với chuỗi số liệu thực đo năm 2007.
  • Lượng hóa chi tiết hai kịch bản vỡ đập, chỉ ra lưu tốc cực đại vượt 8,0 m/s sau đập, thời gian truyền sóng lũ toàn tuyến từ 28 đến 36 giờ và mức độ ngập lụt tăng 25% đến 35% khi kết hợp triều cường 3,0 m đến 4,0 m.
  • Đề xuất 4 nhóm giải pháp đồng bộ kết hợp giữa gia cố công trình đê bao và giải pháp phi công trình về cảnh báo sớm, quy hoạch thoát lũ và lựa chọn các vùng gò đồi cao từ 5 m đến 15 m làm nơi sơ tán.
  • Định hướng nghiên cứu giai đoạn 2015–2018 cần tập trung tích hợp mô đun bùn cát 2D Sisyphe và xây dựng mô hình dự báo ngập lụt trực tuyến theo thời gian thực.

Hãy liên hệ ngay với các viện nghiên cứu chuyên ngành tài nguyên nước để tiếp nhận chuyển giao công nghệ mô hình toán TELEMAC-2D và ứng dụng bản đồ ngập lụt vào công tác quản trị rủi ro thiên tai tại địa phương!