Tổng quan nghiên cứu

Vùng Đồng Tháp Mười là vùng trũng ngập lũ lớn thuộc đồng bằng sông Cửu Long với tổng diện tích tự nhiên đạt 7.088,6 km², trong đó diện tích đất canh tác nông nghiệp chiếm tới 72% (tương đương 510.379,2 ha) trải rộng trên địa bàn ba tỉnh Long An, Đồng Tháp và Tiền Giang. Hằng năm, khu vực này phải đối mặt với các đợt ngập lụt kéo dài từ 3 đến 5 tháng (thường từ tháng 7 đến tháng 11) do lũ tràn từ biên giới Việt Nam - Campuchia kết hợp với lượng mưa nội vùng, khiến đỉnh lũ có thời điểm làm ngập úng tới 90% diện tích toàn vùng. Trận lũ lịch sử năm 2000 với lưu lượng đỉnh tràn biên giới lên tới 12.000 m³/s (vượt xa mức 8.270 m³/s của trận lũ năm 1996) và mực nước đo được tại trạm Tân Châu đạt mức kỷ lục 5,15 m đã gây thiệt hại hàng trăm triệu USD cho cơ sở hạ tầng và sản xuất nông nghiệp. Vấn đề cốt lõi đặt ra là tình trạng ngập lụt sớm vào tháng 8 đe dọa trực tiếp đến việc thu hoạch lúa vụ Hè - Thu, trong khi lũ rút chậm vào cuối mùa lại làm trì hoãn lịch xuống giống lúa vụ Đông - Xuân vào đầu tháng 11.

Mục tiêu cụ thể của luận văn là ứng dụng mô hình toán thủy lực kết hợp công nghệ không gian để mô phỏng diễn biến lũ, đánh giá mức độ giảm thiểu thiệt hại nông nghiệp của các phương án công trình kiểm soát lũ, từ đó lựa chọn phương án tối ưu về mặt kinh tế - kỹ thuật nhằm bảo vệ bền vững cơ cấu hai vụ lúa mỗi năm. Phạm vi nghiên cứu tập trung vào vùng trũng Đồng Tháp Mười, sử dụng chuỗi dữ liệu thực đo các trận lũ năm 1996, 1999 và 2000. Ý nghĩa của công trình được chứng minh qua việc xác định phương án công trình có tỷ suất hoàn vốn nội bộ (IRR) đạt 34,2% và tỷ số lợi ích trên chi phí (B/C) đạt 2,00 ở mức chiết khấu 10%, tạo cơ sở khoa học tin cậy cho công tác quy hoạch phát triển thủy lợi vùng.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng lý thuyết dòng chảy không ổn định một chiều trong mạng lưới sông ngòi và bãi ngập lũ thông qua hệ phương trình Saint-Venant đầy đủ. Hệ phương trình bao gồm phương trình liên tục phản ánh sự bảo toàn khối lượng và phương trình động lượng phản ánh định luật bảo toàn động lượng, giải quyết sự thay đổi phức tạp của lưu lượng, diện tích mặt cắt ướt và mực nước theo không gian và thời gian. Bên cạnh đó, khung phân tích kinh tế dự án thủy lợi được tích hợp với lý thuyết phân tích chi phí - lợi ích (Cost-Benefit Analysis - CBA), ứng dụng các chỉ tiêu kinh tế kỹ thuật như Giá trị hiện tại thuần (NPV), Tỷ số Lợi ích - Chi phí (B/C) và Suất thu lợi nội tại (IRR) để thẩm định tính khả thi tài chính.

Các khái niệm chuyên ngành then chốt được chuẩn hóa bao gồm: ô chứa ngập lũ nội đồng (storage cells) mô tả các vùng trũng trữ nước được phân cách bởi đê và kênh trục; hệ số nhám Manning (n) đặc trưng cho sức cản lòng dẫn với giá trị hiệu chỉnh tối ưu là 0,05; sơ đồ sai phân hữu hạn ẩn 4 điểm Preissmann đảm bảo độ ổn định và hội tụ cao khi tính toán mạng lưới thủy lực sông kênh phức tạp; và bề mặt mô hình ngập lụt (WSML) phục vụ trích xuất độ sâu ngập tức thời.

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

