ĐÁNH GIÁ TÁC ĐỘNG PHÂN LŨ CỦA VÙNG ĐỒNG THÁP MƯỜI CHO ĐỒNG BẰNG SÔNG CỬU LONG

Kiểm soát lũ khu vực đồng bằng sông Cửu Long (ĐBSCL) là một vấn đề đã và đang được nhiều tổ chức, cá nhân quan tâm nghiên cứu vì tính chất phức tạp và quan trọng của nó. Nội dung nghiên cứu trong bài báo tập trung vào xem xét khả năng và tác động phân lũ vào vùng Đồng Tháp Mười (ĐTM) đối với vùng hạ lưu. Hai phương án vị trí điển hình trên hệ thống đê bao dọc bờ trái sông Tiền vùng thượng lưu được xem xét để cho phép nước tràn qua đỉnh đê. Giả thiết được xem xét là cao trình đê bao là 2,5 m, diện tràn dài 2 km cho cả hai phương án xem xét. Bài toán được mô phỏng bởi mô hình toán số 2D theo phương ngang nhằm đánh giá sự thay đổi mực nước cũng như lưu lượng tại một số vị trí điển hình trong trường hợp có và không có phân lũ vào khu vực vùng ĐTM. Kết quả tính toán cho thấy sự giảm mực nước khá rõ ràng cho vùng ngay sau vị trí phân lũ và giảm dần tác động khi càng về hạ lưu.

