Ứng dụng mô hình SWAT đánh giá tác động của BĐKH đến dòng chảy và bùn cát vào hồ thủy điện Nậm Mức trên sông Nậm Mức

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN PHÙNG THỊ THU TRANG ỨNG DỤNG MÔ HÌNH SWAT ĐÁNH GIÁ TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN DÒNG CHẢY VÀ BÙN CÁT VÀO HỒ THỦY ĐIỆN NẬM MỨC TRÊN SÔNG NẬM MỨC LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội, 2017 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN PHÙNG THỊ THU TRANG ỨNG DỤNG MƠ HÌNH SWAT ĐÁNH GIÁ TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN DỊNG CHẢY VÀ BÙN CÁT VÀO HỒ THỦY ĐIỆN NẬM MỨC TRÊN SÔNG NẬM MỨC Chuyên ngành: Thủy văn học Mã số: 60440224 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Huỳnh Thị Lan Hương Hà Nội, 2017 LỜI CAM ĐOAN Tôi xin cam đoan Những nội dung luận văn thực hướng dẫn PGS.TS Huỳnh Thị Lan Hương Mọi tham khảo dùng luận văn trích dẫn nguồn gốc r ràng C c nội dung nghiên cứu kết luận văn trung thực chưa công bố cơng trình Xin trân trọng cảm ơn! Tác giả Phùng Thị Thu Trang i LỜI CẢM ƠN Luận văn thạc sĩ: “Ứng dụng Mơ hình SWAT đ nh gi t c động biến đổi khí hậu đến dịng chảy bùn cát vào hồ thủy điện Nậm Mức sơng Nậm Mức” hồn thành Trong q trình học tập, nghiên cứu hồn thành luận văn, học viên nhận nhiều giúp đỡ thầy cơ, bạn bè gia đình Trước hết học viên xin gửi lời cảm ơn chân thành đến cô giáo PGS.TS Huỳnh Thị Lan Hương trực tiếp hướng dẫn giúp đỡ trình nghiên cứu hoàn thành luận văn Học viên chân thành cảm ơn tới anh chị đồng nghiệp, bạn bè Trung tâm Nghiên cứu Biến đổi khí hậu Phòng Khoa học, Đào tạo Hợp tác quốc tế, Viện Khoa học Khí tượng Thủy văn Biến đổi khí hậu hỗ trợ chun mơn, thu thập tài liệu liên quan để luận văn hoàn thành Xin gửi lời cảm ơn c c thầy cô giáo Khoa Khí tượng – Thủy văn – Hải dương học – Trường Đại học Khoa học Tự nhiên, Đại học Quốc Gia Hà Nội giảng dạy, truyền đạt kiến thức, tạo điều kiện hướng dẫn học viên hoàn thành chương trình học tập thực luận văn Trong khuôn khổ luận văn, thời gian điều kiện hạn chế nên khơng tránh khỏi thiếu sót Vì vậy, học viên mong nhận ý kiến đóng góp q báu thầy giáo đồng nghiệp Xin trân trọng cảm ơn! Tác giả Phùng Thị Thu Trang 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 MỞ ĐẦU CHƢƠNG I: TỔNG QUAN 1.1 Tổng quan khu vực nghiên cứu 1.1.1 Đặc điểm địa lý tự nhiên 1.1.2 Đặc điểm kinh tế xã hội 1.2 Tổng quan nghiên cứu bùn cát nƣớc 13 1.2.1 Tổng quan nghiên cứu bùn cát giới 13 1.2.2 Tổng quan nghiên cứu bùn cát Việt Nam 14 1.3 Tổng quan biến đổi khí hậu kịch biến đổi khí hậu 17 1.3.1 Kịch biến đổi khí hậu Việt Nam 17 1.3.2 Kịch biến đổi khí hậu khu vực nghiên cứu 20 1.4 Khung nghiên cứu 25 CHƢƠNG II: PHƢƠNG PHÁP NGHIÊN CỨU VÀ NGUỒN SỐ LIỆU 26 2.1 Giới thiệu mơ hình nghiên cứu (mơ hình SWAT) 26 2.2 Tình hình số liệu, liệu khu vực nghiên cứu 44 2.2.1 Tình hình số liệu lưu lượng bùn cát khu vực nghiên cứu 44 2.2.2 Tình hình số liệu, liệu cho thời kì 47 2.2.3 Tình hình số liệu, liệu cho kịch 47 2.3 Thiết lập tính tốn dịng chảy bùn cát khu vực nghiên cứu 48 2.3.1 Thiết lập mơ hình 48 iii 2.3.2 Hiệu chỉnh kiểm địn ộ thơng số mơ hình 53 CHƢƠNG 3: ĐÁNH GIÁ TÁC ĐỘNG CỦA BĐKH ĐẾN DÒNG CHẢY VÀ BÙN CÁT vào hồ thủy điện nậm mức SÔNG NẬM MỨC 64 3.1 Tính tốn dịng chảy bùn cát thời kỳ 1986 - 2005 64 3.1.1 Kết tính tốn dịng chảy thời kỳ 64 3.1.2 Kết tính tốn bùn cát thời kỳ 64 3.2 Tính tốn dịng chảy bùn cát theo kịch biến đổi khí hậu 64 3.2.1 Xây dựng kịch biến đổi khí hậu lượng mưa n iệt độ 64 3.2.2 Kết tính tốn dịng chảy bùn cát theo kịch biến đổi khí hậu 72 3.3 Đánh giá tác động biến đổi khí hậu đến dịng chảy bùn cát 80 3.3.1 K ung đán giá tác động biến đổi khí hậu 80 3.3.2 Tác động biến đổi khí hậu đến dịng chảy 81 3.3.3 Tác động biến đổi khí hậu đến bùn cát 84 KẾT LUẬN VÀ KIẾN NGHỊ 87 TÀI LIỆU THAM KHẢO 89 iv DANH MỤC BẢNG Bảng 1.1: Lưu lượng nước trung bình tháng trung bình nhiều năm cửa Bảng 1.2: Các trị số đặc trưng năm 2016 Bảng 1.3: Biến đổi nhiệt độ trung bình năm (oC) so với thời kỳ 1986 – 2005 21 Bảng 1.4: Biến đổi lượng mưa năm (%) so với thời kỳ 1986-2005 24 Bảng 2.1: Diện tích % loại đất khu vực nghiên cứu 52 Bảng 2.2: Các loại hình sử dụng đất lưu vực nghiên cứu 52 Bảng 2.3: Tiêu chuẩn phân loại mức độ xác kết mô theo số NSI PBIAS 54 Bảng 2.4: Nhóm thơng số sau hiệu chỉnh 58 Bảng 2.5: Đ nh gi kết mô dòng chảy giai đoạn hiệu chỉnh kiểm định khu vực nghiên cứu 60 Bảng 2.6: Nhóm thơng số sau hiệu chỉnh dịng chảy bùn cát 61 Bảng 2.7: Đ nh gi kết mô bùn c t giai đoạn hiệu chỉnh kiểm định khu vực nghiên cứu 63 Bảng 3.1:C c mơ hình sử dụng tính tốn cập nhật KBBĐKH 66 Bảng 3.2: Sự thay đổi lượng mưa (%) thời kì theo kịch so với thời kì 70 Bảng 3.3: Sự thay đổi nhiệt độ trung bình (oC) thời kì theo kịch so với thời kì 70 Bảng 3.4: Sự thay đổi nhiệt độ tối cao trung bình (oC) thời kì theo kịch so với thời kì 71 Bảng 3.5: Sự thay đổi nhiệt độ tối thấp trung bình (oC) thời kì theo kịch so với thời kì 71 Bảng 3.6: Thay đổi dịng chảy trung bình tháng thời đoạn ứng với kịch so với thời kì (m3/s) 77 Bảng 3.7: Tỷ lệ thay đổi dịng chảy trung bình tháng thời đoạn ứng với kịch so với thời kì (%) 77 v Bảng 3.8: Thay đổi tổng lượng bùn cát trung bình tháng thời đoạn ứng với kịch so với thời kì (m3/s) 81 Bảng 3.9: Tỷ lệ thay đổi tổng lượng bùn cát trung bình tháng thời đoạn ứng với kịch so với thời kì (%) 81 Bảng 3.10: Dòng chảy Tỷ lệ thay đổi dịng chảy trung bình năm thời đoạn ứng với kịch so với thời kì (%) 81 Bảng 3.11: Dòng chảy Tỷ lệ thay đổi dòng chảy mùa lũ thời đoạn ứng với kịch so với thời kì (%) 82 Bảng 3.12: Dòng chảy Tỷ lệ thay đổi dòng chảy trung bình mùa kiệt thời đoạn ứng với kịch so với thời kì (%) 83 Bảng 3.13: Dòng chảy Tỷ lệ thay đổi bùn c t trung bình năm thời đoạn ứng với kịch so với thời kì (%) 84 Bảng 3.14: Dòng chảy Tỷ lệ thay đổi bùn c t mùa lũ thời đoạn ứng với kịch so với thời kì (%) 85 Bảng 3.15: Bùn cát Tỷ lệ thay đổi dịng chảy trung bình mùa kiệt thời đoạn ứng với kịch so với thời kì (%) 86 vi DANH MỤC HÌNH Hình 1.1: Bản đồ tỉnh Điện Biên Hình 1.2: Biến đổi nhiệt độ trung bình năm (oC) theo kịch RCP4.5 21 Hình 1.3: Biến đổi nhiệt độ trung bình năm (oC) theo kịch RCP8.5 21 Hình 1.4: Biến đổi nhiệt độ tối cao trung bình năm (oC) theo kịch RCP4.5 22 Hình 1.5: Biến đổi nhiệt độ tối cao trung bình năm (oC) theo kịch RCP8.5 22 Hình 1.6: Biến đổi nhiệt độ tối thấp trung bình năm (oC) theo kịch RCP4.5 23 Hình 1.7: Biến đổi nhiệt độ tối thấp trung bình năm (oC) theo kịch RCP8.5 23 Hình 1.8: Biến đổi lượng mưa năm theo kịch RCP4.5 24 Hình 1.9: Biến đổi lượng mưa năm theo kịch RCP8.5 25 Hình 1.10: Sơ đồ khối tính dịng chảy bùn cát gia nhập hồ chứa 25 Hình 2.1: Sơ đồ phát triển mơ hình SWAT 28 Hình 2.2: Sự khác phân phối độ ẩm theo chiều sâu mô theo Green Ampt thực tế 32 Hình 2.3: Diễn biến bùn cát khu vực nghiên cứu 44 Hình 2.4: Tiến trình mơ SWAT 48 Hình 2.5: Bản đồ DEM khu vực nghiên cứu 49 Hình 2.6: C c trạm khí tượng 50 Hình 2.7: C c trạm thủy văn 50 Hình 2.8: Bản đồ phân chia lưu vực thành tiểu lưu vực 50 Hình 2.9: Bản đồ thổ nhưỡng 51 Hình 2.10: Bản đồ sử dụng đất 51 Hình 2.11: Sự thay đổi lưu lượng với số CH_N1 khác 55 Hình 2.12: Sự thay đổi lưu lượng với số ALPHA_BF khác 56 Hình 2.13: Sự thay đổi lưu lượng với số OV_N khác 56 Hình 2.14: Sự thay đổi lưu lượng với số CN2 khác 57 Hình 2.15: Sự thay đổi lưu lượng với số GW_DELAY khác 57 Hình 2.16: Lưu lượng thực đo mô giai đoạn hiệu chỉnh 58 Hình 2.17: Tương quan lưu lượng thực đo mô 59 Hình 2.18: Lưu lượng thực đo mơ giai đoạn kiểm định 59 Hình 2.19:Tương quan lưu lượng thực đo mơ 60 vii Hình 2.20: Tổng lượng bùn cát thực đo mô giai đoạn hiệu chỉnh 61 Hình 2.21: Tương quan Tổng lượng bùn cát thực đo mô 61 Hình 2.22:Tổng lượng bùn cát thực đo mô giai đoạn kiểm định 62 Hình 2.23:Tương quan tổng lượng bùn cát thực đo mơ 62 Hình 3.1: Sơ đồ chi tiết hóa động lực 65 Hình 3.2: Kịch RCP 4.5 72 Hình 3.3: Kịch RCP 8.5 72 Hình 3.4: Kết tính tốn dòng chảy tháng-Kịch RCP4.5 75 Hình 3.5: Kết tính to n thay đổi dịng chảy-Kịch RCP4.5 75 Hình 3.6: Kết tính tốn dịng chảy tháng-Kịch RCP8.5 76 Hình 3.7: Kết tính to n thay đổi dòng chảy-Kịch RCP8.5 76 Hình 3.9: Kết tính tốn tổng lượng bùn cát tháng - Kịch RCP4.5 79 Hình 3.10: Kết tính to n thay đổi tổng lượng bùn cát - Kịch RCP4.5 79 Hình 3.11: Kết tính tốn tổng lượng bùn cát tháng - Kịch RCP8.5 80 Hình 3.12: Kết tính to n thay đổi tổng lượng bùn cát - Kịch RCP8.5 80 Hình 3.12: Sơ đồ khối đ nh gi t c động BĐKH lên dòng chảy, bùn cát 80 Hình 3.13: Dịng chảy trung bình năm ứng với kịch biến đổi khí hậu 81 Hình 3.14: Dịng chảy trung bình mùa lũ ứng với kịch biến đổi khí hậu 82 Hình 3.15: Dịng chảy trung bình mùa kiệt ứng với kịch biến đổi khí hậu 83 Hình 3.16: Tổng lượng bùn c t trung bình năm ứng với kịch biến đổi khí hậu 84 Hình 3.17: Tổng lượng bùn c t trung bình mùa lũ ứng với kịch biến đổi khí hậu 85 Hình 3.18: Tổng lượng bùn cát trung bình mùa kiệt ứng với kịch biến đổi khí hậu 86 viii 2) Kết tính tốn bùn cát Các kết mô bùn cát thời kỳ theo kịch so sánh với trạng bùn cát thời kỳ (Hình 3.9 Hình 3.12) cho thấy: Theo kịch RCP 4.5 Sự biến động tăng bùn cát rõ xu thời kỳ khác kịch RCP 4.5 Nhưng nhìn chung, lượng bùn c t có xu tăng từ th ng VIII đến tháng XI theo xu tăng dịng chảy Trong đó, tháng X có xu hướng tăng nhiều so với tháng lại: giai đoạn đầu kỷ 85%, giai đoạn kỷ 193%, giai đoạn cuối kỷ 337% Th ng th ng có lượng dịng chảy tăng mạnh năm Tuy nhiên, so với biến động dòng chảy, bùn cát có biến động mạnh c c giai đoạn khác Trong biến động tăng dòng chảy thời kỳ khoảng từ 20-30%, bùn c t dao động mạnh khoảng từ 100-150% thời kỳ đầu, cuối kỷ Sự biến động giảm bùn c t thể tháng I-VII tháng XII Trong đó, tổng lượng bùn cát th ng II, III có xu hướng giảm nhiều so với tháng lại, điều phù hợp với xu giảm dòng chảy tháng Dòng chảy giảm vào VII, 57% (2016-2035), 16% (20462065), chí có xu hướng tăng nhẹ 4% (2080-2099) Dựa theo đặc điểm thủy văn lưu vực th ng thường có xuất lũ kép, lũ lớn (ví dụ trận lũ lịch sử xuất vào ngày 17/VII/1994) Theo kịch RCP 8.5 Cũng tương tự kịch RCP 4.5 tổng lượng bùn cát trung bình tháng tháng VIII-XI có xu hướng tăng so với thời kỳ Trong đó, th ng X th ng có lượng bùn cát tăng lớn nhất: 238% (2016-2035), 291% (2046-2065), 455% (2080-2099) Sự biến động tăng bùn cát dao động mạnh thời kỳ từ 50-170% Tổng lượng bùn cát trung bình tháng tháng I-VII tháng XII có xu hướng giảm so với thời kỳ Trong đó, th ng có lượng dịng chảy giảm nhiều tháng II: 100% (2016-2035), 100% (2046-2065), 99% (2080-2099) Tháng 78 VII có lượng bùn cát giảm nhất: 50% (2016-2035), 5% (2046-2065), 22% (20802099) Với kết tính tốn tổng lượng bùn cát trung bình tháng theo kịch RCP 4.5 RCP 8.5 ta thấy xu hướng biến đổi bùn cát thời kì đầu kỉ (2016-2035), kỷ (2036-2065) cuối kỳ (2080-2095) kh tương đồng với tương đồng với biến đổi dịng chảy Bùn cát có xu hướng tăng vào tháng từ tháng VIII đến tháng XI giảm vào tháng I-IV tháng XII: tăng nhiều vào tháng X, tăng vào tháng VIII, chí giảm 7% 2% vào giai đoạn đầu kỷ tương ứng với kịch RCP 4.5 8.5; giảm nhiều vào tháng II vào tháng VII Thời kì cuối kỉ, biến động tăng bùn cát lớn so với thời kỳ đầu kỷ vào c c th ng mùa lũ, biến động giảm bùn cát lại rõ nét kịch Hình 3.8: Kết tính toán tổng ƣợng bùn cát tháng - Kịch RCP4.5 Hình 3.9: Kết tính tốn thay đổi tổng ƣợng bùn cát - Kịch RCP4.5 79 Hình 3.10: Kết tính tốn tổng ƣợng bùn cát tháng - Kịch RCP8.5 Hình 3.11: Kết tính tốn thay đổi tổng ƣợng bùn cát - Kịch RCP8.5 80 Đánh giá tác động biến đổi khí hậu đến dịng chảy bùn cát 3.3 3.3.1 K ung đán giá tác động biến đổi khí hậu Sơng ngịi sản phẩm khí hậu Mọi thay đổi khí hậu tác động mạnh mẽ lên nguồn nước sông Dựa Thông báo quốc gia lần thứ Việt Nam cho Khung công ước Liên hợp quốc biến đổi khí hậu, hai tài liệu “Hướng dẫn đ nh gi t c động Biến đổi khí hậu” “T c động biến đổi khí hậu đến tài nguyên nước Việt Nam” Viện Khoa học Khí tượng Thủy văn Biến đổi khí hậu việc đ nh gi t c động biến đổi khí hậu đến tài ngun nước nói chung dòng chảy, bùn cát đến hồ thủy điện Nậm Mức sơng Nậm Mức nói riêng mơ tả Hình 3.12 Hình 3.12: Sơ đồ khối đánh giá tác động BĐKH lên dòng chảy, bùn cát Các nội dung đ nh gi t c động biến đổi khí hậu lên dịng chảy bùn cát khu vực nghiên cứu chia thành khối lớn mối liên kết chặt chẽ thành phần Trong đó, c c đặc trưng dịng chảy đ nh gi bao gồm: dịng chảy trung bình năm, dịng chảy trung bình mùa kiệt, dịng chảy trung bình mùa lũ Tương tự vậy, bùn c t đ nh gi thông qua: tổng lượng bùn c t năm, tổng lượng bùn cát mùa lũ, tổng lượng bùn cát mùa kiệt Từ kết tính tốn dịng chảy bùn c t tương ứng với kịch biến đổi khí hậu thời kỳ tương lai (2016-2035, 2046-2065, 2080-2099) trạm thủy văn Nậm Mức sơng Nậm Mức rút kết luận t c động biến 80 đổi khí hậu theo hai kịch RCP4.5 RCP8.5 lên dòng chảy bùn cát khu vực nghiên cứu 3.3.2 Tác động biến đổi khí hậu đến dịng chảy 1) Tác động biến đổi khí hậu đến dịng chảy trung bình năm Dịng chảy trung bình năm có xu tăng giai đoạn cuối kỷ khơng có chênh lệch nhiều hai kịch RCP4.5 RCP8.5 Điều giải thích lượng mưa năm tăng nhiệt độ tăng nên lượng nước tổn thất bốc tho t nước lưu vực tăng lên, dẫn đến lượng dịng chảy tăng khơng đ ng kể Dịng chảy trung bình năm có xu giảm giai đoạn đầu kỷ Nhìn chung, t c động biến đổi khí hậu đến dịng chảy năm có xu giống hai kịch RCP4.5 RCP8.5, mức độ biến đổi dòng chảy trung bình năm khơng đ ng kể so với thời kỳ nền, dao động khoảng từ 226% Bảng 3.10: Dòng chảy Tỷ lệ thay đổi dòng chảy trung bình năm ứng với kịch biến đổi khí hậu Thời đoạn 1986-2005 2016-2035 2046-2065 2080-2099 Dịng chảy trung bình năm (m3/s) RCP4.5 RCP8.5 95 95 71 82 92 101 97 107 Tỷ ệ thay đổi (%) RCP4.5 RCP8.5 -26% -3% 2% -13% 6% 13% Hình 3.13: Dịng chảy trung bình năm ứng với kịch biến đổi khí hậu 81 2) Tác động biến đổi khí hậu đến dịng chảy trung bình mùa ũ Tương tự dịng chảy trung bình năm, dịng chảy trung bình mùa lũ khu vực nghiên cứu có xu tăng giai đoạn cuối kỷ khơng có chênh lệch nhiều hai kịch RCP4.5 RCP8.5 Với kịch RCP4.5, dòng chảy trung bình mùa lũ tăng khoảng 9-16% Với kịch RCP8.5, dịng chảy trung bình mùa lũ tăng khoảng 3-26% Trong giai đoạn đầu kỷ dịng chảy có xu giảm Bảng 3.11: Dịng chảy Tỷ lệ thay đổi dịng chảy trung bình mùa ũ ứng với kịch biến đổi khí hậu Thời đoạn Dịng chảy trung bình mùa ũ (m3/s) RCP4.5 RCP8.5 1986-2005 151,58 151,58 2016-2035 127,97 2046-2065 2080-2099 Tỷ ệ thay đổi (%) RCP4.5 RCP8.5 147,45 -16% -3% 165,24 181,13 9% 19% 175,36 190,60 16% 26% Hình 3.14: Dịng chảy trung bình mùa ũ ứng với kịch biến đổi khí hậu 82 3) Tác động biến đổi khí hậu đến dịng chảy trung bình mùa kiệt Dưới t c động biến đổi khí hậu, dịng chảy trung bình mùa kiệt khu vực nghiên cứu có xu giảm so với kịch Tuy nhiên, xu giảm dịng chảy khơng có chênh lệch nhiều hai kịch RCP4.5 RCP8.5 Với kịch RCP4.5, dịng chảy trung bình mùa kiệt giảm khoảng 51-66% Với kịch RCP8.5, dịng chảy trung bình mùa kiệt giảm khoảng 37-55% Bảng 3.12: Dòng chảy Tỷ lệ thay đổi dịng chảy trung bình mùa kiệt ứng với kịch biến đổi khí hậu Thời đoạn Dịng chảy trung bình mùa kiệt (m3/s) RCP4.5 RCP8.5 1986-2005 38,01 38,01 2016-2035 13,04 2046-2065 2080-2099 Tỷ ệ thay đổi (%) RCP4.5 RCP8.5 16,92 -66% -55% 18,59 20,54 -51% -46% 18,22 23,79 -52% -37% Hình 3.15: Dịng chảy trung bình mùa kiệt ứng với kịch biến đổi khí hậu 83 3.3.3 Tác động biến đổi khí hậu đến bùn cát 1) Tác động biến đổi khí hậu đến tổng ƣợng bùn cát trung bình năm T c động biến đổi khí hậu đến tổng lượng bùn cát trung bình năm có xu giống hai kịch Tổng lượng bùn cát trung bình năm có xu tăng giai đoạn cuối kỷ, tương đồng với biến động dịng chảy năm, khơng có chênh lệch nhiều hai kịch RCP4.5 RCP8.5 Trong đó, tổng lượng bùn cát trung bình năm tăng nhiều vào giai đoạn cuối kỷ 55% kịch RCP 8.5 42% kịch RCP 4.5 Bảng 3.13: Tổng ƣợng Tỷ lệ thay đổi tổng ƣợng bùn cát trung bình năm ứng với kịch biến đổi khí hậu Thời đoạn 1986-2005 2016-2035 2046-2065 2080-2099 Tổng ƣợng bùn cát trung bình Tỷ ệ thay đổi (%) năm (tấn tháng) RCP4.5 RCP8.5 RCP4.5 RCP8.5 2.491.157 2.491.157 1.843.912 3.257.737 3.530.073 2.218.845 3.311.790 3.854.169 -26% 31% 42% -11% 33% 55% Hình 3.16: Tổng ƣợng bùn cát trung bình năm ứng với kịch biến đổi khí hậu 84 2) Tác động BĐKH đến tổng ƣợng bùn cát trung bình mùa ũ Dưới t c động biến đổi khí hậu, tổng lượng bùn cát trung bình mùa lũ khu vực nghiên cứu có xu tăng giai đoạn cuối kỷ, giảm giai đoạn đầu kỷ khơng có chênh lệch nhiều hai kịch RCP 4.5 RCP 8.5 Với kịch RCP 4.5, tổng lượng bùn cát trung bình mùa lũ dao động khoảng -24-45%, tăng lớn 45% Với kịch RCP 8.5, tổng lượng bùn c t trung bình mùa lũ dao động khoảng -10-57% Như thấy rõ tương đồng xu biến đổi bùn cát dịng chảy sơng mùa lũ khu vực nghiên cứu Bảng 3.14: Tổng ƣợng Tỷ lệ thay đổi tổng ƣợng bùn cát trung bình mùa ũ ứng với kịch biến đổi khí hậu Thời đoạn Tổng ƣợng bùn cát trung bình mùa ũ (tấn tháng) RCP4.5 RCP8.5 Tỷ ệ thay đổi (%) RCP4.5 RCP8.5 1986-2005 4.785.904 4.785.904 2016-2035 3.618.046 4.331.134 -24% -10% 2046-2065 6.382.978 6.482.285 33% 35% 2080-2099 6.933.591 7.503.751 45% 57% Hình 3.17: Tổng ƣợng bùn cát trung bình mùa ũ ứng với kịch biến đổi khí hậu 85 3) Tác động BĐKH đến tổng ƣợng bùn cát trung bình mùa kiệt Tương tự dịng chảy, t c động biến đổi khí hậu, tổng lượng bùn cát trung bình mùa kiệt khu vực nghiên cứu có xu giảm mùa kiệt Đối với mùa kiệt có khác biệt rõ nét hai kịch RCP 4.5 RCP 8.5 Với kịch RCP 4.5, giảm giai đoạn từ 33-64%, tổng lượng bùn cát trung bình mùa kiệt giảm mạnh vào giai đoạn đầu kỷ Với kịch RCP 8.5, tổng lượng bùn cát trung bình mùa kiệt giảm mạnh vào giai đoạn đầu kỷ 46%, kỷ 28%, lại có xu tăng nhẹ vào cuối kỷ 4% Bảng 3.15: Tổng ƣợng Tỷ lệ thay đổi tổng ƣợng bùn cát trung bình mùa kiệt ứng với kịch biến đổi khí hậu Thời đoạn 1986-2005 2016-2035 2046-2065 2080-2099 Tổng ƣợng bùn cát trung bình mùa kiệt (tấn tháng) RCP4.5 RCP8.5 196.409 69.779 132.496 126.555 196.409 106.556 141.294 204.588 Tỷ ệ thay đổi (%) RCP4.5 RCP8.5 -64% -33% -36% -46% -28% 4% Hình 3.18: Tổng ƣợng bùn cát trung bình mùa kiệt ứng với kịch biến đổi khí hậu 86 KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Luận văn sử dụng thành cơng mơ hình SWAT để mơ diễn biến lưu lượng dòng chảy bùn cát vào hồ thủy điện Nậm Mức lưu vực sông Nậm Mức tính thay đổi dịng chảy bùn cát theo kịch biến đổi khí hậu Điều có ý nghĩa việc tính tốn mơ dòng chảy bùn c t trường hợp lưu vực khơng có số liệu đo đạc đầy đủ dòng chảy, bùn cát Tại c c lưu vực sử dụng mơ hình SWAT để bổ sung số liệu cần thiết Tuy nhiên, để nhận kết với độ x c cao cần khảo s t, đo đạc bùn c t, địa hình, thành phần hạt c c biên mơ hình số ví trí dọc bờ để kiểm định mơ hình Đối với dịng chảy: Sau phân tích độ nhạy thơng số mơ hình b ng phần mềm SWAT – CUP ta tìm thơng số có ảnh hưởng tới lưu lượng dịng chảy CN2, ALPHA_ BE, GW_ DEL Y, CH_N1, OV_ N Trong thơng số CN2 nhạy ảnh hưởng mạnh tới dòng chảy Xu biến đổi dòng chảy theo kịch RCP 4.5 RCP 8.5 kh tương đồng nhau, dịng chảy có xu hướng tăng vào c c th ng VIII đến tháng XI giảm vào c c th ng I đến th ng VII th ng XII Điều thấy r ng biến động tăng giảm dịng chảy thì phân bố dịng chảy theo mùa năm có chuyển dịch so với thời kỳ Đối với bùn cát: Tìm thơng số có ảnh hưởng tới lưu lượng dịng chảy : CN2, SOL_K, SOL_ROCK, USLE_P, SPCON, SPEXP, SLSUBBSN, HRU_SLP, SOL_AWC Có tương đồng lớn xu biến đổi dòng chảy bùn cát khu vực nghiên cứu Dòng chảy bùn c t trung bình năm có xu tăng Tuy nhiên xu không đồng c c mùa năm Mùa lũ dòng chảy bùn cát có xu tăng cịn mùa kiệt ngược lại Sự biến đổi dịng chảy bùn 87 cát hai kịch biến đổi khí hậu RCP 4.5 RCP 8.5 khơng có khác biệt rõ rệt nhiều Đ nh gi t c động biến đổi khí hậu đến dịng chảy bùn cát khu vực nghiên cứu nhận thấy có tương đồng lớn biến động c c đặc trưng dòng chảy bùn c t: năm, mùa lũ, mùa kiệt hai kịch biến đổi khí hậu RCP 4.5 RCP 8.5 Kiến nghị Trong luận văn xét tới thay đổi yếu tố nhiệt độ lượng mưa mà chưa đề cập đến thay đổi các yếu tố khí tượng thủy văn kh c thay đổi thảm phủ thực vật, số nắng, tốc độ gió Ngồi ra, lưu vực Nậm Mực khơng có số liệu nhiệt độ nên luận văn phải sử dụng số liệu nhiệt độ trạm lưu vực gần kề Đây hạn chế lớn luận văn q trình tính tốn đ nh gi t c động biến đổi khí hậu đến dịng chảy bùn c t lưu vực Do đó, để mơ hình có kết mơ tốt cần phải tăng độ xác liệu đầu vào, đồ trạng sử dụng đất thổ nhưỡng phải cập nhật thường xuyên thay đổi; bổ sung thêm trạm quan trắc mưa nhiệt độ, có số liệu địa hình chi tiết, chuẩn hóa cập nhật đồ thổ nhưỡng, đồ thảm phủ rừng, đồ sử dụng đất lưu vực 88 TÀI LIỆU THAM KHẢO [1] Nguyễn Duy Liêm, Bài giảng Mơ hình SWAT [2] Bộ giáo dục đào tạo trường đại học Nông Lâm- thành phố HCM, Bài giàng đ nh gi đất đai [3] Bộ tài nguyên Môi trường (2009), Kịch biến đổi khí hậu, nước biển dâng cho Việt Nam Hà Nội [4] Bộ tài nguyên Môi trường (2012), Kịch biến đổi khí hậu, nước biển dâng cho Việt Nam Hà Nội [5] Bộ tài nguyên mơi trường (2016), Kịch biến đổi khí hậu nước biển dâng 2016 [6] Bộ Tài nguyên Môi trường (2010), Báo cáo tổng hợp “Ảnh hưởng biến đổi khí hậu đến tài nguyên nước Việt Nam” [7] Viện Khoa học Khí tượng Thủy văn Môi trường (2010), Đ nh gi t c động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng Lưu vực sơng Hồng - Thái Bình [8] Viện Khoa học Khí tượng Thủy văn Môi trường (2010), Đ nh gi t c động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng Lưu vực sơng Đồng Nai” [9] Viện Khoa học Khí tượng Thủy văn Mơi trường (2010), Đ nh gi t c động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng Lưu vực sơng Cả [10] Viện Khoa học Khí tượng Thủy văn Mơi trường (2010), Đ nh gi t c động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng Lưu vực sơng Thu Bồn [11] Viện Khoa học Khí tượng Thủy văn Môi trường (2010), Đ nh gi t c động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng Lưu vực sơng Ba 89 [12] Viện Khoa học Khí tượng Thủy văn Môi trường (2010), Đ nh gi t c động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng Đồng b ng sơng Cửu Long [13] Viện Khí tượng Thủy văn (1985), Đặc trưng hình th i lưu vực sơng Việt Nam Hà Nội [14] Trần Thanh Xuân.(2008), Đặc điểm thủy văn tài nguyên nước Việt Nam [15] Viện Khoa học Khí tượng Thủy văn Môi trường (2012), T c động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng 90