Nguồn dữ liệu đầu vào được thu thập đồng bộ từ các cơ quan khí tượng thủy văn và quy hoạch thủy lợi miền Nam, bao gồm: chuỗi số liệu quan trắc mưa 34 năm tại trạm Tân Châu (giai đoạn 1970–2004), số liệu lưu lượng ngày tại 4 trạm biên giới và sông chính (Kratie, Biển Hồ, Vàm Cỏ Đông, Vàm Cỏ Tây), số liệu mực nước giờ tại 8 trạm thủy văn nội vùng, và số liệu mực nước triều tại 3 cửa sông lớn (Vàm Kênh, Bình Đại, Bến Trại). Cỡ mẫu không gian được sơ đồ hóa thành 15 ô đồng đại diện cho 510.379,2 ha đất canh tác nông nghiệp. Phương pháp phân vùng ô đồng dựa trên địa hình tự nhiên, ranh giới các tuyến đường giao thông huyết mạch và hướng tiêu thoát lũ tự nhiên ra sông Tiền và sông Vàm Cỏ.

Lý do lựa chọn mô hình thủy lực VRSAP (Vietnam River System and Plain Model) kết hợp phần mềm ArcView GIS là do tính tương thích cao với đặc thù thủy văn châu thổ sông Mekong. Mô hình được hiệu chỉnh dựa trên số liệu trận lũ năm 1996 với chỉ số hiệu quả Nash-Sutcliffe (EI) tại các trạm kiểm tra dao động từ 0,91 đến 0,97 và sai số đỉnh lũ (∆Hmax) chỉ từ 0,016 m đến 0,107 m. Quá trình kiểm định trên trận lũ năm 1999 tiếp tục ghi nhận kết quả rất tốt với chỉ số EI từ 0,89 đến 0,96 và ∆Hmax từ 0,025 m đến 0,077 m. Sau khi kiểm định đạt chuẩn, mô hình tiếp tục mô phỏng kịch bản ngập lũ lịch sử năm 2000 kết hợp mô hình số độ cao DEM độ phân giải 100x100 m để tính toán độ sâu ngập, diện tích ngập và thời gian ngập lụt cho từng phương án.

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

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

Thứ nhất, mô phỏng thủy lực kịch bản hiện trạng (không có công trình) cho thấy mực nước đỉnh lũ năm 2000 tại các nút dọc sông Tiền giảm dần từ thượng lưu về hạ lưu: tại Tân Châu đạt 5,15 m, Chợ Mới đạt 3,58 m, Cao Lãnh đạt 2,65 m và Mỹ Thuận đạt 1,95 m. Mực nước dâng cao đột ngột từ ngày 20/8 và duy trì ngập sâu suốt tháng 9 và tháng 10, khiến 90% diện tích đất nông nghiệp bị ngập sâu từ 1,0 m đến 3,0 m, gây thiệt hại nghiêm trọng nếu lúa vụ Hè - Thu chưa thu hoạch xong trước ngày 30/8.

Thứ hai, trong nhóm các giải pháp tiêu thoát lũ qua hệ thống kênh trục ngang, phương án nạo vét mở rộng kênh An Phong - Mỹ Hòa - Bắc Đông (bề rộng đáy B=14 m, cao trình đáy H=-3,0 m, mái dốc m=1,5) kết hợp cống kiểm soát cuối kênh mang lại hiệu quả thủy lực vượt trội. Mực nước đỉnh lũ dọc tuyến kênh An Phong giảm từ 18 cm đến 22 cm, kênh Đồng Tiến giảm 14 cm đến 19 cm, kênh Hồng Ngự giảm 13 cm đến 14 cm và kênh Nguyễn Văn Tiếp giảm từ 10 cm đến 15 cm so với khi chưa nạo vét.

Thứ ba, việc hạ thấp mực nước nội đồng vào thời điểm then chốt ngày 30/8 giúp giảm hơn 65% diện tích ngập úng nguy hiểm tại các ô đồng trọng điểm, bảo vệ an toàn cho năng suất lúa vụ Hè - Thu đạt mức trung bình 4,5 tấn/ha. Đồng thời, vào thời điểm ngày 01/11, mực nước tại các ô đồng số 1, 3, 5 và 7 rút nhanh hơn từ 10 đến 25 cm, cho phép nông dân xuống giống vụ Đông - Xuân sớm hơn từ 10 đến 15 ngày, loại bỏ hoàn toàn nguy cơ chậm thời vụ.

