Ứng dụng mô hình thủy động lực 2 chiều mô phỏng xâm nhập mặn tại sông tiền và sông hậu thuộc đồng bằng sông cửu long

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THUỶ LỢI - NGUYỄN THỊ THU ỨNG DỤNG MÔ HÌNH THỦY ĐỘNG LỰC CHIỀU MƠ PHỎNG XÂM NHẬP MẶN TẠI SÔNG TIỀN VÀ SÔNG HẬU THUỘC ĐỒNG BẰNG SÔNG CỬU LONG LUẬN VĂN THẠC SĨ HÀ NỘI, NĂM 2021 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THUỶ LỢI - NGUYỄN THỊ THU ỨNG DỤNG MƠ HÌNH THỦY ĐỘNG LỰC CHIỀU MƠ PHỎNG XÂM NHẬP MẶN TẠI SÔNG TIỀN VÀ SÔNG HẬU THUỘC ĐỒNG BẰNG SƠNG CỬU LONG Chun ngành: Khoa học mơi trường Mã số: 8440301 NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Phạm Văn Chiến PGS.TS Bùi Quốc Lập HÀ NỘI, NĂM 2021 LỜI CAM ĐOAN Tác giả xin cam đoan cơng trình nghiên cứu thân tác giả Đề tài nghiên cứu không trùng lặp với đề tài luận văn trước Các kết nghiên cứu kết luận luận văn trung thực, không chép từ nguồn hình thức Việc tham khảo nguồn tài liệu (nếu có) thực trích dẫn ghi nguồn tài liệu tham khảo quy định Tác giả luận văn Nguyễn Thị Thu i LỜI CẢM ƠN Trước hết, học viên xin trân trọng cảm ơn PGS.TS Phạm Văn Chiến – Khoa Kỹ thuật tài nguyên nước PGS.TS Bùi Quốc Lập Khoa Hóa Mơi trường, Trường Đại học Thuỷ lợi tận tình hướng dẫn, định hướng tạo điều kiện cho em hoàn thành luận văn Em xin chân thành cảm ơn thầy, Khoa Hóa Môi trường, Khoa Kỹ thuật tài nguyên nước, phòng Đào tạo Đại học Sau đại học, Trường Đại học Thuỷ lợi động viên, khích lệ đóng góp ý kiến quý báu cho em việc soạn thảo, hướng dẫn thủ tục để em hoàn thành luận văn thuận lợi Trong trình thực hoàn thành luận văn, thời gian kiến thức còn hạn chế nên tránh khỏi thiếu sót Em mong nhận đóng góp ý kiến, bảo tận tình q thầy, bạn để luận văn hoàn thiện Em xin chân thành cảm ơn! Hà nội, ngày tháng năm 2021 Tác giả luận văn Nguyễn Thị Thu ii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC BẢNG v DANH MỤC CHỮ VIẾT TẮT vi DANH MỤC HÌNH vii MỞ ĐẦU CHƯƠNG TỔNG QUAN VỀ XÂM NHẬP MẶN VÀ KHU VỰC NGHIÊN CỨU 1.1 Tổng quan xâm nhập mặn vùng châu thổ lớn Việt Nam 1.1.1 1.1.2 1.1.3 1.1.4 Khái niệm xâm nhập mặn Tổng quan xâm nhập mặn giới Tổng quan xâm nhập mặn Việt Nam Tổng quan xâm nhập mặn đồng sông Cửu Long 10 1.2 Tổng quan vùng ĐBCSL 13 1.2.1 1.2.2 1.2.3 1.2.4 1.2.5 Vị trí địa lý 13 Địa hình-địa mạo 14 Thổ nhưỡng, đất đai 16 Đặc điểm thời tiết khí hậu .16 Đặc điểm thủy văn tài nguyên nước 19 1.3 Hiện trạng xâm nhập mặn vùng ĐBSCL .22 1.4 Ảnh hường xâm nhập mặn vùng ĐBSCL .24 1.4.1 1.4.2 Ảnh hưởng mặn đến hoạt động trồng trọt, nuôi trồng thủy hải sản 24 Ảnh hưởng hạn mặn đến hoạt động dân sinh, kinh tế, xã hội 26 CHƯƠNG PHƯƠNG PHÁP NGHIÊN CỨU VÀ THIẾT LẬP MƠ HÌNH MƠ PHỎNG XÂM NHẬP MẶN 28 2.1 Phân tích lựa chọn mơ hình tốn 28 2.1.1 2.1.2 2.1.3 Một số mơ hình mơ xâm nhập mặn giới 28 Một số mơ hình mơ xâm nhập mặn Việt Nam 31 Lựa chọn mơ hình sử dụng .32 2.2 Giới thiệu mơ hình tốn 33 iii 2.2.1 2.2.2 Giới thiệu module thủy động lực 33 Giới thiệu module lan truyền mặn 35 2.3 Các tiêu đánh giá sai số 38 2.4 Xây dựng mơ hình tính toán 38 2.4.1 2.4.2 2.4.3 2.4.4 Thiết lập mơ hình mô thủy động lực 39 Kết mô thủy động lực 42 Thiết lập mơ hình mơ lan truyền xâm nhập mặn 53 Kết mô xâm nhập mặn 57 CHƯƠNG KẾT QUẢ MÔ PHỎNG MẶN VÀ ĐỀ XUẤT CÁC BIỆN PHÁP GIẢM THIỂU 64 3.1 Kết mô lan truyền mặn theo không gian vùng nghiên cứu 64 3.2 Xác định kịch mơ lan truyền xâm nhập mặn có xét đến ảnh hưởng BĐKH 68 3.3 Mô xâm nhập mặn theo kịch 72 3.3.1 3.3.2 Kết mô theo kịch BĐKH 72 So sánh kết mô xâm nhập mặn với nghiên cứu trứớc 74 3.4 Đề xuất giải pháp giảm thiểu, thích ứng 76 3.4.1 3.4.2 Giải pháp phi cơng trình 79 Giải pháp cơng trình 83 KẾT LUẬN VÀ KIẾN NGHỊ 89 TÀI LIỆU THAM KHẢO 92 PHỤ LỤC: KẾT QUẢ MÔ PHỎNG PHÂN BỐ VẬN TỐC THEO HƯỚNG DÒNG CHẢY 99 iv DANH MỤC BẢNG Bảng 1.1: Nhiệt độ khơng khí trung bình (oC) [31] 18 Bảng 1.2: Tổng lượng mưa trung bình (mm) [31] 18 Bảng 2.1: Bảng thống kê vị trí biên thượng hạ lưu vùng tính tốn 42 Bảng 2.2: Bảng tổng hợp tiêu sai số cho mực nước số vị trí 45 Bảng 2.3: Bảng tổng hợp tiêu sai số kiểm định với n = 0.035 52 Bảng 2.4: Bảng thống kê vị trí biên thượng hạ lưu vùng tính toán 54 Bảng 3.1: Các tiêu sai số cho hiệu chỉnh thống sô module lan truyền mặn ứng với hệ số khuyếch tán K= 300 m2/s .58 Bảng 3.2: Các tiêu sai số cho kiểm định mô hình lan truyền mặn 62 Bảng 3.3: Kết mô chiều dài đoạn sông dài với nồng độ mặn 4g/l mùa khô năm (tháng 4) 2015 2016 67 Bảng 3.4: Kết mô chiều dài xâm nhập mặn nhóm tác giả Viện khoa học Thủy lợi miền Nam [40] 68 Bảng 2.5: Mực nước biển dâng theo kich RCP4.5 khu vực Mũi Kê Gà – Cà Mau72 Bảng 2.6: Các kịch mô 72 Bảng 3.5: Diễn biến xâm nhập mặn mùa khô qua năm 73 Bảng 3.6: Diện tích xâm nhập mặn năm tăng so với năm 2015 (%) 75 Bảng 3.7: Khả ni số lồi thủy sản theo vùng nhiễm mặn [45] [46] 81 Bảng 3.8: Bảng tổng hợp giai đoạn triển khai cống đề xuất tỉnh [44] 84 v DANH MỤC CHỮ VIẾT TẮT STT Từ viết tắt Nội dung BĐCM Bán đảo Cà Mau ĐBSCL Đồng sông Cửu Long ĐTM Đồng Tháp Mười HSH Hữu sông Hậu NTTS Nuôi trồng thuỷ sản SXXNN Sản xuất nơng nghiệp TBNN Trung bình nhiều năm TGHT Tứ giác Hà Tiên TGLX Tứ giác Long Xuyên 10 XNM Xâm nhập mặn vi DANH MỤC HÌNH Hình 1.1: Khu vực nghiên cứu vị trí trạm thủy văn 14 Hình 1.2: Dạng triều biển Đơng 21 Hình 1.3: Dạng triều biển Tây 22 Hình 1.4: Độ mặn lớn đầu tháng 3/2016 so với kỳ năm 2015 vùng cửa sông Cửu Long [32] .24 Hình 2.1: Lưới tính tốn 37 Hình 2.1: Giới hạn vùng tính tốn mơ hình 39 Hình 2.2: Địa hình vùng tính tốn .40 Hình 2.3: Lưới dung tính tốn .41 Hình 2.4: Đường trình mực nước thực đo trạm Vàm Nao 43 Hình 2.5: Đường trình mực nước thực đo trạm Mỹ Thuận 43 Hình 2.6: Đường trình mực nước thực đo trạm Cần Thơ 44 Hình 2.11: Mực nước tính toán thực đo trạm Vàm Nao .45 Hình 2.12: Mực nước tính tốn thực đo trạm Cần Thơ 46 Hình 2.13: Mực nước tính tốn thực đo trạm Mỹ Thuận 46 Hình 2.4: Đường trình lưu lượng ngày trạm thủy văn Châu Đốc 48 Hình 2.5: Đường trình lưu lượng ngày trạm thủy văn Tân Châu .49 Hình 2.6: Đường trình mực nước trạm Vàm Kênh 49 Hình 2.7: Đường trình mực nước trạm Bình Đại .50 Hình 2.8: Đường trình mực nước trạm An Thuận 50 Hình 2.9: Đường trình mực nước trạm Bến Trại 51 Hình 2.10: Đường trình mực nước trạm Mỹ Thanh 51 Hình 2.14: Mực nước tính tốn thực đo trạm Vàm Nao .53 Hình 2.15: Mực nước tính tốn thực đo trạm Mỹ Thuận 53 vii Hình 2.16: Mực nước tính toán thực đo trạm Cần Thơ 53 Hình 2.17:Giá trị đo mặn trạm Bình Đại năm 2015 54 Hình 2.18:Giá trị đo mặn trạm Bến Trại năm 2015 55 Hình 2.19:Giá trị đo mặn trạm An Định năm 2015 56 Hình 2.20:Giá trị đo mặn trạm Bình Đại năm 2016 56 Hình 2.21:Giá trị đo mặn trạm Bến Trại năm 2016 57 Hình 2.22:Giá trị đo mặn trạm An Định năm 2016 57 Hình 3.1: Độ mặn tính tốn thực đo trạm Hồ Bình 59 Hình 3.2: Độ mặn tính tốn thực đo trạm Đại Ngải 60 Hình 3.3: Độ mặn tính tốn thực đo trạm Hồ Bình 61 Hình 3.4: Độ mặn tính tốn thực đo trạm Đại Ngải 61 Hình 3.5: Độ mặn thời điểm 5h ngày 19/3/2016 64 Hình 3.6: Độ mặn thời điểm 7h ngày 19/3/2016 65 Hình 3.7: Độ mặn thời điểm 9h ngày 19/3/2016 65 Hình 3.8: Độ mặn thời điểm 13h ngày 17/4/2016 66 Hình 3.9: Độ mặn thời điểm 15h ngày 17/4/2016 66 Hình 3.10: Độ mặn thời điểm 19h ngày 17/4/2016 67 Hình 2.23: Xu biến đổi mực nước biển từ số liệu vệ tinh Biển Đơng [42]) 71 Hình 3.11: Kết mô XNM theo kịch PRC4.5 73 Hình 3.12: Diễn biến xâm nhập mặn mùa khô 2015-2016 [31] 75 Hình 3.13: Các đường đẳng mặn tháng khô hạn (tháng 4, bên trái) tháng bắt đầu mừa mưa (tháng 6, bên phải) mô năm 2050 với mực nước biển dâng 30cm [44] 76 Hình 3.14: Dự báo xâm nhập mặn ĐBSCL năm 2050 80 Hình 3.15: Kế hoạch xây dựng cửa cống năm 2050[44] 85 viii sông nghiên cứu, nước sông Tiền sông Hậu dọc theo hai nhánh có chuyển nước qua lại nhờ hệ thống sông nhánh dầy đặc tác động mạnh mẽ thủy triều cửa sơng sâu vào sơng nội địa, chế động thủy động lực dòng chảy kiệt hệ thống sông phức tạp Bên cạnh, hạn chế thời gian, số liệu thực đo nên kết tính tốn mơ cịn khác biệt với giá trị đo đạc Độ mặn đồng kết cân dòng nước mặn thủy triều điều khiển xâm nhập vào đồng thông qua trộn lẫn dòng nước đẩy nước mặn trở lại Cả hai dịng chảy phụ thuộc mạnh mẽ vào địa hình, độ sâu lịng sơng, lưu lượng nước điều kiện thủy hải văn Tuy nhiên, không đầy đủ đồng liệu nêu nguyên nhân dẫn đến khác biệt kết tính tốn từ mơ hình so với giá trị thực đo Do đó, để kết mô tốt cần cập nhật liệu địa số liệu thủy hải văn liên quan khác 91 TÀI LIỆU THAM KHẢO [1] Nghị số 120/NQ-CP, “Nghị số 120/NQ - CP Chính phủ : Về phát triển bền vững đồng sơng Cửu Long thích ứng với biến đổi khí hậu,” 2017 [2] Trung tâm quy hoạch điều tra tài nguyên nước quốc gia, “Xâm nhập mặn đồng sông Cửu Long giải pháp khắc phục,” 2019 http://nawapi.gov.vn/index.php?option=com_content&view=article&id=4108%3 Axam-nhp-mn-ng-bng-song-cu-long-va-cac-gii-phap-khcphc&catid=70%3Anhim-v-chuyen-mon-ang-thc-hin&Itemid=135&lang=vi (accessed Jul 17, 2020) [3] Phung L M., Trí V P Đ., and Đạt T Q., “ỨNG DỤNG MƠ HÌNH TỐN THỦY LỰC MỘT CHIỀU ĐÁNH GIÁ VÀ DỰ BÁO TÌNH HÌNH XÂM NHẬP MẶN TRÊN HỆ THỐNG SƠNG CHÍNH TRÊN ĐỊA BÀN TỈNH TRÀ VINH,” Tạp Chí Khoa Học Trường Đại Học Cần Thơ, vol Số 25, pp 68–75, May 2013 [4] Lê Anh Tuấn, “Giáo trình Thủy văn mơi trường.” Đại học Cần Thơ, 2008 [Online] Available: http://www.leanhtuan.com/pdf/GT_ThuyVanMoiTruong.pdf [5] Priyantha Ranjan, “Effect of climate change and land use change on saltwater intrusion - The Encyclopedia of Earth,” Jul 11, 2012 https://editors.eol.org/eoearth/wiki/Effect_of_climate_change_and_land_use_cha nge_on_saltwater_intrusion (accessed May 22, 2021) [6] DMC, “Kiến thức về: Xâm nhập mặn.” Trung Tâm Phòng tránh giảm nhẹ thiên tai, 2016 [7] W P Budgell, “A Dynamic-Stochastic Approach for Modelling AdvectionDispersion Processes in Open Channels,” in Developments in Water Science, vol 17, A H El-Shaarawi and S R Esterby, Eds Elsevier, 1982, pp 244–263 doi: 10.1016/S0167-5648(08)70714-7 [8] D W Pritchard, “The Dymanic structure of a coastal plain estuary,” J Mar Res., vol 15, pp 33–42, 1971 92 [9] J J Leendertse, “Aspects of a Computational Model for Long-Period Water-Wave Propagation,” Dec 1966, Accessed: Jan 23, 2021 [Online] Available: https://www.rand.org/pubs/research_memoranda/RM5294.html [10] R C Martyr-Koller et al., “Application of an unstructured 3D finite volume numerical model to flows and salinity dynamics in the San Francisco Bay-Delta,” Estuar Coast Shelf Sci., vol 192, pp 86–107, Jun 2017, doi: 10.1016/j.ecss.2017.04.024 [11] M S Maleki Tirabadi, M E Banihabib, and T O Randhir, “SWAT-S: A SWATsalinity module for watershed-scale modeling of natural salinity,” Environ Model Softw., vol 135, p 104906, Jan 2021, doi: 10.1016/j.envsoft.2020.104906 [12] A Eladawy, A Negm, M Elzeir, O Saavedra, I El-Shinnawy, and K Nadaoka, “Modeling the Hydrodynamics and Salinity of El-Burullus Lake (Nile Delta, Northern Egypt),” J Clean Energy Technol., vol 1, pp 157–163, Jan 2013, doi: 10.7763/JOCET.2013.V1.37 [13] W Zhang et al., “Numerical Simulation and Analysis of Saltwater Intrusion Lengths in the Pearl River Delta, China,” J Coast Res., vol 29, no 2, pp 372– 382, Sep 2012, doi: 10.2112/JCOASTRES-D-12-00068.1 [14] D J Twigt, E D De Goede, F Zijl, D Schwanenberg, and A Y W Chiu, “Coupled 1D–3D hydrodynamic modelling, with application to the Pearl River Delta,” Ocean Dyn., vol 59, no 6, p 1077, Sep 2009, doi: 10.1007/s10236-0090229-y [15] P Rengasamy, “World salinization with emphasis on Australia,” J Exp Bot., vol 57, no 5, pp 1017–1023, Mar 2006, doi: 10.1093/jxb/erj108 [16] S S Kaushal et al., “Increased salinization of fresh water in the northeastern United States,” Proc Natl Acad Sci., vol 102, no 38, pp 13517–13520, Sep 2005, doi: 10.1073/pnas.0506414102 [17] H V Cường and Nguyễn Thị Ngọc Nhẫn, “Nghiên cứu ảnh hưởng phân cấp 93 yếu tố tác động đến xâm nhập mặn vùng hạ du sơng Hồng-Thái Bình,” Tạp Chí Nơng Nghiệp Và Phát Triển Nông Thôn, vol Kỳ 1-Tháng 4, no 382, pp 83–95, 2020 [18] Hạnh N Đ., Hải Ô T., Hạ L Á., Anh T N., and Dư N H., “NGHIÊN CỨU XÂY DỰNG CHƯƠNG TRÌNH MƠ PHỎNG Q TRÌNH XÂM NHẬP MẶN HAI CHIỀU (2D),” p 10 [19] Tấn T T., “NGHIÊN CỨU TÁC ĐỘNG CỦA XÂM NHẬP MẶN ĐẾN ĐẤT TRỒNG LÚA TẠI XÃ PHÚ AN HUYỆN PHÚ VANG, TỈNH THỪA THIÊN HUẾ,” Hue Univ J Sci HU JOS, vol 126, no 7A, Art no 7A, Mar 2017, Accessed: May 22, 2021 [Online] Available: http://222.255.146.83/index.php/TCKHDHH/article/view/4042 [20] Đoàn Văn Hải and Lê Thị Huệ, “Nghiên cứu ứng dụng mô hình hóa xây dựng phần mềm dự báo lũ, xâm nhập mặn sông Cửu Long hiển thị kết dự báo mặn lên Google Earth.” http://tapchikttv.vn/article/492 (accessed Jan 16, 2021) [21] T Q Toan, “9 - Climate Change and Sea Level Rise in the Mekong Delta: Flood, Tidal Inundation, Salinity Intrusion, and Irrigation Adaptation Methods,” in Coastal Disasters and Climate Change in Vietnam, N D Thao, H Takagi, and M Esteban, Eds Oxford: Elsevier, 2014, pp 199–218 doi: 10.1016/B978-0-12800007-6.00009-5 [22] S Eslami et al., “Tidal amplification and salt intrusion in the Mekong Delta driven by anthropogenic sediment starvation,” Sci Rep., vol 9, no 1, Art no 1, Dec 2019, doi: 10.1038/s41598-019-55018-9 [23] D T Vu, T Yamada, and H Ishidaira, “Assessing the impact of sea level rise due to climate change on seawater intrusion in Mekong Delta, Vietnam,” Water Sci Technol., vol 77, no 6, pp 1632–1639, Jan 2018, doi: 10.2166/wst.2018.038 [24] T V Tran et al., “Examining spatiotemporal salinity dynamics in the Mekong River Delta using Landsat time series imagery and a spatial regression approach,” Sci Total Environ., vol 687, pp 94 1087–1097, Oct 2019, doi: 10.1016/j.scitotenv.2019.06.056 [25] P A Duc, N T M Linh, and P Van Mien, “Study on Migration of Marine Organisms into Inland and Changes of Biodiversity at Water Bodies in Mekong Delta for Evaluation of Saline Intrusion of Sea Level Rise,” APCBEE Procedia, vol 1, pp 252–257, Jan 2012, doi: 10.1016/j.apcbee.2012.03.041 [26] D Tran Anh, L P Hoang, M D Bui, and P Rutschmann, “Simulating Future Flows and Salinity Intrusion Using Combined One- and Two-Dimensional Hydrodynamic Modelling—The Case of Hau River, Vietnamese Mekong Delta,” Water, vol 10, no 7, Art no 7, Jul 2018, doi: 10.3390/w10070897 [27] R Wassmann et al., “High-resolution mapping of flood and salinity risks for rice production in the Vietnamese Mekong Delta,” Field Crops Res., vol 236, pp 111– 120, Apr 2019, doi: 10.1016/j.fcr.2019.03.007 [28] Vũ Thùy Linh, N D Liêm, Hồ Minh Dũng, and Nguyễn Kim Lợi, “Nghiên cứu ứng dụng mơ hình hóa đánh giá xu ngập lụt xâm nhập mặn bối cảnh biến đổi khí hậu: nghiên cứu thí điểm thành phố Hồ Chí Minh,” Tạp Chí Khí Tượng Thủ Văn, vol EME2, pp 98–110, 2019 [29] Trần Quốc Đạt, Nguyễn Hiếu Trung, and Kanchit Likitdecharote, “MÔ PHỎNG XÂM NHẬP MẶN ĐỒNG BẰNG SÔNG CỬU LONG DƯỚI TÁC ĐỘNG MỰC NƯỚC BIỂN DÂNG VÀ SỰ SUY GIẢM LƯU LƯỢNG TỪ THƯỢNG NGUỒN,” Tạp Chí Khoa Học, vol 21b, no Đại học Cần Thơ, pp 141–150, 2012 [30] N D Khang, A Kotera, T Sakamoto, and M Yokozawa, “Sensitivity of Salinity Intrusion to Sea Level Rise and River Flow Change in Vietnamese Mekong DeltaImpacts on Availability of Irrigation Water for Rice Cropping,” J Agric Meteorol., vol 64, no 3, pp 167–176, 2008, doi: 10.2480/agrmet.64.3.4 [31] Viện Quy hoạch Thuỷ lợi miền Nam, “Quy hoạch tổng thể thuỷ lợi đồng sông Cửu Long điều kiện biến đổi khí hậu-nước biển dâng.” Viện Quy hoạch Thủy lợi miền Nam, 2017 95 [32] VAWR, “Báo cáo xâm nhập mặn cửa sông vùng ven biển đồng sông Cửu Long đề xuất giải pháp chống hạn.” Viện Khoa học Thủy lợi Việt Nam, 2016 [33] Nguyễn Ân Niên and Nguyễn Văn Lân, “Nghiên cứu xâm nhập mặn Việt Nam.” Viện Khoa học Thủy lợi Việt Nam, 1999 [34] Lê Hữu Thuần, “Nghiên cứu sở khoa học xác định nguyên nhân, đề xuất giải pháp ứng phó với xâm nhập mặn điều kiện Biến đổi khí hậu vùng ĐBSCL.” Cục Quản lý tài nguyên nước, 2013 [35] Trần Thị Tuyết, “Nghiên cứu đề xuất khung sách chủ động thích ứng với biến đổi khí hậu vùng đồng sông Cửu Long từ tiếp cận sinh thái xã hội.” Đề tài khoa học mã số 507.01-2019.301, Quỹ Phát triển Khoa học Công nghệ quốc gia tài trợ., 2019 [36] Hà Huy Ngọc and Trần Thị Tuyết, “Kiến tạo sách nhằm chủ động thích ứng với biến đổi khí hậu vùng đồng sơng Cửu Long,” 2020 https://www.mof.gov.vn/webcenter/portal/vclvcstc/r/m/nckh/ctnc/nckhctnc_chitie t?dDocName=MOFUCM183108&dID=191132&_afrLoop=2523847638204213# %40%3FdID%3D191132%26_afrLoop%3D2523847638204213%26dDocName %3DMOFUCM183108%26_adf.ctrl-state%3Dgc0l4azno_4 (accessed Jan 18, 2021) [37] Nguy K., “MỘT MƠ HÌNH BÀI TỐN HAI CHIỀU CHO DÒNG CHẢY TRÊN CÁC LÒNG DẪN HỞ BỀ MẶT TỰ DO,” p [38] Hậu H P., “Luận án Tiến sĩ ‘Mơ hình tốn học dịng chảy hở chiều suy rộng,’” Đại Học Đà Nẵng, p 193, 2019 [39] C Pham Van, E Deleersnijder, D Bousmar, and S Soares-Frazão, “Simulation of flow in compound open-channel using a discontinuous Galerkin finite-element method with Smagorinsky turbulence closure,” J Hydro-Environ Res., vol 8, no 4, pp 396–409, Dec 2014, doi: 10.1016/j.jher.2014.04.002 [40] Tăng Đức Thắng, Tô Quang Toản, and Nguyễn Quang Kim, “Một số vấn đề 96 ĐBSCL điều kiện phát triển thượng lưu, biến đổi khí hậu-nước biển dâng, lún sụt đất bão tố.” Viện Quy hoạch Thủy lợi miền Nam, 2016 [41] Bộ Tài nguyên môi trường, “Kịch Biến đổi khí hậu nước biển dâng cho Việt Nam.” 2016 [42] Bộ Tài nguyên Mơi trường, “Dự thảo Kịch Biến đổi khí hậu nước biển dâng cho Việt Nam.” 2020 [43] C T Hoanh, K Jirayoot, G Lacombe, and V Srinetr, “Impacts of climate change and development on Mekong flow regimes First assessment - 2009,” International Water Management Institute, IWMI Research Reports H043262, 2010 Accessed: May 01, 2021 [Online] Available: https://econpapers.repec.org/paper/iwtrerpts/h043262.htm [44] Viện Quy hoạch Thuỷ lợi miền Nam, Bộ Nông nghiệp Phát triển nơng thơn, and JICA, “DỰ ÁN THÍCH ỨNG VỚI BIẾN ĐỔI KHÍ HẬU CHO PHÁT TRIỂN BỀN VỮNG NÔNG NGHIỆP VÀ NÔNG THÔN VÙNG VEN BIỂN ĐỒNG BẰNG SƠNG CỬU LONG.” JICA, 2013 [45] Đồn Xn Diệp, “Ảnh hưởng độ mặn oxy hòa tan lên thay đổi sinh lý sinh trưởng tôm sú (Penaeus monodon Fabricius, 1798).” Luận án Tiến sĩ , Đại học Cần Thơ, 2012 [46] Lê Thị Phương Mai, “Nghiên cứu tác động xâm nhập mặn khả thích ứng nuôi trồng thủy hải sản Đồng sông Cửu Long.” Luận án Tiến sĩ , Đại học Cần Thơ, 2017 [47] Lam Mỹ Lan, Trần Ngọc Thảo, and Đỗ Thị Thanh Hương, “Ảnh hưởng độ mặn lên điều hòa áp suất thẩm thấu tăng trưởng cá leo (Wallago attu).” Tạp chí khoa học Đại học Cần Thơ, Số Chuyên đề thuỷ sản Số 1, trang 319-325, 2011 [48] Nguyễn Thanh Phương, Nguyễn Anh Tuấn, Trần Ngọc Hải, Võ Nam Sơn, and Dương Nhựt Long, Giáo trình Ni trồng Thuỷ sản, vol 188 trang Nhà xuất Đại học Cần Thơ, 2014 97 [49] Nguyễn Trọng Hồng Phúc, Trần Thị Kiều Linh, Nguyễn Minh Trí, Thị Thê Phước, Trần Thanh Trang, and Nguyễn Thanh Phương, “Ảnh hưởng tương tác nhiệt độ độ mặn lên tăng trưởng hormone tăng trưởng cá tra giống.,” Tạp Chí Khoa Học Trường Đại Học Cần Thơ, vol Số 36, pp 88–97, 2015 98 PHỤ LỤC: KẾT QUẢ MÔ PHỎNG PHÂN BỐ VẬN TỐC THEO HƯỚNG DỊNG CHẢY Kết mơ vận tốc dịng chảy theo hướng dọc Từ Hình 0.1 đến Hình 0.6 ta thấy phân bố vận tốc hướng dọc hai trạm Mỹ Tho Trà Vinh Tại thời điểm 5h ngày 19/3/2016 vận tốc hướng dọc trạm Mỹ Tho dao động từ -0.36 – -0.15 m/s đến 7h vận tốc hướng dọc nằm khoảng 0.15 – 0.25 m/s Thời điểm 9h ngày 19/3/2016 vận tốc nằm khoảng 0.2 – 0.4 m/s Cùng ngày 19/3/2016, trạm Trà Vinh có vận tốc hướng dọc dao động từ -0.2– -0.1 m/s lúc 5h, đến 7h 9h 0.25 – 0.35 m/s 0.16 – 0.3 m/s Ngày 17/4/2016, vận tốc hướng dọc trạm Mỹ Tho khoảng từ 0.16 – 0.48 m/s lúc 1h dao động khoảng -0.12 – -0.04 m/s lúc 3h Tại lúc 7h ngày 17/4/2016, vận tốc khoảng 0.1 – 0.2 m/s Cũng ngày 17/4/2016, vận tốc hướng dọc trạm Trà Vinh có thay đổi Ở trạm Trà Vinh, vận tốc nằm khoảng 0.2 – 0.5 m/s lúc 1h, giảm xuống 0.08 – 0.16 m/s lúc 3h lúc 7h cịn từ -0.06 – 0.02 m/s Hình 0.1: Trường vận tốc hướng dọc (thành phần 𝑢 ̅ ) thời điểm 5h ngày 19/3/2016 99 Hình 0.2: Trường vận tốc hướng dọc (thành phần 𝑢 ̅ ) thời điểm 7h ngày 19/3/2016 Hình 0.3: Trường vận tốc hướng dọc (thành phần 𝑢 ̅ ) thời điểm 9h ngày 19/3/2016 100 Hình 0.4: Trường vận tốc hướng dọc (thành phần 𝑢 ̅ ) thời điểm 1h ngày 17/4/2016 Hình 0.5: Trường vận tốc hướng dọc (thành phần 𝑢 ̅ ) thời điểm 3h ngày 17/4/2016 101 Hình 0.6: Trường vận tốc hướng dọc (thành phần 𝑢 ̅ ) thời điểm 7h ngày 17/4/2016 Kết mơ vận tốc dịng chảy theo hướng ngang Kết mơ phân bố vận tốc dịng chảy theo hướng ngang vùng nghiên cứu thể hình từ Hình 0.7 đến Hình 0.12 Tại thời điểm 5h ngày 19/3/2016 vận tốc hướng ngang trạm Mỹ Tho dao động từ 0.1 – 0.5 m/s đến 7h vận tốc hướng ngang nằm khoảng 0.08 – 0.5 m/s Thời điểm 9h ngày 19/3/2016 vận tốc nằm khoảng 0.02 – 0.45 m/s Cùng ngày 19/3/2016, trạm Trà Vinh có vận tốc hướng ngang dao động từ 0.2 – 0.5 m/s lúc 5h, đến 7h 9h 0.1 – 0.35 m/s -0.15 – 0.3 m/s Ngày 17/4/2016, vận tốc hướng ngang trạm Mỹ Tho khoảng từ 0.06 – 0.4 m/s lúc 1h dao động khoảng -0.16 – 0.32 m/s lúc 3h Tại lúc 7h ngày 17/4/2016, vận tốc khoảng -0.02 – 0.06 m/s Cũng ngày 17/4/2016, vận tốc hướng ngang trạm Trà Vinh có thay đổi Ở trạm Trà Vinh, vận tốc nằm khoảng 0.01 – 0.3 m/s lúc 1h, giảm xuống 0.08 – 0.25 m/s lúc 3h lúc 7h cịn từ -0.04 – 0.08 m/s 102 Hình 0.7: Trường vận tốc hướng ngang (thành phần 𝑣̅ ) thời điểm 5h ngày 19/3/2016 Hình 0.8: Trường vận tốc hướng ngang (thành phần 𝑣̅ ) thời điểm 7h ngày 19/3/2016 103 Hình 0.9: Trường vận tốc hướng ngang (thành phần 𝑣̅ ) thời điểm 9h ngày 19/3/2016 Hình 0.10: Trường vận tốc hướng ngang (thành phần 𝑣̅ ) thời điểm 1h ngày 17/4/2016 104 Hình 0.11: Trường vận tốc hướng ngang (thành phần 𝑣̅ ) thời điểm 3h ngày 17/4/2016 Hình 0.12: Trường vận tốc hướng ngang (thành phần 𝑣̅ ) thời điểm 7h ngày 17/4/2016 105