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 Phịng chun mơn TS Nguyễn Chính Kiên Thủy động lực Giảm nhẹ thiên tai lưu vực ThS Dương Thị Thanh Hương Thủy động lực Giảm nhẹ thiên tai lưu vực ThS Nguyễn Tuấn Anh PGS TS Trần Thu Hà ThS Nguyễn Hồng Phong 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 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 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 .1 CHƯƠNG I: TỔNG QUAN VỀ MÔ HÌNH THỦY VĂN .3 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 .7 1.1.4 Một số mơ hình thủy văn phổ biến .9 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 .31 3.2 Sơ đồ khối phương pháp số giải tốn ước tính thơng số tối ưu .32 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 ô vuông .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 vng 32 Hình 3.3:Sơ đồ khối thuật toán RosenBrock .34 Hình 3.4: Sơ đồ khối thuật tố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 toán SCE 41 Hình 3.8: Sơ đồ khối thuật toá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 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 đạ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à, Tuyên 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, 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 km2, 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 km2, 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 km2, 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 km2, 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 n Bái có diện tích 27.306 km2 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 n có diện tích 11.860 km2 thuộc phần lưu vực sông Lô Lưu vực Hàm Yê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 km2 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 km2, 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-Pa-Phướ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 km2 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 km2 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 ô vuông, 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 toán Lâm Sơn Bản Chát Huội Quảng RosenBrock HC KĐ HookeJeeves HC KĐ 0.70 NelderMead HC KĐ HC SCE KĐ 0.64 Max 0.74 TGTB 5.6s 8.8s 6.3s 45.5s 29.3s TBC 0.40 0.56 0.35 0.44 0.57 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 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 337s 724s TBC 0.50 Max 0.54 TGTB 1484s TBC 0.19 Max 0.33 TGTB 404s 693s 1005s TBC 0.28 0.57 0.55 Max 0.51 TGTB 235s 372s 760s TBC 0.14 0.50 0.42 Max 0.35 TGTB 329s 534s 987s TBC 0.34 0.53 0.60 Max 0.47 TGTB 335s 588s 1004s TBC 0.12 0.35 0.54 Max 0.19 TGTB 135s 262s TBC 0.43 0.64 Max 0.59 TGTB 85s 174s 95s 194s 493s TBC 0.4 0.52 0.10 0.26 0.56 Thương Max 0.55 Hịa Bình Thác Bà Tun Quang n Bái Hàm Yên Cầu Lục Nam 0.51 0.43 0.24 0.4 0.22 0.4 0.2 0.53 0.51 0.44 0.61 0.56 0.61 0.37 0.74 0.56 0.48 0.40 0.52 0.55 0.60 0.51 0.49 0.55 0.66 0.3 0.57 0.30 0.70 0.53 0.50 0.16 0.50 0.45 0.58 0.50 0.43 0.11 0.70 0.4 0.73 0.59 0.35 0.74 0.57 213s 357s 198s 335s 707s TBC 0.43 0.57 0.22 0.44 0.61 Max 0.56 TGTB 139s 0.64 0.52 623s TGTB 0.5 0.50 0.61 215s 0.39 134s 0.33 0.56 0.51 281s 0.62 0.65 0.52 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 tốn thử nghiệm 16 m3/s Lưu lượng Thực đo Lưu lượng tính toán 14 12 10 11/9/2004 0:00 10/20/2004 0:00 9/30/2004 0:00 9/10/2004 0:00 8/21/2004 0:00 8/1/2004 0:00 7/12/2004 0:00 6/22/2004 0:00 6/2/2004 0:00 5/13/2004 0:00 4/23/2004 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 5000 4000 3000 2000 1000 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 Lưu lượng thực đo m3/s Lưu lượng tính tốn 400 350 300 250 200 150 100 50 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 1600 1400 1200 1000 800 600 400 200 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo 1800 Lưu lượng tính tốn 1600 1400 1200 1000 800 600 400 200 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 200 150 100 50 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 5000 4000 3000 2000 1000 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 3000 2500 2000 1500 1000 500 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 250 200 150 100 50 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 1800 1600 1400 1200 1000 800 600 400 200 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 5000 4000 3000 2000 1000 9/6/2020 0:00 8/27/2020 0:00 8/17/2020 0:00 8/7/2020 0:00 7/28/2020 0:00 7/18/2020 0:00 7/8/2020 0:00 6/28/2020 0:00 6/18/2020 0:00 6/8/2020 0:00 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 m3/s Lưu lượng thực đo Lưu lượng tính tốn 900 800 700 600 500 400 300 200 100 9/6/2020 0:00 9/1/2020 0:00 8/27/2020 0:00 8/22/2020 0:00 8/17/2020 0:00 8/12/2020 0:00 8/7/2020 0:00 8/2/2020 0:00 7/28/2020 0:00 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