Thứ tư, tính toán kinh tế cho thấy phương án nạo vét kênh An Phong - Mỹ Hòa kết hợp cống điều tiết đạt tỷ số Lợi ích / Chi phí B/C = 2,00 ở mức lãi suất chiết khấu 10%, với tỷ suất hoàn vốn nội bộ IRR đạt 34,2%. Khi phân tích độ nhạy ở các mức lãi suất chiết khấu 6%, 8%, 12% và 14%, tỷ số B/C tương ứng đạt mức từ 3,61 xuống 1,52, khẳng định tính bền vững và khả thi tài chính vượt trội so với các phương án nạo vét kênh trục dọc song song sông Tiền.

Thảo luận kết quả

Cơ chế mang lại hiệu quả tiêu thoát nước cao của phương án nạo vét kênh An Phong - Mỹ Hòa bắt nguồn từ việc mở rộng trục thoát ngang chính, tạo hành lang chuyển hướng nhanh dòng nước lũ từ vùng trũng trung tâm ra hệ thống sông Vàm Cỏ Tây và Vàm Cỏ Đông. Việc bố trí công trình cống ở cuối kênh giữ vai trò then chốt: mở tự do để tháo lũ nhanh trong mùa mưa và đóng kín để ngăn mặn xâm nhập vào mùa khô, đồng thời duy trì mực nước trữ cho các ô đồng khi cần thiết.

Dữ liệu mô phỏng không gian được thể hiện trực quan qua chuỗi bản đồ ngập lụt GIS và bảng tổng hợp dòng tiền tài chính theo các giai đoạn đầu tư. Trên các bản đồ phân vùng ngập úng ngày 30/8 và 01/11, sự thu hẹp rõ rệt của các dải màu thể hiện độ sâu ngập trên 1,0 m minh chứng cho khả năng kiểm soát ngập lụt cục bộ của công trình. So với các nghiên cứu trước đây tại đồng bằng sông Cửu Long vốn ưu tiên xây dựng đê bao triệt để, kết quả của luận văn chứng minh rằng giải pháp công trình bán kiểm soát kết hợp tiêu thoát qua kênh trục ngang vừa tối ưu hóa chi phí đầu tư ban đầu, vừa hạn chế tình trạng dâng cao mực nước lũ ở vùng hạ lưu, đồng thời bảo tồn trọn vẹn nguồn lợi phù sa và thủy sản mùa lũ tự nhiên.

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

Dựa trên các kết quả tính toán thủy lực và hiệu quả kinh tế, bốn nhóm giải pháp chiến lược được đề xuất triển khai:

Thứ nhất, ưu tiên bố trí nguồn vốn đầu tư công thực hiện nạo vét mở rộng tuyến kênh trục ngang An Phong - Mỹ Hòa - Bắc Đông đạt thông số kỹ thuật chuẩn (bề rộng đáy B=14 m, cao trình đáy H=-3,0 m, hệ số mái dốc m=1,5) và xây dựng cống điều tiết cuối kênh. Mục tiêu hạ thấp mực nước lũ nội đồng từ 18 đến 22 cm, hoàn thành trong giai đoạn 2006–2010 do Bộ Nông nghiệp và Phát triển Nông thôn chủ trì phối hợp cùng Ủy ban nhân dân tỉnh Đồng Tháp và Long An.

Thứ hai, nâng cấp và gia cố toàn diện tuyến đê bao kiểm soát lũ sớm Tân Thành - Lò Gạch đạt cao trình đỉnh +6,5 m tại khu vực Hồng Ngự và +5,5 m tại khu vực Vĩnh Hưng, kết hợp nạo vét mở rộng lòng kênh biên giới với bề rộng đáy từ 24 m đến 32 m. Mục tiêu kiểm soát an toàn dòng lũ tràn qua biên giới vượt mốc 12.000 m³/s trước thời điểm ngày 30/8 hằng năm, thực hiện trong giai đoạn 2008–2012 do Ban Quản lý Đầu tư và Xây dựng Thủy lợi 10 đảm trách.