ĐÁNH GIÁ TÁC ĐỘNG PHÂN LŨ CỦA VÙNG ĐỒNG THÁP MƯỜI CHO ĐỒNG BẰNG SÔNG CỬU LONG Nguyễn Thống*, Lưu Xuân Lộc*, Hồ Tuấn Đức* * Khoa Kỹ thuật Xây Dựng- Trường Đại học Bách Khoa, TP.HCM TÓM TẮT Kiểm soát lũ khu vực đồng sông Cửu Long (ĐBSCL) vấn đề nhiều tổ chức, cá nhân quan tâm nghiên cứu tính chất phức tạp quan trọng Nội dung nghiên cứu báo tập trung vào xem xét khả tác động phân lũ vào vùng Đồng Tháp Mười (ĐTM) vùng hạ lưu Hai phương án vị trí điển hình hệ thống đê bao dọc bờ trái sông Tiền vùng thượng lưu xem xét phép nước tràn qua đỉnh đê Giả thiết xem xét cao trình đê bao 2,5 m, diện tràn dài km cho hai phương án xem xét Bài toán mô mô hình toán số 2D theo phương ngang nhằm đánh giá thay đổi mực nước lưu lượng số vị trí điển hình trường hợp có phân lũ vào khu vực vùng ĐTM Kết tính toán cho thấy giảm mực nước rõ ràng cho vùng sau vị trí phân lũ giảm dần tác động hạ lưu Từ khoá: Đồng Tháp Mười, Phân lũ, Telemac2D, Thủy động lực học, Đồng sông Cửu Long IMPACT ASSESSMENT OF FLOOD DIVERSION OF DONG THAP MUOI FOR THE MEKONG DELTA Abstract Flood management in the Mekong Delta draws attention to many organizations and individuals because of its complexity and importance This paper focuses on the capacity and impact of flood diversion of Dong Thap Muoi for the Mekong Delta Along the dykes on the left of Tien River upstream, two typical locations allowing floodwater overtopping were considered, with an assumption of 2.5 m dyke‘s height and km long overtoppable dam The problem is modeled by horizontal 2D numerical modeling to estimate the variation of stage and discharge on some typical positions in the computational area, with and without flood diversion Calculation results shows a significant decline of water level in the area near the dykes, but this effect decreases downstream Keywords: Dong Thap Muoi, Flood diversion, Telemac2D, Hydrodynamic, Mekong River’s Delta ĐẶT VẤN ĐỀ Kiểm soát mực nước lưu lượng vị trí quan trọng vùng đồng sông Cửu Long mùa lũ quan trọng nhằm giảm thiểu tác hại lũ gây Hệ thống công trình thủy lợi nhằm chủ động thoát lũ biển Tây biện pháp xem xét thực bước theo quy hoạch Bên cạnh xem xét khả phương án tạm trữ nước vị trí hợp lý khu vực ĐBSCL, ví dụ vùng Đồng Tháp Mười (ĐTM), khả xem xét Với phương án này, việc giảm đỉnh lũ cho vùng hạ lưu, dung tích tạm trữ khu vực thích hợp vùng ĐBSCL có tác dụng điều hòa dòng chảy vào mùa nước Đây ưu điểm so với phương án xem xét thoát lũ trực tiếp triệt để biển Tây Với hệ thống đê bao hữu khu vực ĐBSCL có tác dụng bảo vệ chống ngập cho khu vực cần bảo vệ mùa lũ Bên cạnh đó, công trình làm giảm khả phân thoát lũ phương diện tổng thể cho toàn khu vực ĐBSCL Đây tác động xấu công trình đê bao cục tiếp tục xem xét xây dựng vùng mà quy hoạch tổng thể cần xem xét lưu ý Trong khuôn khổ nghiên cứu này, giới hạn xem xét tác động phân lũ vùng tạm trữ tiềm lên thay đổi cao trình mực nước lưu lượng cho vùng hạ lưu Sự phân lũ từ sông Tiền vào vùng ĐTM mô thông qua chảy tràn qua tuyến đê dọc bờ trái sông Tiền mực nước sông vượt qua cao trình 2,5m (do giới hạn nội dung báo, việc xem xét với cao trình khác chưa xem xét đây) Về tác dụng phân lũ vai trò điều tiết nước mùa lũ mùa khô khu vực ĐBSCL khu tạm trữ xem xét nghiên cứu khác Nghiên cứu chế độ thủy động lực học miền tính sử dụng mô hình toán số 2D theo phương ngang Telemac2d Một số đê bao miền tính xác định vị trí địa lý tương ứng đỉnh đê định nghĩa nhờ vào đường cong khép kín tạo lưới miền tính đỉnh đê có tác dụng tách miền tạo lưới khu vực tính (trên đỉnh đê lưới miền tính) Các vị trí dự kiến nước chảy tràn cao trình mực nước vượt qua cao trình xác định mô dòng chảy qua đập tràn đa giác với mặt thượng lưu hạ lưu xác định cặp điểm đối diện đường cong khép kín Chế độ thủy lực qua đê đập ngăn (mực nước trước sau đê thấp cao trình đỉnh đê), chảy ngập không ngập phụ thuộc vào cao trình mực nước hai bên đê cao trình đỉnh đê Đập tràn xác định với thông số: hệ số dòng chảy, chiều rộng tràn, cao trình ngưỡng tràn thông qua lập trình từ mã nguồn mở tập tin điều khiển Telemac2D CƠ SỞ LÝ THUYẾT Telemac2D [3], [4] phần mềm thủy lực Tập đoàn Điện lực Pháp (EDF) chủ trì với tham gia nhiều tổ chức nghiên cứu giới dùng để mô dòng chảy chiều theo phương nằm ngang (trung bình theo phương thẳng đứng) mô tả hệ phương trình Saint Venant sau:  h  div hU  q t Z u u u u v  g s  Fx  divh e grad (u ) t x y x h Z v v v u v  g s  Fy  divh e grad ( v) t x y y h   (1) (2) (3)  Trong đó: h(m) – chiều sâu, u & v(m/s) – thành phần vận tốc theo phương ngang x & y vận tốc U , q(m/s) – lưu lượng đơn vị nguồn, Zs (m) – cao độ mặt thoáng, Fx,y(m/s2) – ngoại lực (không kể trọng lực, ví dụ lực Coriolis, ) tác dụng đơn vị khối lượng chiếu theo phương ngang x & y,  e (m2/s)- hệ số khuếch tán Mô hình Telemac2D lập trình lựa chọn phương pháp phần tử hữu hạn thể tích hữu hạn Bài toán số lập trình xử lý dạng điều kiện biên như: Biên mực nước Z, biên lưu lượng Q, biên lưu lượng mực nước (Q & Z), biên vận tốc (u,v), biên vận tốc mực nước (u, v & Z) dạng biên sóng đến Các dạng biên khác có người sử dụng lập trình riêng nhờ vào phần mềm tổ chức dạng mã nguồn mở Chế độ thủy lực chảy qua đê mô chảy tràn qua đập đa giác với lưu lượng đơn vị xác định bởi: 3/ Q   2g Z tl  Zđt  Trường hợp chảy không ngập: (4) 3 1/ (5)  2g Z hl  Z đt Z tl  Z hl  Z  2Z tl Tiêu chuẩn chảy không ngập: (6) Z hl  đt Q(m3/s/m) – lưu lượng tràn đơn vị, –hệ số lưu lượng (0,40,5), Ztl, Zhl (m) – mực nước thượng lưu, hạ lưu, Zđt(m) – cao trình đỉnh tràn Q Trường hợp chảy ngập: MÔ PHỎNG DÒNG CHẢY 2D VÙNG VEN BIỂN VÀ KHU VỰC ĐBSCL 3.1: Miền tính Miền nghiên cứu lấy toàn vùng ĐBSCL Giới hạn phía thượng lưu sông Mékong có vị trí khoảng 200 km kể từ Châu Đốc Phía hạ lưu miền tính lấy rộng khơi cách bờ trung bình khoảng 60 km để giảm thiểu ảnh hưởng đáy biển vùng ven bờ lên triều thiên văn Tại vị trí khơi miền tính, giá trị triều thiên văn áp dụng xem điều kiện biên cho mô Diện tích miền tính khoảng 80850 km2 (hình 1) Toàn miền tính mô hình hoá 716 nghìn phần tử tam giác phi cấu trúc với phần tử tam giác nhỏ có cạnh khoảng 20 m mô tả hai chiều vùng sông rạch phần tử lớn có cạnh khoảng 6000 m mô tả phần tử vùng biển miền tính Hình 1: Lưới miền tính 2D (716 nghìn phần tử tam giác) 3.2: Điều kiện biên địa hình miền nghiên cứu Miền nghiên cứu giới hạn biên thượng lưu biên hạ lưu cụ thể sau: - Biên Q(t) nhập lưu tràn từ biên giới Campuchia qua vùng Đồng Tháp Mười - Biên Q(t) sông Mékong - Biên Q(t) nhập lưu tràn từ biên giới Campuchia qua ngã kênh Vĩnh Tế - Biên Z(t) & U, V biển (biên hở) xác định từ sở liệu triều thiên văn khu vực tương ứng Địa hình khu vực đất liền lấy từ đồ DEM có khoảng cách điểm đo 25 m Địa hình đáy biển lấy từ sở liệu ETOPO1 [7] có khoảng cách trung bình điểm đo 1800 m 3.3: Hiệu chỉnh mô hình 1.50 1.00 Z (m) Nội dung nghiên cứu đánh giá xem xét chế độ thủy động lực học khu vực xem xét ảnh hưởng mực nước dâng biển lên chế độ mực nước khu vực Nội dung báo giới hạn xem xét thay đổi mực nước mùa lũ, mô hình hiệu chỉnh kiểm định cho thời đoạn mùa nhiều nước Trên sở số liệu quan trắc được, mô hình hiệu chỉnh chu kỳ có mực nước triều tương đối cao từ ngày 23/9/2011 đến 26/9/2011 Các trạm đo mực nước Z lưu lượng Q vị trí đo Cần Thơ, Mỹ Thuận, Tân Châu Châu Đốc so sánh với kết tính từ mô hình thời kỳ Điều kiện biên cho thời kỳ mô sau: - Q(t) =400 m3/s nhập lưu tràn từ biên giới Campuchia qua vùng Đồng Tháp Mười; - Q(t) =32000 m3/s sông Mékong; - Q(t) =600 nhập lưu tràn từ biên giới Campuchia qua ngã kênh Vĩnh Tế Mực nước Z biên hở vị trí khơi biển Đông lấy theo triều thiên văn thời kỳ lấy với sóng (S2, N2, K2, K1, O1, P1 Q1) từ sở liệu OTIS [8] có độ phân giải 1/300 Đồ thị sau trình bày mực nước vị trí điển hình A1 A2 biên miền tính (hình 2) sau: A1 A2 0.50 0.00 -0.50 -1.00 13 19 25 31 37 43 49 55 61 67 73 T (h) -1.50 -2.00 Hình 2: Biên triều thiên văn vị trí A1 A2 Quy luật ma sát đáy lấy theo quy luật Chézy Giá trị thay đổi theo thời gian không gian phụ thuộc vào chiều sâu nước vị trí xem xét Đồ thị sau trình bày phân bố giá trị hệ số Chézy thời điểm điển hình mô (hình 3): Hình 3: Hệ số Chézy điển hình cho miền tính Mô hình hiệu chỉnh thông qua chuỗi giá trị lưu lượng mực nước thực đo từ 23/9/2011 đến 26/9/2011 vị trí Cần Thơ, Mỹ Thuận, Tân Châu Châu Đốc Kết trình bày đồ thị từ Hình đến Hình 11 sau: 2.0 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 20000 16000 12000 8000 4000 25 37 49 61 T (h) 73 Hình 4: Lưu lượng Cần Thơ 20000 16000 Quan trắc 12000 8000 4000 13 25 37 49 61 T (h) 73 2.0 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 Mô Q (m3/s) 8000 Quan trắc 4000 2000 T (h) 13 25 37 49 61 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 20000 Q (m3/s) 24000 Mô Quan trắc 12000 8000 4000 T (h) 13 25 37 49 Mô 13 25 37 49 61 73 Quan trắc Mô 13 25 37 49 61 73 Hình 9: Mực nước Châu Đốc 16000 Quan trắc T (h) 73 61 T (h) 73 Hình 8: Lưu lượng Châu Đốc 28000 49 Hình 7: Mực nước Mỹ Thuận Hình 5: Mực nước Cần Thơ 6000 37 T (h) Hình 6: Lưu lượng Mỹ Thuận Hình 4: Lưu lượng Cần Thơ 10000 25 Mô Hình 5: Mực nước Cần Thơ Mô Q (m3/s) 24000 13 Z (m) 13 Z (m) Quan trắc 61 73 5.0 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 Z (m) Q (m3/s) Quan trắc Z (m) Mô 24000 Quan trắc Mô T (h) 13 25 37 49 61 73 Hình 11: Mực nước Tân Châu Hình 10: Lưu lượng Tân Châu Để so sánh đánh giá kết tính từ mô hình số liệu quan trắc, số tiêu chí đánh giá độ thích hợp (GOF) kết tính từ mô hình giá trị quan trắc thực bảng 3: Bảng 1: Một số tiêu so sánh GOF cho lưu lượng Q QTB (mô phỏng) Sai biệt đỉnh TB Sai biệt đáy TB 3 (m /s) (m /s) (m /s) (m3/s) 17248,9 17589,8 1307,6 1423,5 14786,0 14713,5 2658,3 884,8 8291,9 8930,0 562,6 784,1 22561,8 23294,9 638,5 894,6 QTB (quan trắc) Q Cần Thơ Mỹ Thuận Châu Đốc Tân Châu Bảng 2: Một số tiêu so sánh GOF cho mực nước Z ZTB (quan trắc) ZTB (mô phỏng) Sai biệt đỉnh TB Sai biệt đáy TB (m) (m) (m) (m) 1,02 1,18 0,11 0,06 0,95 1,21 0,21 0,04 3,86 4,20 0,32 0,36 4,42 4,42 0,65 0,62 Z Cần Thơ Mỹ Thuận Châu Đốc Tân Châu So sánh kết cho mô hình giá trị quan trắc nêu cho thấy kết tính từ mô hình có khả mô tốt chế độ dòng chảy khu vực chu kỳ tính toán xem xét 20000 Mô Q (m3/s) 24000 Quan trắc 16000 12000 8000 4000 13 25 37 49 Hình 12: Lưu lượng Cần Thơ 61 T (h) 73 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 Z (m) 3.4: Kiểm định mô hình Trên sở mô hình hiệu chỉnh nêu trên, thực bước kiểm định để đánh giá khả mô lại chế độ thủy động lực học khu vực Giai đoạn kiểm định mô hình thực khoảng thời gian ngày từ 27/9/2013 đến 30/9/2013 Điều kiện biên cho thời kỳ mô sau: - Q(t) =200 m3/s nhập lưu tràn từ biên giới Campuchia qua vùng Đồng Tháp Mười; - Q(t) =36000 m3/s sông Mékong tràn bờ kế cận; - Q(t) =300 nhập lưu tràn từ biên giới Campuchia qua ngã kênh Vĩnh Tế - Biên Z biển lấy theo triều thiên văn thời kỳ Số liệu so sánh ghi nhận vị trí Cần Thơ, Mỹ Thuận, Tân Châu Châu Đốc Kết trình bày đồ thị từ Hình 12 đến Hình 19 sau: Quan trắc Mô T (h) 13 25 37 49 61 Hình 13: Mực nước Cần Thơ 73 16000 Quan trắc 12000 8000 4000 13 25 37 49 61 T (h) 73 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 -0.2 5000 Mô Q (m3/s) 6000 Quan trắc 3000 2000 1000 T (h) 13 25 37 49 61 73 3.6 3.2 2.8 2.4 2.0 1.6 1.2 0.8 0.4 0.0 16000 Quan trắc 8000 4000 13 25 37 Hình 18: Lưu lượng Tân Châu 49 61 73 T (h) Mô 13 25 37 49 61 73 Hình 17: Mực nước Châu Đốc 12000 49 Quan trắc 61 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 T (h) 0.0 73 Quan trắc Z (m) Q (m3/s) 20000 Mô 37 T (h) Hình 16: Lưu lượng Châu Đốc 24000 25 Mô Hình 15: Mực nước Mỹ Thuận 4000 13 Z (m) Hình 14: Lưu lượng Mỹ Thuận 7000 Quan trắc Z (m) Mô Q (m3/s) 20000 13 25 37 Mô 49 61 T (h) 73 Hình 19: Mực nước Tân Châu Để so sánh đánh giá kết tính từ mô hình số liệu quan trắc, số tiêu chí đánh giá độ thích hợp (GOF) kết tính từ mô hình giá trị quan trắc thực bảng 3: Q Cần Thơ Mỹ Thuận Châu Đốc Tân Châu Bảng 3: Một số tiêu so sánh GOF cho lưu lượng Q QTB (quan trắc) QTB (mô phỏng) Sai biệt đỉnh Sai biệt đáy (m3/s) (m3/s) TB (m3/s) TB (m3/s) 15851,0 14995,3 1832,4 1215,0 13225,6 13322,3 1068,4 65,7 6174,4 6169,3 11,2 90,2 22106,3 21886,3 290,7 16,1 Bảng 4: Một số tiêu so sánh GOF cho mực nước Z Z ZTB (quan trắc) ZTB (mô phỏng) (m) 1,07 0,92 3,29 3,81 (m) 0,99 0,88 3,22 3,42 Cần Thơ Mỹ Thuận Châu Đốc Tân Châu Sai biệt đỉnh TB (m) 0,19 0,09 0,10 0,42 Sai biệt đáy TB (m) 0,03 0,03 0,03 0,37 TÁC ĐỘNG PHÂN LŨ CỦA VÙNG ĐTM LÊN KHU VỰC ĐBSCL Mùa lũ ĐBSCL thường kéo dài khoảng tháng, từ tháng VII tháng XI hàng năm Nhìn chung, lũ ĐBSCL hiền hòa với biên độ Châu Đốc Tân Châu khoảng 3,5-4,0 m, lũ lên xuống chậm với cường suất trung bình 5-7 cm/ngày Trong trận lũ lớn sớm, cường suất đạt 10-12 cm/ngày cao từ 20-30 cm/ngày (lũ 1984) Lũ ĐBSCL thường có đỉnh, lớn vào khoảng cuối tháng IX đến nửa đầu tháng X Trong năm có hai đỉnh lũ, đỉnh phụ thường xuất sớm vào cuối tháng VIII đầu tháng IX, đỉnh thường cao có phần muộn trung bình nhiều năm Nhằm xem xét tác động phân lũ vào khu vực ĐTM lũ về, nghiên cứu xem xét phương án có vị trí thuộc đê bờ trái sông Tiền (hình 20) Hình 20: Vị trí hai phương án tràn PA không tràn PA tràn 1.8 PA tràn 1.5 Z (m) 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 Z (m) Chiều dài đỉnh tràn điển hình km Cao trình ngưỡng tràn 2,5 m Thời gian mô ngày mùa lũ điển hình từ 27/9/2013 đến 30/9/2013 Kết so sánh mực nước lưu lượng vị trí điển hình sông chính: Tân Châu, Cao Lãnh, Mỹ Thuận, Châu Đốc, Cần Thơ trình bày đồ thị bảng tóm tắt sau: Lũ ĐBSCL thường có đỉnh, lớn vào khoảng cuối tháng IX đến nửa đầu tháng X PA không tràn PA tràn PA tràn 1.2 0.9 0.6 0.3 T (h) 13 25 37 49 61 Hình 21:Mực nước Tân Châu 73 T (h) 0.0 13 25 37 49 61 Hình 22:Mực nước Cao Lãnh 73 0.9 0.6 0.3 0.0 1.2 13 25 37 49 61 Z (m) -0.3 1.5 0.9 73 T (h) Hình 23:Mực nước Mỹ Thuận 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 PA không tràn PA tràn Z (m) 1.2 PA không tràn PA tràn Z (m) 1.5 T (h) 13 25 37 49 61 73 Hình 24:Mực nước Châu Đốc 0.6 0.3 T (h) 0.0 -0.3 13 25 37 PA không tràn PA tràn 49 61 73 PA tràn Hình 25:Mực nước Cần Thơ Một số giá trị thay đổi điển hình mực nước đỉnh mực nước đáy tính trung bình cho ngày mô vị trí trạm đo nói tổng hợp bảng sau: Bảng 4: Thay đổi mực nước trung bình theo phương án tràn Mực nước trung bình (m) Không phân lũ PA tràn PA tràn Vị trí Tân Châu 3,83 2,44 3,67 Cao Lãnh 1,40 1,14 1,26 Mỹ Thuận 0,89 0,76 0,82 Châu Đốc 3,08 2,82 2,96 Cần Thơ 0,95 0,80 0,87 Bảng 5: Thay đổi mực nước lớn nhỏ trung bình theo phương án tràn Vị trí Tân Châu Cao Lãnh Mỹ Thuận Châu Đốc Cần Thơ Mực nước đỉnh trung bình (m) Không phân lũ PA tràn 3,82 1,50 1,23 3,08 1,33 2,44 1,32 1,15 2,82 1,24 Mực nước đáy trung bình (m) PA tràn Không phân lũ PA tràn PA tràn 3,67 1,40 1,19 2,96 1,29 3,81 0,93 -0,03 3,07 0,21 2,40 0,48 -0,26 2,81 -0,04 3,66 0,68 -0,16 2,95 0,08 Sự thay đổi lưu lượng hai vị trí điển hình phương án trình bày đồ thị sau: PA không tràn PA tràn 24000 PA tràn PA không tràn Q (m3/s) 16000 Q (m3/s) 20000 PA tràn 20000 PA tràn 16000 12000 12000 8000 8000 4000 4000 13 25 37 49 61 T (h) 73 Hình 26: Lưu lượng Mỹ Thuận 13 25 37 49 61 T 73 (h) Hình 27: Lưu lượng Cần Thơ Một số giá trị thay đổi điển hình lưu lượng cho ngày mô vị trí trạm đo nói tổng hợp bảng sau: Bảng 6: Thay đổi lưu lượng trung bình theo phương án tràn Lưu lượng trung bình m3/s) Không phân lũ PA tràn PA tràn Vị trí Mỹ Thuận 12394 10357 11271 Cần Thơ 16752 14326 15473 Bảng 7: Thay đổi lượng lớn nhỏ trung bình theo phương án tràn Vị trí Mỹ Thuận Cần Thơ Lưu lượng đỉnh trung bình m3/s Không phân lũ PA tràn 15353 21405 13277 19471 Lưu lượng đáy trung bình m3/s PA tràn Không phân lũ PA tràn PA tràn 14200 20375 5716 6419 1615 1133 3560 3699 Kết tính toán trình bày nêu cho phép rút số nhận xét sau đây: - Khả phân lũ ĐTM làm cho giảm mực nước vùng hạ lưu rõ ràng Với phương vị trí tràn 1, mực nước trung bình Mỹ Thuận giảm 0,12 m vị trí tràn giảm 0,07 m - Do tăng khả thoát nước có phương án tràn, mực nước trung bình Tân Châu giảm nhiều (1,39 m với phương án tràn 1) so với trường hợp phương án tràn Kết cho thấy tác động yếu tố tăng khả thoát lũ nhanh có ảnh hưởng lớn đến việc giảm mực nước toàn vùng ĐBSCL - Phương án vị trí tràn tràn tính với cao trình ngưỡng tràn (2,5 m) chiều rộng diện tràn (2 km), nhiên khả thoát lưu lượng phân lũ qua vùng ĐTM phương án tràn lớn mực nước sông cao nên cột nước tràn qua đỉnh lớn - Tác động làm hạ mực nước lưu lượng giảm dần tiến gần hạ lưu ảnh hưởng mạnh mực nước biển lên yếu tố gần biên hạ lưu KẾT LUẬN Phương án phân lũ vào vùng ĐTM khả xem xét nhằm mục đích giảm mực nước, lưu lượng cho số vị trí hệ thống sông Mékong Với lựa chọn thiết kế hợp lý cho cao trình đê phân lũ chiều dài tràn vị trí thích hợp cho phép thực chủ trương thiết kế công trình thủy lợi theo hướng “sống chung với lũ”, thay xây dựng hệ thống đê bao vượt lũ triệt để Vì với hệ thống đê bao triệt để làm giảm khả thoát lũ bình diện toàn ĐBSCL tăng áp lực lên dòng sông Tác động làm gia tăng nguy không an toàn cho vùng đô thị quan trọng tuyến sông TÀI LIỆU THAM KHẢO Viện nghiên cứu thủy lợi Miền Nam (2009): Quy hoạch kiểm soát lũ Đồng sông Cửu Long Viện nghiên cứu thủy lợi Miền Nam (2009): Quy hoạch đê biển-đê cửa sông Đồng sông Cửu Long HERVOUET Jean Michel (2007) Hydrodynamics of Free Surface Flows modelling with the finite element method WILEY LANG Pierre and all (2010) Telemac2d_manuel_utilisateur_v6p0 EDF MEISSNER Loren P (1995) Fortran 90 PWS Publishing Company NOAA National Geophysical Center http://www.ngdc.noaa.gov/mgg/global/global.html OTIS Regional Tidal Solutions http://volkov.oce.orst.edu/tides/region.html Lời cám ơn: - Nghiên cứu thực khuôn khổ đề tài đặt hàng Trung tâm CARE-RESCIF, HCMUT-VNUHCM Mã số T-KTXD-2014-86 - Các mô thực hệ thống máy tính hiệu cao (HPCC), trường Đại học Bách Khoa Tp HCM Người phản biện: TS Trương Ngọc Tường Ngày nhận bài: 25/7/2015 Ngày phản biện xong: 31/08/2015 [...].. .sông chính Tác động này làm gia tăng nguy cơ không an toàn cho các vùng đô thị quan trọng trên các tuyến sông này TÀI LIỆU THAM KHẢO 1 Viện nghiên cứu thủy lợi Miền Nam (2009): Quy hoạch kiểm soát lũ Đồng bằng sông Cửu Long 2 Viện nghiên cứu thủy lợi Miền Nam (2009): Quy hoạch đê biển-đê cửa sông Đồng bằng sông Cửu Long 3 HERVOUET Jean Michel (2007) Hydrodynamics... http://www.ngdc.noaa.gov/mgg/global/global.html 7 OTIS Regional Tidal Solutions http://volkov.oce.orst.edu/tides/region.html Lời cám ơn: - Nghiên cứu này được thực hiện trong khuôn khổ đề tài đặt hàng của Trung tâm CARE-RESCIF, HCMUT-VNUHCM Mã số T-KTXD-2014-86 - Các mô phỏng đã được thực hiện trên hệ thống máy tính hiệu năng cao (HPCC), trường Đại học Bách Khoa Tp HCM Người phản biện: TS Trương Ngọc