NGHIÊN CỨU XÂY DỰNG MÔ HÌNH THỦY VĂN THÔNG SỐ TẬP TRUNG TRONG DỰ BÁO LŨ CHO CÁC LƯU VỰC SÔNG Ở VIỆT NAM

VIỆN HÀN LÂM KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM VIỆN CƠ HỌC - BÁO CÁO KẾT QUẢ THỰC HIỆN ĐỀ TÀI CƠ SỞ CẤP VIỆN CƠ HỌC NĂM 2020 ĐỀ TÀI NGHIÊN CỨU XÂY DỰNG MƠ HÌNH THỦY VĂN THƠNG SỐ TẬP TRUNG TRONG DỰ BÁO LŨ CHO CÁC LƯU VỰC SÔNG Ở VIỆT NAM Chủ nhiệm Đề tài: TS Nguyễn Chính Kiên HÀ NỘI – 2020 DANH SÁCH CÁN BỘ THAM GIA THỰC HIỆN ĐỀ TÀI TT Họ tên TS Nguyễn Chính Kiên ThS Dương Thị Thanh Hương Phịng chun mơn Thủy động lực Giảm nhẹ thiên tai lưu vực Thủy động lực Giảm nhẹ thiên tai lưu vực ThS Nguyễn Tuấn Anh Thủy động lực Giảm nhẹ thiên tai lưu vực PGS TS Trần Thu Hà Thủy động lực Giảm nhẹ thiên tai lưu vực ThS Nguyễn Hồng Phong Thủy động lực Giảm nhẹ thiên tai lưu vực KS Nguyễn Thị Hằng Thủy động lực Giảm nhẹ thiên tai lưu vực CN Trần Thị Thanh Huyền Thủy động lực Giảm nhẹ thiên tai lưu vực i MỤC LỤC DANH MỤC KÝ HIỆU, CHỮ VIẾT TẮT THƯỜNG SỬ DỤNG iii DANH MỤC CÁC BẢNG iv DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ v MỞ ĐẦU CHƯƠNG I: TỔNG QUAN VỀ MƠ HÌNH THỦY VĂN 1.1 Mơ hình thủy văn 1.1.1 Khái niệm vai trò 1.1.2 Phân loại 1.1.3 Tiêu chuẩn lựa chọn mơ hình 1.1.4 Một số mơ hình thủy văn phổ biến 1.2 Mơ hình thủy văn thơng số tập trung Nielsen-Hansen 13 1.2.1 Giới thiệu 13 1.2.2 Cấu trúc mơ hình 14 1.2.3 Các thông số mơ hình 16 CHƯƠNG II: PHƯƠNG PHÁP SỐ GIẢI BÀI TỐN ƯỚC TÍNH THƠNG SỐ TỐI ƯU CHO MƠ HÌNH THỦY VĂN THƠNG SỐ TẬP TRUNG 18 2.1 Tổng quan phương pháp giải tốn ước tính thơng số tối ưu 18 2.2 Phương pháp ô vuông 19 2.3 Phương pháp Rosenbrock 22 2.4 Phương pháp Nelder-Mead 26 2.5 Phương pháp Hooke-Jeeves 27 2.6 Giải thuật di truyền 28 2.7 Phương pháp SCE (Shuffled Complex Evolution) 30 CHƯƠNG III: XÂY DỰNG PHẦN MỀM DỰ BÁO LŨ BẰNG MƠ HÌNH THỦY VĂN THƠNG SỐ TẬP TRUNG CĨ SỬ DỤNG PHƯƠNG PHÁP ƯỚC TÍNH THƠNG SỐ TỐI ƯU 31 3.1 Sơ đồ khối mơ hình Nielsen-Hansen 3.2 Sơ đồ khối phương pháp số giải tốn ước tính thơng số tối ưu 3.2.1 Sơ đồ khối thuật tốn vuông 32 3.2.2 Sơ đồ khối thuật toán RosenBrock 34 ii 3.2.3 Sơ đồ khối thuật toán Nelder-Mead 35 3.2.4 Sơ đồ khối thuật toán Hooke-Jeeves 38 3.2.5 Sơ đồ khối giải thuật di truyền 39 3.2.6 Sơ đồ khối thuật toán SCE 41 3.3 Xây dựng phần mềm 44 CHƯƠNG IV: THỬ NGHIỆM CÔNG NGHỆ TRONG TÍNH TỐN DỰ BÁO LŨ TRÊN MỘT SỐ LƯU VỰC SÔNG 47 4.1 Thử nghiệm dự báo 47 4.2 Kết tính tốn đánh giá 49 KẾT LUẬN VÀ KIẾN NGHỊ 58 TÀI LIỆU THAM KHẢO 59 iii DANH MỤC KÝ HIỆU, CHỮ VIẾT TẮT THƯỜNG SỬ DỤNG NST Nhiễm sắc thể SCE Shuffled Complex Evolution WMO Tổ chức khí tượng giới NSE Nash–Sutcliffe Efficiency SCE Shuffled Complex Evolution KTTV Khí tượng – Thủy văn LV Lưu vực CSDL Cơ sở liệu BTNMT Bộ Tài nguyên môi trường KHCN Khoa học công nghệ VN Việt Nam iv DANH MỤC CÁC BẢNG Bảng 1.1: Bảng tiêu chuẩn lựa chọn mơ hình WMO Bảng 1.2: Các thơng số hiệu chỉnh mơ hình Nielsen-Hansen 17 Bảng 4.1: Bảng kết hệ số NSE sử dụng phương pháp tính tốn dị tìm thơng số tối ưu thử nghiệm cho 12 lưu vực 49 v DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ Hình 1.1: Sơ đồ phân loại mơ hình dịng chảy Hình 1.2: Sơ đồ mơ mơ hình Nielsen-Hansen 14 Hình 2.1: Mơ phương pháp giải vng 21 Hình 2.2: Các ánh xạ gương, phép co, phép dãn 26 Hình 3.1: Sơ đồ khối mơ hình Nielsen-Hansen 31 Hình 3.2: Sơ đồ khối thuật tốn vuông 32 Hình 3.3:Sơ đồ khối thuật tốn RosenBrock 34 Hình 3.4: Sơ đồ khối thuật toán Nelder-Mead 35 Hình 3.5: Sơ đồ khối thuật toán Hooke-Jeeves 38 Hình 3.6: Sơ đồ khối thuật toán giải thuật di truyền 40 Hình 3.7: Sơ đồ khối thuật tốn SCE 41 Hình 3.8: Sơ đồ khối thuật tốn CCE 43 Hình 3.9: Giao diện chung phần mềm 45 Hình 3.10: Thơng tin tham số phương án tính tốn 45 Hình 3.11: Một số phương pháp giải có phần mềm 46 Hình 3.12: Lựa chọn loại hình tính tốn cho mơ hình thủy văn 46 Hình 4.1: Vị trí lưu vực nghiên cứu 47 Hình 4.2: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Lâm Sơn 51 Hình 4.3: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực hồ Bản Chát 51 Hình 4.4: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực sơng Cầu 52 Hình 4.5: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Hàm Yên 52 Hình 4.6: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Hịa Bình 53 Hình 4.7: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Huội Quảng 53 Hình 4.8: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Sơn La 54 Hình 4.9: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Thác Bà 54 vi Hình 4.10: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực sông Thương 55 Hình 4.11: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Tuyên Quang 55 Hình 4.12: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Yên Bái 56 Hình 4.13: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Lục Nam 56 MỞ ĐẦU Tính cấp thiết đề tài Dự báo thuỷ văn, thủy lực công việc dự báo trước cho tương lai cách có khoa học trạng thái biến đổi yếu tố thuỷ văn, nhiên biến đổi trình tự nhiên phức tạp, chịu tác động nhiều yếu tố Tính biến động yếu tố phụ thuộc vào không gian thời gian nên gây khó khăn lớn cho q trình dự báo Thêm vào đó, thiếu trạm quan trắc cần thiết thiếu kết hợp ngành liên quan liệu quan trắc thực tế thường khơng đầy đủ, khơng mang tính chất đại diện Dự báo thủy văn có ý nghĩa lớn điều tiết dòng chảy, khai thác nguồn nước cho phát điện, giao thông thủy, tưới (nhất vùng khô hạn), cấp nước, quản lý nguồn nước, chất lượng nước Dự báo thủy văn đặc biệt quan trọng đối phó với tượng nguy hiểm sông lũ, lụt, hạn hán, quan trọng thiết kế thi cơng, vận hành cơng trình thủy lợi, cơng trình thủy nói chung Hiện nay, có nhiều phương pháp dự báo đưa dựa mơ hình vật lý tốn học, kết mơ hình nói đạt số thành công đáng ghi nhận Tuy nhiên, mơ hình có ưu nhược điểm định, vấn đề tìm kiếm phương pháp đủ tốt, đáp ứng yêu cầu thực tế giải toán dự báo thuỷ văn, thủy lực nội dung nghiên cứu thời Các mơ hình thủy văn thường bao gồm hàng chục thông số đặc trưng lưu vực Thông qua việc sử dụng mơ hình thủy văn với thơng số này, ta có chuỗi giá trị lưu lượng chảy Cần hiệu chỉnh thơng số mơ hình cho giá trị lưu lượng tính tốn gần với giá trị lưu lượng thực đo biết trước Trải qua nhiều lần hiệu chỉnh cho nhiều số liệu (nhiều năm), ta có thơng số phù hợp để sử dụng dự báo cho lưu vực Tuy nhiên, để hiệu chỉnh thông số lúc kinh nghiệm người dự báo viên khó khăn, nhiều cơng sức kết đầu (lưu lượng dòng chảy) hàm phi tuyến thơng số Ngồi ra, hàm quan hệ phi tuyến khơng thể đạo hàm nên cần tìm phương pháp để dị tìm tham số cần thiết Tìm kiếm phương pháp giải tốn ước tính thơng số tối ưu mơ hình thủy văn giúp giải tốn hiệu quả, tìm thông số phù hợp cho vùng phục vụ cho trình tác nghiệp dự báo thủy văn cách xác Trong đề tài này, tác giả chọn lựa hướng đề xuất “Nghiên cứu phát triển mơ hình thủy văn thơng số tập trung dự báo lũ cho lưu vực sông Việt Nam” Mục đích nghiên cứu Nghiên cứu xây dựng mơ hình thủy văn thơng số tập trung sử dụng phương pháp số giải tốn ước tính thơng số tối ưu Phương pháp nghiên cứu - Phân tích, thống kê tính tốn tư liệu thu thập - Mơ hình hóa số trị tính tốn mơ hình thủy văn thơng số tập trung có sử dụng ước tính tối ưu việc phát triển chương trình tính tốn ngơn ngữ Fortran, C# 46 Jeeves, phương pháp Nelder-Mead, phương pháp giải thuật di truyền phương pháp SCE Hình 3.11: Một số phương pháp giải có phần mềm Tiếp đến, ngồi phương pháp dị tìm thơng số thủ cơng người dùng có thêm lựa chọn khác dị tìm thơng số tự động Hộp thoại ComboBox giúp người dùng dễ dàng lựa chọn thay đổi kiểu tính tốn cho mơ hình thủy văn Hình 3.12: Lựa chọn loại hình tính tốn cho mơ hình thủy văn 47 CHƯƠNG IV THỬ NGHIỆM CƠNG NGHỆ TRONG TÍNH TỐN DỰ BÁO LŨ TRÊN MỘT SỐ LƯU VỰC SÔNG 4.1 Thử nghiệm dự báo Nhóm tác giả thử nghiệm số phương pháp giải tốn ước tính thơng số tối ưu cho mơ hình thủy văn, phục vụ dự báo lũ cho số lưu vực sơng Việt Nam – nhóm tác giả thành viên nhóm tư vấn cho Cục Phòng tránh thiên tai dự báo lũ lụt đồng sơng Hồng – Thái Bình hàng năm Trong nghiên cứu này, thử nghiệm 12 lưu vực thủy văn bao gồm: Bản Chát, Huội Quảng, Sơn La, Hịa Bình, Thác Bà, Tun Quang, Hàm Yên, Yên Bái, Cầu, Thương, Lục Nam, Lâm Sơn Trong 12 lưu vực trên, ngồi Lâm Sơn, 11 lưu vực lưu vực hệ thống lưu vực sơng Hồng – Thái Bình, hiệu chỉnh dị tìm thơng số với số liệu năm 2016 đến 2019 lấy thông số tốt kiểm định lại với số liệu năm 2020 Lâm Sơn lưu vực nhỏ bao gồm trạm mưa trạm đo lưu lượng, lựa chọn để thử nghiệm phương pháp ô vng, hiệu chỉnh dị tìm số liệu năm 2004 kiểm định với số liệu năm 2006 Hình 4.1: Vị trí lưu vực nghiên cứu - Lưu vực thủy văn Bản Chát có diện tích gần 1.929 km , phía giáp với biên giới Việt Nam – Trung Quốc, phía giáp với lưu vực thủy văn Huội Quảng Nước lưu vực dồn thủy điện Bản Chát, thủy điện bậc thang phía thủy điện Huội Quảng, lưu lượng xả hồ Bản Chát hồ Huội Quảng khoảng thời gian ngắn 48 - Lưu vực thủy văn Thác Bà có tổng diện tích 6.430 km , có nhánh dịng sông Chảy thuộc lưu vực sông Lô Lưu vực thủy văn Thác Bà lưu vực thủy văn độc lập, khơng chịu ảnh hưởng q trình điều tiết hồ thủy điện khác - Lưu vực thủy văn Tun Quang có tổng diện tích 5.868 km , có nhánh dịng sơng Gâm thuộc lưu vực sông Lô Lưu vực thủy văn Tuyên Quang lưu vực thủy văn độc lập, không chịu ảnh hưởng trình điều tiết hồ thủy điện khác - Lưu vực thủy văn Huội Quảng có diện tích gần 2.825 km2, phía giáp với đập thủy điện Bản Chát, phía giáp với lưu vực thủy văn Sơn La Lưu vực Huội Quảng chịu ảnh hưởng trực tiếp thủy điện Bản Chát gây ảnh hưởng tới lưu vực thủy văn Sơn La Cần thu thập số liệu điều tiết thủy điện Bản Chát để tính tốn thủy văn cho lưu vực hồ Huội Quảng - Lưu vực hồ Sơn La có diện tích gần 43.760 km2, phía giáp với đập thủy điện Lai Châu, phía giáp lưu vực hồ Hịa Bình Để tính tốn thủy văn cho lưu vực Sơn La cần thu thập số liệu điều tiết thủy điện Lai Châu thủy điện Huội Quảng thủy điện bậc thủy điện Sơn La Lưu vực Sơn La nhận lưu lượng nước trực tiếp xả từ thủy điện Lai Châu thủy điện Huội Quảng sau khoảng thời gian chảy truyền - Lưu vực hồ Hịa Bình có diện tích gần 53.600 km , phía giáp với đập thủy điện Sơn La, phía kết nối với hạ lưu sơng Đà qua đập thủy điện Hịa Bình Để tính tốn thủy văn cho lưu vực Hịa Bình cần phải thu thập số liệu điều tiết hồ chứa Sơn La thủy điện bậc thủy điện Hịa Bình lưu vực Hịa Bình nhận lưu lượng nước trực tiếp thủy điện Sơn La xả - Lưu vực Yên Bái có diện tích 27.306 km thuộc phần lưu vực sơng Thao Lưu vực sơng Thao dịng sông Hồng bắt nguồn từ dãy núi Ngụy Sơn cao 2000 m thuộc tỉnh Vân Nam Trung Quốc (ở Trung Quốc sơng Thao cịn có tên sơng Ngun) Sơng Thao có chiều dài sơng 902 km (tính đến Việt Trì) chiều dài nước 332 km - Lưu vực Hàm Yên có diện tích 11.860 km thuộc phần lưu vực sơng Lơ Lưu vực Hàm n tính từ thượng nguồn sông Lô – Hà Giang trải dài xuống Bắc Quang, Vĩnh Tuy Hàm Yên - Lưu vực sơng Cầu có diện tích khoảng 2.760 km có dịng sơng Cầu bắt nguồn từ vùng núi Tam Tao (1.326 m), phía tây tỉnh Bắc Kạn, chảy 49 theo hướng tây bắc - đông nam, qua thị xã Bắc Kạn, chuyển hướng gần bắc - nam chảy qua Chợ Mới, thành phố Thái Nguyên, Đáp Cầu tiếp nhận sông Thương thượng lưu Phả Lại khoảng km - Lưu vực sơng Thương có diện tích khoảng 2.330 km , có tên nhánh sơng sơng Thương có chiều dài 157 km bắt nguồn từ dãy núi Na-PaPhước cao 600 m huyện Chi Lăng tỉnh Lạng Sơn, chảy theo hướng đông bắc - tây nam qua thành phố Bắc Giang, tiếp nhận thêm sông Lục Nam chảy vào sông Cầu thượng lưu Phả Lại - Lưu vực sông Lục Nam có diện tích lưu vực vào khoảng 2.090 km Dịng sơng Lục Nam dài khoảng 175 km bắt nguồn từ vùng núi Kham, cao 700 m huyện Đình Lập tỉnh Lạng Sơn, chảy theo hướng đông bắc tây nam qua huyện Sơn Động, Lục Ngạn, Lục Nam tỉnh Bắc Giang đổ vào sông Thương Tứ Yên, cách Phả Lại 10 km phía thượng lưu - Lưu vực Lâm Sơn có diện tích lưu vực khoảng 33 km Dịng sơng Bùi có chiều dài 91 km, bắt nguồn từ xã Lâm Sơn, huyện Lương Sơn, Hịa Bình chảy qua tỉnh Hà Nội, Hồ Bình hợp lưu với sông Đáy xã Phúc Lâm, huyện Chương Mỹ, Hà Nội 4.2 Kết tính tốn đánh giá Ngồi phương pháp vng, phươnfg pháp khác lựa chọn ngẫu nhiên hay nhiều nghiệm ban đầu, việc lựa chọn dẫn đến cực trị địa phương, tác giả thực chạy phương pháp 1000 lần (100 lần với số lưu vực lớn) để tính tốn nhằm hạn chế việc kết rơi vào cực trị địa phương Trong bảng 4.1 kết hệ số NSE sau tính tốn thử nghiệm phương pháp số dị tìm thông số tối ưu cho lưu vực Bảng 4.1: Bảng kết hệ số NSE sử dụng phương pháp tính tốn dị tìm thơng số tối ưu thử nghiệm cho 12 lưu vực Bảng kết NSE tính tốn RosenBrock HC KĐ HookeJeeves HC KĐ HC KĐ HC SCE KĐ 0.64 Sơn Max 0.74 TGTB 5.6s 8.8s 6.3s 45.5s 29.3s Bản TBC 0.40 0.56 0.35 0.44 0.57 Chát Max 0.51 TGTB 187s TBC 0.28 Max 0.34 0.44 0.74 0.59 0.61 0.52 346s 0.19 0.25 0.28 0.75 0.49 198s 0.21 0.70 HC TBC 0.54 0.70 GA Lâm Huội Quảng 0.70 NelderMead 0.59 0.30 0.74 0.5 0.74 0.60 0.35 345s 0.30 0.36 KĐ 0.75 0.58 0.62 0.51 755s 0.25 0.50 0.52 0.5 50 TGTB Sơn La 201s 355s 0.50 Max 0.54 TGTB 1484s TBC 0.19 Bình Max 0.33 0.24 0.44 TGTB 404s 693s Thác TBC 0.28 Bà Max 0.51 TGTB 235s Tuyên TBC 0.14 Quang Max 0.35 TGTB 329s Yên TBC 0.34 Bái Max 0.47 TGTB 335s Hàm TBC 0.12 Yên Max 0.19 TGTB 135s TBC 0.43 Max 0.59 TGTB 85s TBC 0.4 Max 0.55 TGTB 213s Lục TBC 0.43 Nam Max 0.56 TGTB 139s Thương 724s TBC Hòa Cầu 337s 0.43 0.48 0.40 0.52 0.61 0.55 0.55 0.60 372s 0.56 0.42 0.51 0.49 534s 0.61 0.60 0.55 0.66 588s 0.37 0.54 0.3 0.57 262s 0.74 0.30 0.70 174s 0.56 0.53 0.16 198s 0.57 0.64 0.43 0.22 0.61 215s 0.70 0.73 0.59 194s 0.10 357s 0.5 0.50 0.50 95s 0.52 0.51 0.50 623s 0.64 0.53 0.58 1004s 0.35 0.2 0.45 987s 0.53 0.4 0.52 760s 0.50 0.22 0.50 1005s 0.57 0.4 0.51 0.4 0.56 0.35 335s 0.33 0.56 134s 281s 0.57 0.52 707s 0.44 0.39 0.65 493s 0.26 0.11 0.74 0.61 0.51 0.62 0.59 652s Trong đó: - TBC: giá trị trung bình cộng NSE nhiều lần tính - Max: Giá trị NSE lớn đạt - TGTB: Thời gian trung bình 01 lần tính - Các có giá trị trống: NSE < 0, phương pháp không hội tụ Đối với lưu vực Lâm Sơn, phương pháp ô vuông thử nghiệm nhiều lần thay đổi số ô lưới chia Khi ô lưới tăng từ đến 10, số NSE đạt 0.39 thời gian tính tốn lớn Cần phải chia ô lưới mịn để cực trị địa phương đủ rộng mắt lưới, nhiên ứng dụng thực tế khó thời gian tính tốn tăng lên nhiều lần 51 Từ hình 4.2 đến hình 4.13 đồ thị so sánh kết lưu lượng thực đo lưu lượng nước tính tốn kiểm định cho 12 lưu vực thủy văn lựa chọn để tính toán thử nghiệm 16 Lưu lượng Thực đo m /s Lưu lượng tính tốn 14 12 10 0:0011/9/2004 0:009/30/2004 0:009/10/2004 0:008/21/2004 0:008/1/2004 0:007/12/2004 0:006/22/2004 0:006/2/2004 0:005/13/2004 0:004/23/2004 0:0010/20/2004 Hình 4.2: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Lâm Sơn 6000 Lưu lượng thực đo m /s Lưu lượng tính tốn 5000 4000 3000 2000 1000 0:009/6/ 2020 0:08/27/ 2020 0:08/17/ 2020 0:008/7/ 2020 0:07/28/ 2020 0:07/18/ 2020 0:007/8/ 2020 0:06/28/ 2020 0:06/18/2 020 0:006/8/2 020 Hình 4.3: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực hồ Bản Chát 52 450 400 Lưu lượng thực đo m /s Lưu lượng tính tốn 350 300 250 200 150 100 50 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/20 20 Hình 4.4: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực sơng Cầu 1800 1600 Lưu lượng thực đo m /s Lưu lượng tính tốn 1400 1200 1000 800 600 400 200 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/202 Hình 4.5: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Hàm Yên 53 2000 1800 Lưu lượng thực đo m /s Lưu lượng tính tốn 1600 1400 1200 1000 800 600 400 200 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/202 Hình 4.6: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Hịa Bình 250 Lưu lượng thực đo m /s Lưu lượng tính tốn 200 150 100 50 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/202 Hình 4.7: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Huội Quảng 54 6000 Lưu lượng thực đo m /s Lưu lượng tính tốn 5000 4000 3000 2000 1000 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/202 Hình 4.8: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Sơn La 3500 Lưu lượng thực đo m /s Lưu lượng tính tốn 3000 2500 2000 1500 1000 500 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/202 Hình 4.9: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Thác Bà 55 300 Lưu lượng thực đo m /s Lưu lượng tính tốn 250 200 150 100 50 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/202 Hình 4.10: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực sông Thương 2000 1800 Lưu lượng thực đo m /s Lưu lượng tính tốn 1600 1400 1200 1000 800 600 400 200 0:009/6/2020 0:008/27/2 020 0:008/17/2 020 0:008/7/20 20 0:007/28/2 020 0:007/18/2 020 0:007/8/20 20 0:006/28/2 020 0:006/18/2 020 0:006/8/202 Hình 4.11: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Tuyên Quang 56 6000 Lưu lượng thực đo m /s Lưu lượng tính tốn 5000 4000 3000 2000 1000 0:008/17/2020 0:008/7/2020 0:007/28/2020 0:007/18/2020 0:007/8/2020 0:006/28/2020 0:006/18/2020 0:006/8/2020 0:009/6/2020 0:008/27/2020 Hình 4.12: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Yên Bái 1000 Lưu lượng thực đo m /s Lưu lượng tính tốn 900 800 700 600 500 400 300 200 100 0:009/6/2020 0:009/1/2 020 0:008/27/202 0:008/22/20 20 0:008/12/2020 0:008/7/2020 0:008/2/2020 0:007/28/2020 0:008/17/20 20 Hình 4.13: Đồ thị lưu lượng thực đo lưu lượng tính kiểm định lưu vực Lục Nam 57 Đánh giá kết quả: - Dựa tiêu chí thời gian tính: Phương pháp RosenBrock nhanh phương pháp SCE thường lâu - Dựa tiêu chí kết tốt nhất: Trong 12 lưu vực, phương pháp SCE thường cho kết tính tốn tốt (chỉ số NSE cao nhất), đặc biệt lưu vực tiếp nhận lượng nước từ lưu vực thượng nguồn đổ vào Huội Quảng, Sơn La, Hồ Bình, phương pháp khác không hội tụ (chỉ số NSE