Thứ ba, thiết lập hệ thống quan trắc tự động mực nước, lưu lượng và vận hành thời gian thực trên nền tảng tích hợp mô hình VRSAP-GIS cho toàn bộ hệ thống cống đầu mối tại Đồng Tháp Mười. Mục tiêu nâng độ chính xác của công tác cảnh báo ngập lụt và điều tiết cống lên trên 95%, triển khai trong giai đoạn 2010–2015 do các Trung tâm Quản lý và Khai thác Công trình Thủy lợi địa phương thực hiện.

Thứ tư, tái cơ cấu lịch mùa vụ nông nghiệp thích ứng linh hoạt với lũ, duy trì ổn định mô hình 2 vụ lúa (Hè - Thu và Đông - Xuân) kết hợp mở cống đón lũ lấy phù sa trong tháng 9 và tháng 10. Mục tiêu giảm 40% chi phí phân bón hóa học đầu vào và gia tăng 15% thu nhập nông hộ từ đánh bắt thủy sản, do Sở Nông nghiệp và Phát triển Nông thôn các tỉnh phối hợp cùng Hội Nông dân hướng dẫn thực hiện thường niên.

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

Luận văn là tài liệu tham khảo khoa học có giá trị cao cho bốn nhóm đối tượng chính sau:

Thứ nhất, các kỹ sư thủy lợi và chuyên viên quy hoạch tài nguyên nước tại các Viện Quy hoạch Thủy lợi, Ban Quản lý Dự án Nông nghiệp. Luận văn cung cấp quy trình chi tiết về việc ứng dụng hệ phương trình Saint-Venant, kỹ thuật phân chia ô đồng và phương pháp phân tích chi phí - lợi ích (CBA) phục vụ lập báo cáo tiền khả thi cho các công trình tiêu thoát lũ quy mô lớn.

Thứ hai, các nhà hoạch định chính sách và cơ quan quản lý nông nghiệp cấp tỉnh, huyện tại Đồng Tháp, Long An và Tiền Giang. Dữ liệu ngập lụt theo từng mốc thời gian ngày 30/8 và 01/11 cung cấp luận cứ thực tiễn để ban hành khung lịch thời vụ xuống giống an toàn, tránh thất thoát mùa vụ do thiên tai.

Thứ ba, giảng viên, nghiên cứu sinh và học viên cao học thuộc các chuyên ngành Kỹ thuật Tài nguyên nước, Thủy văn học, Quản lý Môi trường. Công trình là hình mẫu chuẩn mực trong việc kết hợp mô hình thủy lực 1D (VRSAP) với hệ thống thông tin địa lý (ArcView GIS) và mô hình số độ cao DEM để trực quan hóa không gian ngập lụt.

Thứ tư, các chuyên gia thuộc các tổ chức phi chính phủ và đối tác quốc tế (như Danida, CIDA, Ủy hội sông Mekong quốc tế MRC). Luận văn mang lại nguồn dữ liệu kiểm chứng độc lập về quản lý tổng hợp tài nguyên nước và quản trị rủi ro ngập lụt thích ứng với biến đổi khí hậu ở hạ lưu sông Mekong.

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

1. Mô hình VRSAP có ưu điểm gì nổi bật khi mô phỏng mạng lưới thủy lực vùng Đồng Tháp Mười? Mô hình VRSAP giải hệ phương trình Saint-Venant một chiều đầy đủ bằng sơ đồ ẩn 4 điểm Preissmann, cho phép mô phỏng chính xác dòng chảy không ổn định trong mạng lưới sông kênh chằng chịt và tương tác trao đổi nước với các ô đồng trũng. Trong nghiên cứu này, mô hình đạt chỉ số hiệu quả Nash-Sutcliffe từ 0,91 đến 0,97 và sai số đỉnh lũ dưới 0,107 m.

2. Vì sao phương án nạo vét kênh An Phong - Mỹ Hòa lại được đánh giá tối ưu hơn các phương án khác? Phương án này tạo hành lang thoát lũ ngang chủ đạo giúp hạ mực nước nội đồng từ 18 cm đến 22 cm, bảo vệ an toàn vụ lúa Hè - Thu trước 30/8 và rút nước sớm trước 01/11 cho vụ Đông - Xuân. Về kinh tế, phương án đạt tỷ số B/C = 2,00 ở mức lãi suất 10% và suất hoàn vốn nội bộ IRR lên tới 34,2%.

