NGHIÊN CỨU ĐÁNH GIÁ TÁC ĐỘNG XÂM NHẬP MẶN CHO HỆ THỐNG KÊNH MƯƠNG NỘI ĐỒNG

BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI LUẬN VĂN THẠC SĨ NGHIÊN CỨU ĐÁNH GIÁ TÁC ĐỘNG XÂM NHẬP MẶN CHO HỆ THỐNG KÊNH MƯƠNG NỘI ĐỒNG LƯU VỰC SƠNG SỊ TỈNH NAM ĐỊNH ĐÀO ANH VĂN CHUYÊN NGÀNH: THỦY VĂN HỌC MÃ SỐ: 80440224 NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS LÃ VĂN CHÚ HÀ NỘI, NĂM 2019 CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC TÀI NGUN VÀ MƠI TRƯỜNG HÀ NỘI Cán hướng dẫn chính: PGS.TS Lã Văn Chú Cán chấm phản biện 1: TS Đoàn Minh Trí Cán chấm phản biện 2: PGS.TS Nguyễn Thanh Hùng Luận văn thạc sĩ bảo vệ tại: HỘI ĐỒNG CHẤM LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI Ngày 20 tháng 01 năm 2019 i LỜI CAM ĐOAN Tôi xin cam đoan nội dung, số liệu, kết nêu luận văn trung thực chƣa đƣợc cơng bố cơng trình khác TÁC GIẢ LUẬN VĂN (Ký ghi rõ họ tên) Đào Anh Văn ii LỜI CẢM ƠN Sau hai năm học tập nghiên cứu khoa học anh chị em học viên Lớp Cao học Thủy văn CH3A.T (khóa 2017 - 2019) Trƣờng Đại học Tài nguyên Mơi trƣờng Hà Nội, đến tơi hồn thành khố luận tốt nghiệp Với lòng chân thành nhất, xin chân thành gửi lời cảm ơn: - Các thầy cô cán làm công tác quản lý Khoa Khí tƣợng Thủy văn, Phòng Đào tạo, Trƣờng Đại học Tài nguyên Môi trƣờng Hà Nội Viện Khí tƣợng Thủy văn Biến đổi khí hậu tạo điều kiện thuận lợi để em hoàn thành chƣơng trình Cao học luận văn tốt nghiệp - Lãnh đạo Trung tâm Quan trắc khí tƣợng thủy văn, lãnh đạo Tổng cục Khí tƣợng Thủy văn tạo điều kiện cho tơi đƣợc tham gia khóa học Cảm ơn đồng nghiệp tạo điều kiện giúp đỡ tơi cơng tác để tơi có thời gian học tập hoàn thành luận văn tốt nghiệp Trƣờng Đại học Tài nguyên Môi trƣờng Hà Nội - Các cán thuộc Viện Khoa học Khí tƣợng Thủy văn Biến đổi khí hậu, Phòng Quản lý Mạng lƣới thơng tin khí tƣợng thủy văn thuộc Đài Khí tƣợng Thủy văn khu vực Đồng Bằng Bắc Bộ, Phòng Chỉnh lý Bảo quản tƣ liệu khí tƣợng thủy văn thuộc Trung tâm Thơng tin Dữ liệu khí tƣợng thủy văn giúp đỡ tơi thu thập liệu trình thực luận văn - Các anh chị em học viên Lớp Cao học Thủy văn CH3A.T, trƣờng Đại học Tài nguyên Môi trƣờng Hà Nội quan tâm, động viên, giúp đỡ vững tâm phấn đấu suốt q trình học tập hồn thành luận văn tốt nghiệp - Đặc biệt PGS.TS Lã Văn Chú, Viện Khoa học Khí tƣợng Thủy văn Biến đổi khí hậu, ngƣời Thầy tận tình hƣớng dẫn suốt q trình tìm hiểu, nghiên cứu hồn thành luận văn tốt nghiệp; kết đạt đƣợc luận văn kiến thức khoa học quý báu mà thầy giành nhiều thời gian tâm huyết để hƣớng dẫn, bảo thời gian qua Do thời gian có hạn, số liệu thực đo chƣa đƣợc đầy đủ nhƣ mong muốn khả nghiên cứu hạn chế nên luận văn khơng thể tránh đƣợc thiếu sót Rất mong nhận đƣợc ý kiến đóng góp quý báu thầy cô, đồng nghiệp ngƣời quan tâm iii THÔNG TIN LUẬN VĂN Họ tên học viên: Đào Anh Văn Lớp: CH3A.T Khóa: III Năm học: 2017 - 2019 Cán hƣớng dẫn: PGS.TS Lã Văn Chú Tên đề tài: Nghiên cứu đánh giá tác động xâm nhập mặn cho hệ thống kênh mƣơng nội đồng lƣu vực sơng Sò tỉnh Nam Định Luận văn đƣợc thực 87 trang bao gồm ba chƣơng chính: Chƣơng 1: Tổng quan nghiên cứu xâm nhập mặn giới thiệu khu vực nghiên cứu Trong chƣơng 1, đƣa số nghiên cứu xâm nhập mặn giới Việt Nam Giới thiệu vị trí, đặc điểm tự nhiên lƣu vực nghiên cứu Đánh giá trạng xâm nhập mặn khu vực Chƣơng 2: Xây dựng chƣơng trình tính tốn xâm nhập mặn hệ thống sơng Hồng - Thái Bình sơng Sò Phân tích lựa chọn mơ hình tính tốn Xây dựng sơ đồ mạng sông Thiết lập điều kiện biên, điều kiện ban đầu Hiệu chỉnh kiểm định thông số mơ hình MIKE 11 cho hệ thống sơng Chƣơng 3: Xây dựng kịch xâm nhập mặn lƣu vực sơng Sò dƣới tác động nƣớc biển dâng đề xuất lấy nƣớc tƣới phục vụ sản xuất nông nghiệp Xây dựng đồ mặn theo kịch tháng mùa cạn cho khu vực kênh mƣơng nội đồng lƣu vực sơng Sò Đánh giá khả lấy nƣớc tƣới cho sơng Sò cống sông Hồng – Ninh Cơ dƣới ảnh hƣởng nƣớc biển dâng iv MỤC LỤC LỜI CAM ĐOAN I TÁC GIẢ LUẬN VĂN I LỜI CẢM ƠN II THÔNG TIN LUẬN VĂN III MỤC LỤC IV MỞ ĐẦU 1 Ý NGHĨA KHOA HỌC VÀ THỰC TIỄN CỦA ĐỀ TÀI LUẬN VĂN MỤC TIÊU CỦA LUẬN VĂN 3 PHƢƠNG PHÁP NGHIÊN CỨU BỐ CỤC CỦA LUẬN VĂN CHƢƠNG TỔNG QUAN CÁC NGHIÊN CỨU VỀ XÂM NHẬP MẶN VÀ GIỚI THIỆU KHU VỰC NGHIÊN CỨU 1.1 TỔNG QUAN CÁC PHƢƠNG PHÁP NGHIÊN CỨU VÀ TÍNH TỐN XÂM NHẬP MẶN VÙNG CỬA SƠNG 1.1.1 Các nghiên cứu nƣớc 1.1.2 Các nghiên cứu nƣớc GIỚI THIỆU KHU VỰC NGHIÊN CỨU 15 1.2.1 Vị trí địa lý [1] 15 1.2.2 Đặc điểm địa hình [1] 16 1.2.3 Đặc điểm khí tƣợng thủy văn 18 1.2.4 Đặc điểm thủy văn sơng ngòi 19 1.3 ĐÁNH GIÁ THỰC TRẠNG XÂM NHẬP TRIỀU, MẶN Ở KHU VỰC ĐỒNG BẰNG SƠNG HỒNG – THÁI BÌNH 21 1.3.1 Đánh giá trạng xâm nhập triều mặn vùng đồng sông Hồng – Thái Bình 21 1.3.2 Đánh giá số liệu tình hình xâm nhập mặn vùng sơng Sò 26 CHƢƠNG 34 THIẾT LẬP MÔ HÌNH TÍNH TỐN XÂM NHẬP MẶN HỆ THỐNG SƠNG HỒNG - THÁI BÌNH VÀ SƠNG SỊ 34 2.1 Lựa chọn công cụ tính tốn 34 2.1.1 Cơ sở lý thuyết mô hình MIKE 11 34 2.1.2 Phương pháp giải: 37 v a) Phương pháp giải mơ hình thủy đọng lực 37 b) Phương pháp giải mơ hình khuếch tán 39 2.2 Hiệu chỉnh kiểm định thơng số mơ hình MIKE 11 cho hệ thống sơng Hồng – Thái Bình 39 2.2.1 Xây dựng sơ đồ mạng sông 40 2.2.2 Thiết lập điều kiện biên 42 2.2.3 Thiết lập điều kiện ban đầu 45 2.2.4 Hiệu chỉnh thông số cho mơ hình 46 2.2.5 Kiểm định thông số cho mô hình MIKE 11 cho mạng sơng Hồng – Thái Bình52 2.2.6 Kiểm định thơng số cho mơ hình MIKE 11 cho mạng sơng Sò 55 2.3 Nhận xét chung kết hiệu chỉnh kiểm định mơ hình cho hệ thống sơng Hồng – Thái Bình hệ thống sơng Sò 60 CHƢƠNG 63 XÂY DỰNG KỊCH BẢN XÂM NHẬP MẶN LƢU VỰC SƠNG SỊ DƢỚI TÁC ĐỘNG CỦA NƢỚC BIỂN DÂNG SLR VÀ ĐỀ XUẤT LẤY NƢỚC TƢỚI PHỤC VỤ SẢN XUẤT NÔNG NGHIỆP 63 3.1 Sơ đồ thủy lực 63 3.1.1 Sơ đồ mạng sông 63 3.1.2 Điều kiện tính tốn 63 3.2 Xây dựng đồ mặn theo kịch tháng mùa cạn cho khu vực kênh mƣơng nội đồng lƣu vực sơng Sò 64 3.2.1 ựa chọn ịch biến đ i h hậu ĐK t nh toán m nhập mặn vào nội đồng sơng Sò năm mặn xâm nhập sâu vào hệ thống sơng Sò 64 3.2.2 Tính tốn xâm nhập mặn vào nội đồng sơng Sò năm mặn xâm nhập sâu ứng với kịch nước biển dâng – RCP4.5 65 3.3 Đánh giá khả lấy nƣớc tƣới cho sơng Sò cống sông Hồng – Ninh Cơ dƣới ảnh hƣởng nƣớc biển dâng 71 3.3.1 Mục đ ch 71 3.3.2 Xây dựng kịch tính tốn 72 3.3.3 Điều kiện tính tốn 72 3.3.4 Sơ đồ mạng sông 72 3.3.5 Kết tính tốn 73 KẾT LUẬN VÀ KIẾN NGHỊ 83 DANH MỤC TÀI LIỆU THAM KHẢO 86 vi DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT Chữ viết tắt KTTV BĐKH Khí tƣợng thủy văn Biến đổi khí hậu SLR Nƣớc biển dâng (Sea Level Rise) % Phần trăm Phần nghìn Mực nƣớc đỉnh triều trung bình nhiều năm Mực nƣớc chân triều trung bình nhiều năm ‰ Hđtb Hctb RCP AD NASH Giải thích đƣờng nồng độ khí nhà kính (Representative Concentration Pathways) Mơ đun khuếch tán bình lƣu Sai số lệch đỉnh vii DANH MỤC BẢNG ảng 1.1 Tóm tắt số mơ hình tốn thường sử dụng Việt Nam 13 Bảng 1.2 Độ mặn lớn thời gian xuất [1], [8] 22 Bảng 1.3 Khoảng cách xâm nhập mặn số sơng [10] 24 Bảng Triết giảm độ mặn triền sông [10] 25 Bảng 1.5 Độ mặn thực đo lớn tháng mùa cạn 28 Bảng Số liệu đặc trưng triều trung bình cửa sơng Sò trung bình trạm Ba Lạt Phú Lễ chuỗi số liệu trung bình năm 31 Bảng ệ số nhám sông hệ thống sơng ồng – Thái ình 49 Bảng 2 Kết đánh giá sai số tính tốn thực đo số trạm kiểm tra 50 Bảng Kết đánh giá sai số độ mặn tính tốn thực đo vị trí kiểm tra hệ thống sơng Hồng – Thái ình cho tháng năm 2006 51 Bảng Kết đánh giá sai số tính toán thực đo số trạm kiểm tra hệ thống sơng Hồng – Thái Bình tháng năm 2007 54 Bảng Kết đánh giá sai số độ mặn t nh toán thực đo vị tr iểm tra hệ thống sơng ồng – Thái ình cho tháng năm 2007 55 Bảng Các vị tr biên lưu lượng gia nhập vào hệ thống kênh, cống lưu vực sơng Sò 58 Bảng Kết độ mặn tính tốn thực đo cống Nhất Đỗi 59 Bảng Mực nước biển dâng theo kịch RCP4.5 64 Bảng Diện t ch m2 huyện bị ảnh hưởng m nhập mặn theo thời ỳ Kịch RCP4.5 67 Bảng 3 Tỷ lệ diện t ch % huyện bị ảnh hưởng m nhập mặn theo thời ỳ Kịch RCP4.5 67 Bảng Độ mặn lớn vị tr cống lấy nước sông ồng 76 Bảng Độ mặn lớn vị tr cống lấy nước sông Ninh Cơ 76 Bảng T ng thời gian lấy nước trạng phương án ĐK vị tr cống lấy nước điển hình sông ồng – Ninh Cơ 78 Bảng T ng lượng nước lấy qua cống đầu mối từ tháng I-IV (103m3) 79 Bảng Tỷ lệ biến đ i t ng lượng nước ịch biến đ i h hậu so với ịch trạng % 80 viii DANH MỤC HÌNH Hình Bản đồ lưu vực sơng Sò, tỉnh Nam Định [2] ình 2.1 a Sơ đồ sai phân hữu hạn điểm ẩn Abbott; b Sơ đồ sai ph n điểm ẩn Abbott mặt phẳng x~t 37 ình 2 Sơ đồ tính thủy lực mô mặn mạng sông Hồng – Thái Bình 41 ình Sơ đồ mạng sơng Hồng – Thái Bình giao diện Mike 11 42 ình Đường trình lưu lượng thực đo trạm Sơn T y 1/2006 44 ình Đường trình lưu lượng thực đo trạm Chũ 1/2006 44 ình Đường trình lưu lượng thực đo trạm Gia Bảy 1/2006 45 ình Đường trình lưu lượng thực đo trạm Chũ 1/2006 46 ình Đường trình lưu lượng thực đo trạm ưng Thi 1/2006 47 ình Đường trình mực nước triều trạm thủy văn cửa sông thuộc hệ thống sơng Hồng – Thái ình năm 2006 47 ình 2.10 Đường q trình độ mặn cửa sơng thuộc hệ thống sơng Hồng – Thái Bình 48 Hình 11 Giao diện hiệu chỉnh thơng số mơ hình 49 ình 12 Đường trình lưu lượng thực đo trạm Sơn T y 1/2007 53 ình 13 Đường trình lưu lượng thực đo trạm Gia Bảy 1/2007 53 ình 14 Đường trình lưu lượng thực đo trạm Chũ 1/2007 54 ình 2.15 Sơ đồ thủy lực cho hệ thống sơng Sò 56 Hình 16 Mặt cắt ngang sơng lưu vực sơng Sò giao diện mơ hình MIKE 11 57 Hình Bản đồ vùng ảnh hưởng mặn lưu vực sơng Sò trạng – Kịch A1 độ gia tăng S R=0cm 68 Hình ản đồ vùng ảnh hưởng mặn lưu vực sơng Sò ứng với năm 2030 – Kịch A2 độ gia tăng S R=13cm) 69 Hình 3 ản đồ vùng ảnh hưởng mặn lưu vực sơng Sò ứng với năm 2050 – Kịch A3 độ gia tăng S R=22cm 69 Hình ản đồ vùng ảnh hưởng mặn lưu vực sông Sò ứng với năm 2070 – Kịch A4 độ gia tăng S R=33cm) 70 Hình iễn biến q trình mặn lớn dọc sơng ồng phương án trạng Kịch A1- độ gia tăng S R 0cm 73 Hình iễn biến trình mặn lớn dọc sông ồng phương án ĐK năm 2030 – Kịch A2- độ gia tăng S R 13cm 74 74 Hình iễn biến trình mặn lớn dọc sông ĐK năm 2030 – Kịch A2- độ gia tăng S R 13cm) ồng phương án Hình iễn biến q trình mặn lớn dọc sơng ĐK năm 2050– Kịch A3- độ gia tăng S R 22cm ồng phương án Hình iễn biến trình mặn lớn dọc sơng ĐK năm 2070– Kịch A4- độ gia tăng S R 33cm ồng phương án 75 Hình iễn biến trình mặn lớn dọc sông Ninh Cơ phương án trạng– Kịch A1- độ gia tăng S R 0cm Hình 10 iễn biến trình mặn lớn dọc sông Ninh Cơ phương án ĐK năm 2030– Kịch A2- độ gia tăng S R 13cm 76 Hình 11 iễn biến q trình mặn lớn dọc sơng Ninh Cơ phương án ĐK năm 2050– Kịch A3- độ gia tăng S R 22cm Hình 12 iễn biến q trình mặn lớn dọc sơng Ninh Cơ phương án ĐK năm 2070– Kịch A4- độ gia tăng SLR 33cm) Bảng Độ mặn lớn TT Cống đầu mối vị tr cống lấy nước sông ồng KB A1 Hiện trạng KB A2 2030 KB A3 2050 KB A4 2070 Cống số 4.69 4.91 5.13 5.44 Cống Hạ Miêu 7.99 8.30 8.62 9.01 Cống Hạ Miêu 8.52 8.81 9.11 9.58 Cống Cát Xuyên 9.85 10.21 10.58 11.04 Cống Liêu Đông 10.73 11.16 11.59 12.16 Cống Tài 11.50 11.89 12.30 12.86 Cống Ngô Đồng 14.09 14.52 14.96 15.54 Bảng Độ mặn lớn TT Cống đầu mối vị tr cống lấy nước sông Ninh Cơ KB A1 Hiện trạng KB A2 2030 KB A3 2050 KB A4 2070 Cống Chợ Đê 1.65 1.77 1.90 1.95 Cống Tây Khu 1.85 1.98 2.12 2.23 Cống Đông Nê 2.09 2.23 2.39 2.52 Cống Bắc Cầu 2.13 2.27 2.43 2.57 Cống Tra Thƣợng 2.15 2.30 2.46 2.59 Cống Rộc 2.18 2.40 2.54 2.68 Cống Kèo 3.28 3.62 3.80 4.11 77 Hình 13 iễn biến trình mặn lớn vị tr cống dọc sơng ồng theo ịch ĐK Hình 14 Diễn biến trình mặn lớn vị tr cống dọc sông Ninh Cơ theo ịch ĐK 78 Bảng T ng thời gian lấy nước trạng phương án ĐKH vị trí cống lấy nước điển hình sơng Hồng – Ninh Cơ (giờ) TT 10 11 12 13 14 Tên Cống Cống Chợ Đê Cống Tây Khu Cống Đông Nê Cống Bắc Cầu Cống Tra Thƣợng Cống Rộc Cống Kèo Cống số Cống Cát Xuyên Cống Hạ Miêu Cống Hạ Miêu Cống Liêu Đông Cống Tài Cống Ngô Đồng Hiện trạng BĐKH 2030 BĐKH 2050 BĐKH 2070 Sông Tháng I Tháng II Tháng III Tháng IV Tháng I Tháng II Tháng III Tháng IV Tháng I Tháng II Tháng III Tháng IV Tháng I Tháng II Tháng III Tháng IV Ninh Cơ 744 672 710 720 744 672 675 718 744 672 639 713 744 672 590 696 Ninh Cơ 744 672 710 718 744 668 682 715 744 664 645 709 741 660 619 693 Ninh Cơ 744 667 685 713 739 663 654 705 737 659 606 694 731 650 541 682 Ninh Cơ 649 636 655 711 635 620 641 703 622 603 632 689 602 568 629 666 Ninh Cơ 631 612 655 707 611 598 636 704 586 583 615 701 557 525 603 697 Ninh Cơ Ninh Cơ Hồng 585 502 497 598 531 527 644 638 637 700 659 656 577 483 475 586 514 405 603 617 393 694 653 444 538 470 464 558 487 387 578 610 277 670 646 431 513 420 451 551 451 348 545 592 251 657 637 409 Hồng 506 389 310 431 330 273 164 304 313 216 143 280 283 177 108 262 Hồng 321 198 119 287 304 160 98 280 288 149 87 254 259 132 72 232 Hồng 286 183 90 256 237 145 75 246 218 130 70 237 213 112 66 222 Hồng 230 122 51 223 213 104 46 215 204 100 42 211 176 83 37 206 Hồng 186 93 55 206 175 81 41 200 163 74 36 194 152 60 31 183 Hồng 163 72 27 201 154 66 25 194 141 60 23 189 128 50 22 178 79 Hình 15 iễn biến thời gian lấy nước cống sơng ồng Hình 16 iễn biến thời gian lấy nước cống sông Ninh Cơ Bảng T ng lư ng nước lấy qua cống đầu mối từ tháng I-IV (103m3) KB A1 Hiện trạng KB A2 2030 KB A3 2050 KB A4 2070 Sông Hồng 9796 10477 11477 12434 Cống Hạ Miêu Sông Hồng 15042 15935 16915 17993 Cống Hạ Miêu Sông Hồng 18155 19399 20752 22261 Cống Cát Xuyên Sông Hồng 10709 11411 12078 12817 TT Cống đầu mối Cống số Tên sông 80 KB A1 Hiện trạng KB A2 2030 KB A3 2050 KB A4 2070 Cống Liêu Đông Sông Hồng 8274 8690 9167 9657 Cống Tài Sông Hồng 9093 9442 9832 10235 Cống Ngô Đồng Sông Hồng 42833 45730 48967 52620 Cống Chợ Đê Sông Ninh Cơ 8618 9137 9694 10341 Cống Tây Khu Sông Ninh Cơ 6722 7232 7679 8219 10 Cống Đông Nê Sông Ninh Cơ 12928 13071 13229 13406 11 Cống Bắc Cầu Sông Ninh Cơ 9425 9857 10538 11104 12 Cống Tra Thƣợng Sông Ninh Cơ 19320 20132 20987 22144 13 Cống Rộc Sông Ninh Cơ 7535 7851 8185 8636 14 Cống Kèo Sông Ninh Cơ 9081 9620 10276 11020 TT Cống đầu mối Tên sông Bảng Tỷ lệ biến đ i t ng lư ng nước kịch biến đ i khí hậu so với kịch trạng (%) Cống số Sông Hồng KB A2 2030 6.95 Cống Hạ Miêu Sông Hồng 5.94 11.75 17.45 Cống Hạ Miêu Sông Hồng 6.85 13.39 19.79 Cống Cát Xuyên Sông Hồng 6.56 12.00 17.45 Cống Liêu Đông Sông Hồng 5.03 10.28 15.09 Cống Tài Sông Hồng 3.84 7.83 11.62 Cống Ngô Đồng Sông Hồng 6.76 13.41 19.99 Cống Chợ Đê Sông Ninh Cơ 6.02 11.78 17.77 Cống Tây Khu Sông Ninh Cơ 7.59 13.23 19.49 10 Cống Đồng Nê Sông Ninh Cơ 1.11 2.30 3.61 11 Cống Bắc Cầu Sông Ninh Cơ 4.58 11.29 15.93 12 Cống Tra Thƣợng Sông Ninh Cơ 4.20 8.28 13.46 13 Cống Rộc Sông Ninh Cơ 4.19 8.28 13.45 14 Cống Kèo Sông Ninh Cơ 5.94 12.42 18.87 TT Cống đầu mối Tên sông KB A3 2050 16.04 KB A4 2070 22.99 81 Hình 3.17 iễn biến trình t ng lượng nước lấy qua cống đầu mối từ tháng I-IV theo ịch ĐK Nhận xét: Diễn biến xâm nhập mặn cho thấy nhiều cống lấy nƣớc đầu mối sông Hồng Ninh Cơ bị nhiễm mặn - Trên sông Hồng: Phƣơng án trạng cho thấy độ mặn lớn cống Cát Xuyên xấp xỉ 10‰, cống Ngơ Đồng độ mặn lớn tƣơng ứng 14,09% Phƣơng án nƣớc biển dâng 2030, 2050, 2070 độ mặn lớn có xu tăng lên qua thời kỳ, đến năm 2070 đạt 15,54‰ cống Ngô Đồng Trên hầu hết cống lấy nƣớc dọc sông Hồng xuất độ mặn lớn > ‰ - Trên sông Ninh Cơ: Ở phƣơng án trạng độ mặn lớn nhỏ 4‰ cống nằm phần gần sông Ninh Cơ nên ảnh hƣởng xâm nhập mặn khơng nhiều, có xu tăng dần theo kịch nƣớc biển dâng thời kỳ từ năm 2030 đến 2070 nhƣng mức dƣới 4‰ Chế độ đóng - mở cống đầu mối sông Hồng Ninh Cơ cung cấp nƣớc tƣới cho hệ thống sơng Sò biến đổi dọc theo sơng, gần cửa sơng thời gian đóng cống có xu tăng dần 82 - Trên sông Hồng: Các cống từ cửa sông đến cống Cồn Giữa gần nhƣ đóng hồn tồn, Cống Ngô Đồng mở lấy nƣớc 1-2 ngày vào thời kỳ triều Phía thƣợng lƣu, số lấy nƣớc có lớn hơn, nhiên Cống số có số ngày cống lấy đƣợc nƣớc hồn tồn ngày - Trên sơng Ninh Cơ: Các cống có lấy nƣớc lớn Trên cống nhƣ Chợ Đê, Tây Khu, Đồng Nê có khả lấy nƣớc hầu hết thời gian mùa cạn mực nƣớc cho phép để phục vụ nông nghiệp nuôi trồng thủy sản Các kịch nƣớc biển dâng không làm tăng đáng kể chiều dài xâm nhập mặn sông nhƣng độ mặn vị trí tƣơng ứng có xu tăng lên đáng kể Với quy trình vận hành cống nhƣ phƣơng án trạng, lƣợng nƣớc lấy vào nội đồng sơng Sò qua cống lấy nƣớc ứng với kịch BĐKH lớn so với kịch trạng A1, dao động từ 1-7% kịch A2 (2030), từ đến 16% kịch A3 (2050) từ 3.61% đến 22.99% kịch A4 (2070) 83 KẾT LUẬN VÀ KIẾN NGHỊ 1- KẾT LUẬN Đề tài luận văn thực đƣợc nội dung sau đây: 1- Phân tích lựa chọn đƣợc chƣơng trình MIKE 11 để tính tốn q trình xâm nhập mặn vào hệ thống: Hệ thống sơng Hồng - Thái Bình hệ thống sơng Sò thuộc tỉnh Nam Định Đây chƣơng trình đƣợc áp dụng cho nhiều hệ thống vùng cửa sơng Việt Nam có đủ độ tin cậy áp dụng cho khu vực nghiên cứu 2- Đã tiến hành hiệu chỉnh kiểm định thông số mơ hình MIKE 11 theo tiêu chuẩn hành Việc mơ hình hóa q trình thủy lực mùa cạn đƣợc thực cho mùa cạn điển hình năm 2006 Kết hiệu chỉnh kiểm định cho thấy sai số vị trí kiểm tra dao động từ 2-11% Nhƣ vậy, thông số sử dụng mơ hình ổn định đƣợc sử dụng tiếp phần tính tốn thời kỳ biến đổi khí hậu 3- Lập đồ xâm nhập mặn cho lƣu vực sơng Sò dƣới tác động biến đổi khí hậu theo kịch nƣớc biển dâng RCP4.5 cho thời kỳ 2030, 2050 2070, báo cáo tiến hành thu thập, xử lý đánh giá phân tích tài liệu khí tƣợng thủy văn trạm lƣu vực sông phục vụ làm đầu vào cho mô hình thủy lực tính tốn xâm nhập mặn Qua q trình mơ dự báo thấy tƣơng lai dƣới tác động biến đổi khí hậu, tình hình xâm nhập mặn lƣu vực sơng Sò ngày trở nên nghiêm trọng Kết tính tốn cho thấy, huyện Hải Hậu có nguy xâm nhập mặn cao nhất, huyện Giao Thủy với mức độ thấp nhƣng tình trạng đáng báo động Nguy xâm nhập mặn huyện tƣơng lai gây nhiều khó khăn cản trở cho công tác lấy nƣớc phục vụ tƣới nhƣ yêu cầu dân sinh kinh tế khu vực Đặc biệt huyện ven biển vấn đề cấp nƣớc cho ngành trở nên nghiêm trọng tình trạng xâm nhập mặn ngày lấn sâu vào đất liền 84 4- Đã tính tốn xây dựng kịch lấy nƣớc theo Kịch Biến đổi khí hậu RCP 4.5 cho thấy, hầu nhƣ tồn diện tích trƣớc cống ngăn mặn từ sơng Sò bị nhiễm mặn Về kết luận lấy nƣớc tƣới cho khu vực nghiên cứu: - Lƣợng nƣớc qua cống lấy nƣớc dọc sông Hồng – Ninh Cơ tăng theo kịch BĐKH với quy trình vận hành nhƣ phƣơng án trạng, nhiên khả nguy nhiễm mặn vào hệ thống nội đồng sơng Sò cao - Số mở cống sông Ninh Cơ sông Hồng tập trung nhiều vào tháng I, IV; số mở cống tháng cống sông Ninh Cơ cao số mở cống sông Hồng, thời gian lấy nƣớc vào tháng sông Hồng cho thời kỳ 2050, 2070 Trên sông Ninh Cơ, cống nằm xa biển nên khả nhiễm mặn thấp nên thời gian lấy nƣớc tháng mùa cạn có xu không thay đổi 5- Một số tồn kết nghiên cứu: + Trong luận văn sử dụng tài liệu địa hình mặt cắt hệ thống sơng HồngThái Bình từ nhiều nguồn nhƣng với thời kỳ đo đạc khơng liên tục từ trƣớc năm 2000 có thay đổi điều kiện lòng dẫn nên có ảnh hƣởng đến kết tính tốn thủy lực + Số liệu biên (Q, H) hệ thống mạng sơng Hồng – Thái Bình số liệu thực đo nhƣng chế độ đo khơng đồng nên tính toán phải sử dụng phƣơng pháp nội suy trình nên có ảnh hƣởng đến kết hiệu chỉnh kiểm định mơ hình + Do khơng có trạm đo triều cửa sơng Sò (biên dƣới) nên sử dụng kết nội suy từ trạm mực nƣớc triều từ cửa sông Hồng Ninh Cơ nên có ảnh hƣởng đến độ tin cậy kết tính tốn 2- KIẾN NGHỊ 1) Tổ chức đo đạc, khảo sát số liệu địa hình bao gồm mặt cắt dọc, ngang hệ thống sông Hồng – Thái Bình số liệu cũ (từ năm 1999- 2000) không phản ánh sát với diễn biến địa hình 85 2) Bổ sung thêm trạm đo đạc thủy văn độ mặn tồn hệ thống sơng, chủ yếu khu vực cửa sơng để thu thập số liệu giúp công tác hiệu chỉnh cho mơ hình tính tốn dự báo xâm nhập mặn 3) Bổ sung thêm kịch Biến đổi khí hậu với mức biển dâng khác để có nhìn tổng qt khả ảnh hƣởng triều, mặn đến vùng cửa sơng Sò 4) Tăng cƣờng công tác dự báo xâm nhập triều mặn để phục vụ lấy nƣớc cho sản xuất nông nghiệp vùng cửa sông 86 DANH MỤC TÀI LIỆU THAM KHẢO Công ty TNHH thành viên khai thác CTTL Xuân Thủy, Báo cáo quy trình vận hành hệ thống thủy nơng Xn Thủy trường h p có mưa lớn gây ngập úng, 10/2014 Nguyễn Lan Châu, Đánh giá tác động hệ thống hồ chứa Tun Quang, Thác à, Hồ ình đến chế độ dòng chảy mùa cạn hạ lưu sơng Hồng đề xuất giải pháp đảm bảo mực nước hạ du Đề cƣơng đề tài nghiên cứu cấp Bộ, 2008-2010) Nguyễn Tất Đắc, Ảnh hưởng gió chướng lưu lư ng nguồn tới xâm nhập mặn đồng sơng Cửu Long, Tạp chí KTTV tháng số 463, năm 1999 Nguyễn Tất Đắc, Nguyễn Văn Điệp, Nguyễn Minh Sơn, 1988, Mơ hình tính tốn dòng chảy chất lư ng nước hệ thống kênh, sông (WFQ87) kỹ thuật chương trình, Uỷ ban Quốc gia Chƣơng trình Thủy văn Quốc tế Việt Nam Lã Thanh Hà, Đỗ Văn Tuy, 1999, Tính tốn lập phương án dự báo xâm nhập mặn mô hình SALHO cho vùng cửa sơng TP Hải Phòng, Đề tài NCKH cấp Thành phố Lã Thanh Hà, Nghiên cứu khả dự báo xâm nhập mặn vùng đồng sơng Hồng – sơng Thái Bình mơ hình tốn, Tạp chí KTTV tháng số 523 năm 2003 Lã Thanh Hà (2013), Đánh giá đặc điểm tài ngun nước mặt hệ thống sơng địa bàn tỉnh Nam Định, Đề tài NCKH cấp Tỉnh, Sở Khoa học Công nghệ, tỉnh Nam Định Lã Thanh Hà (2004), “Tính tốn thủy lực phục vụ dự án nạo vét, cải tạo hệ thống tiêu sông Sò tỉnh Nam Định”Dự án:” Phục vụ thiết kế chỉnh trị sơng Sò làm sở thiết kế trị sông, tăng cƣờng khả tiêu cho hệ thống thủy nông Xuân Thủy, tỉnh Nam Định”, Bộ NN&PTNT Lã Thanh Hà, 2014, “Xây dựng mơ hình dự báo xâm nhập mặn hạ lưu Sông Mã, Sông Hoạt, Sông Yên, tỉnh Thanh Hóa” Nhiệm vụ 01 thuộc Đề án: “Xây dựng mơ hình dự báo xâm nhập mặn hạ lưu Sơng Mã, Sơng oạt, Sơng n, tỉnh Thanh óa”, UBND tỉnh Thanh Hóa 10 Đồn Thanh Hằng & Lã Thanh Hà (2006-2009), Xây dựng chương trình dự báo xâm nhập mặn cho khu vực đồng sông Hồng – Thái Bình, Đề tài NCKH cấp Bộ, Bộ TN&MT 11 Hà Thanh Hƣơng, ứng dụng mơ hình tốn số trị tính tốn xâm nhập mặn hạ lưu sơng Văn Úc, Tạp chí KTTV tháng số 465 năm 1999 87 12 Lê Bắc Huỳnh, Thực trạng suy giảm nguồn nước hạ lưu lưu vực sông vấn đề đặt quản lý, Tạp chí KTTV tháng số 557, năm 2007 13 Trần Văn Hƣng, 1996, Nghiên cứu môi trường tài nguyên nước phục vụ phát triển kinh tế- xã hội dải ven biển Đồng Sơng Hồng Luận án Phó tiến sỹ khoa học Địa lý- Địa chất, Hà Nội 14 Nguyễn Nhƣ Khuê, 1986,Modelling of tidal propagation and salility intrusion in the Mekong main estuarine system - Technical paper, Mekong Secretariat, December 15 Trịnh Đình Lƣ, Nghiên cứu đánh giá ảnh hưởng điều tiết hồ Hòa ình đến xâm nhập mặn vùng cửa sơng Hồng sơng Thái Bình, Đề tài NCKH cấp Tổng cục, năm 2001 16 Nguyễn Hữu Nhân, 2002, Phần mềm thủy lực HydroGis-Thuyết minh kỹ thuật, Dự án án tiến kỹ thuật tiến „Xây dựng phần mềm mô tả lũ lụt xâm nhập mặn vùng đồng sông Cửu Long”, Đề tài NCKH cấp Tổng cục KTTV 17 Hoàng Niêm, 1993, Báo cáo đề tài: “Đánh giá trạng sử dụng tài nguyên nước quan điểm sinh thái phát triển lâu bền” Đề tài NCKH cấp Nhà nƣớc (KT-02-10) 18 Trần Văn Phúc, Sơ đánh giá khả ảnh hưởng hồ Hồ Bình hồ Trị An tình hình xâm nhập mặn vào hạ lưu sơng Hồng sơng Đồng Nai, Tạp chí KTTV tháng số 356 năm 1990 19 Trần Văn Phúc, 1990, Mơ hình hố q trình xâm nhập mặn đồng sông Hồng sông Cửu Long, Đề tài NCKH cấp Tổng cục 20 Trần Văn Phúc, Ảnh hưởng hồ Hồ ình đến xâm nhập mặn vùng hạ lưu sơng Hồng – sơng Thái Bình, Đề tài NCKH cấp Tổng cục, năm 1991 21 Tăng Đức Thắng, Sử dụng toán lan truyền nguồn nước thành phần: cách tiếp cận thứ hai nghiên cứu xâm nhập mặn theo nguồn nước khơng mặn, Tạp chí KTTV tháng số 537 năm 2005 22 Trần Thục (2005), Xây dựng công nghệ tính tốn dự báo lũ lớn hệ thống sơng Hồng - Thái Bình Đề tài NCKH cấp Bộ, Bộ TN & MT 23 Lê Sâm, 2004, Nghiên cứu xâm nhập mặn phục vụ phát triển kinh tế - xã hội vùng ven biển đồng sông Cửu Long, Đề tài NCKH cấp Nhà nƣớc, KC-08.18 thuộc Chƣơng trình Bảo vệ Mơi trƣờng Phòng tránh thiên tai, mã số KC-08 24 Trần Thanh Xuân, Ranh giới triều mặn đồng sơng Hồng, Tạp chí KTTV tháng 11 số 347 năm 1989 25 Hoàng Minh Tuyển,2014 “Nghiên cứu diễn biến, xác định nguyên nhân thay đ i tỷ lệ phân phối dòng chảy sơng Hồng sang sông Đuống đề xuất định 88 hướng giải pháp nhằm đảm bảo tỷ lệ phân phối dòng chảy h p lý”, Đề tài NCKH cấp Bộ, Bộ TN & MT 26 Vi Văn Vỵ, Xâm nhập mặn đồng Bắc Bộ, Viện Khí tƣợng Thủy văn, năm 1986 27 Bộ Tài nguyên Môi trƣờng 2007 Báo cáo hiệu xả nước làm vụ Đông – Xuân 2006-2007 giải pháp khắc phục hạn hạ lưu sơng Hồng Báo cáo Chính phủ Bộ Tài ngun Môi trƣờng số 91/BTNMT, ngày 31 tháng năm 2007 28 Brett F Sanders1 Michael Piasecki, Mitigation of Salinity Intrusion in Wellmixed Estuaries by Optimization of Freshwater Diversion Rates, J Hydr Engrg Volume 128, Issue 1, pp 64-77; 29 Eric Wolanski Bernard Cassagne, Salinity intrusion and rice farming in the mangrove-fringed Konkoure River delta, Guinea, Wetlands Ecology and Management, Volume 8, 28 – 36; 30 Muralikrishna, Rathna Devanathan, Circulation and salinity distribution in coastal inlets, Coastal Engineering, Volume 2, 1978, Pages 119-131;