3. Vai trò của việc kết hợp mô hình VRSAP và phần mềm ArcView GIS trong nghiên cứu là gì? Mô hình VRSAP tính toán mực nước chính xác tại các nút và ô đồng, sau đó ArcView GIS sử dụng mô hình số độ cao DEM (lưới 100x100 m) để nội suy bề mặt mực nước (WSML) và trích xuất độ sâu ngập theo không gian. Cách tiếp cận này giúp xác định chính xác diện tích ngập để định lượng thiệt hại nông nghiệp.

4. Lũ lụt tại vùng Đồng Tháp Mười mang lại những lợi ích tích cực nào cho sản xuất nông nghiệp? Lũ mang lại nguồn phù sa màu mỡ bồi đắp đồng ruộng, giúp rửa trôi độc chất phèn nhôm trong đất chua phèn và tiêu diệt mầm mống sâu bệnh hại lúa. Ngoài ra, dòng lũ còn cung cấp nguồn lợi thủy sản tự nhiên phong phú, giúp các nông hộ gia tăng đáng kể thu nhập phụ trong khoảng thời gian giữa hai vụ lúa.

5. Độ nhạy kinh tế của dự án thay đổi như thế nào khi biến động mức lãi suất chiết khấu? Phân tích độ nhạy chứng minh tính bền vững tài chính vượt trội của dự án: khi lãi suất chiết khấu thay đổi ở các mức 6%, 8%, 10%, 12% và 14%, tỷ số Lợi ích / Chi phí (B/C) lần lượt đạt 3,61; 2,58; 2,00; 1,71 và 1,52. Tỷ số B/C luôn lớn hơn 1 trong mọi kịch bản lãi suất, khẳng định phương án hoàn toàn khả thi về mặt kinh tế.

Kết luận

  • Luận văn đã thiết lập thành công mô hình tích hợp VRSAP và ArcView GIS (DEM 100x100 m) với độ tin cậy rất cao (chỉ số hiệu quả Nash-Sutcliffe đạt 0,91–0,97) để mô phỏng chính xác diễn biến ngập lụt vùng Đồng Tháp Mười.
  • Đánh giá định lượng tác động của trận lũ lịch sử năm 2000 (lưu lượng đỉnh 12.000 m³/s, mực nước Tân Châu 5,15 m) đến 510.379,2 ha đất sản xuất nông nghiệp tại 15 ô đồng trũng.
  • Khẳng định phương án nạo vét mở rộng hệ thống kênh An Phong - Mỹ Hòa - Bắc Đông kết hợp cống kiểm soát cuối kênh là giải pháp tối ưu nhất, giúp hạ mực nước đỉnh lũ từ 18 cm đến 22 cm và đẩy nhanh tốc độ rút nước đầu tháng 11.
  • Chứng minh tính khả thi kinh tế vượt bậc của công trình với tỷ số Lợi ích / Chi phí B/C đạt 2,00 (ở mức chiết khấu 10%), B/C đạt 3,61 (ở mức 6%) và suất thu lợi nội tại IRR đạt 34,2%.
  • Đề xuất lộ trình đầu tư công trình cụ thể cho giai đoạn 2006–2015 gắn liền với việc chuyển đổi lịch thời vụ nông nghiệp an toàn thích ứng bền vững với lũ lụt châu thổ.

Đóng góp quan trọng nhất của luận văn là cung cấp khung phương pháp luận kết hợp thủy lực học công trình và phân tích kinh tế thực chứng, giải quyết hài hòa bài toán kiểm soát ngập lụt nông nghiệp gắn liền với bảo tồn lợi ích sinh thái tự nhiên. Lộ trình tiếp theo đòi hỏi các địa phương nhanh chóng triển khai lập báo cáo kỹ thuật chi tiết cho tuyến kênh An Phong - Mỹ Hòa và tích hợp công nghệ quan trắc tự động vào quy trình vận hành cống. Các nhà quản lý, kỹ sư thủy lợi và nhà nghiên cứu quan tâm có thể khai thác trực tiếp cơ sở dữ liệu và thuật toán của luận văn để ứng dụng hiệu quả vào công tác quy hoạch tài nguyên nước tại lưu vực sông Mekong.