Quy trình hiệu chỉnh tự động đƣợc sử dụng để trợ giúp tối ƣu hóa chỉ tiêu NSE (hệ số Nash–Sutcliffe) trong khi giữ sai số tổng nƣớc tích lũy thấp. Do dòng chảy từ các lƣu vực cơ sở hạ lƣu và trung lƣu không đƣợc đo đạc một cách độc lập, hay nói cách khác trên lƣu vực sông Nhuệ Đáy chỉ có Ba Thá đƣợc xây dựng để thực hiện quan trắc lƣu lƣợng, vì thế không có hiệu chỉnh riêng của mô hình NAM đƣợc thực hiện cho từng lƣu vực cơ sở này. Bộ thông số hiệu chỉnh và kiểm định tại lƣu vực ND1 đƣợc áp dụng cho các lƣu vực còn lại.
Quá trình hiệu chỉnh tại Ba Thá trên lƣu vực sông Đáy từ tháng 01/1972 đến tháng 12/1974 cho thấy mức độ phù hợp tƣơng đối tốt giữa giá trị tính toán và thực đo dựa vào sai số tổng lƣợng, lƣu lƣợng đỉnh lũ, hệ số tƣơng quan và chỉ số Nash (hình 3. 3a). Sai số tổng lƣợng tích lũy là 0.3%, dòng chảy trung bình các tháng mùa lũ dƣới ƣớc tính 4.39% trong khi sai số đỉnh lũ vƣợt giá trị thực đo 23.33%. Hệ số tƣơng quan là 0.811 và chỉ số Nash đạt 0.824. Bộ thông số này đƣợc kiểm định với số liệu mƣa và lƣu lƣợng tại mặt cắt cửa ra của lƣu vực ND1 từ tháng 1 năm 1976 đến tháng 12 năm 1978 cũng cho thấy mức độ phù hợp tốt giữa kết quả tính toán và thực đo với sai số tổng lƣợng tích lũy là 0.4%, sai số đỉnh lũ dƣới ƣớc tính 8.23%, hệ số tƣơng quan 0.862 và hệ số Nash 0.858 (hình 3. 3b).
Nhận xét: Ngoại trừ một vài đỉnh lũ tính toán chênh lệch (lớn hơn hoặc nhỏ
hơn) so với quan trắc, nhìn chung các giai đoạn đều thể hiện mức độ phù hợp tƣơng đối tốt giữa chuỗi dòng chảy quan trắc và tính toán. Kết quả hiệu chỉnh và kiểm định đều đạt chỉ tiêu Nash khoảng 85%, thuộc loại khá, và hệ số tƣơng quan cao do đó có thể sử dụng bộ thông số trong mô hình vừa hiệu chỉnh và kiểm định tại trạm Ba Thá (bảng 3. 3) để tính toán dòng chảy đến cho các lƣu vực con thuộc lƣu vực sông Nhuệ Đáy (ND). Đồng thời qua giai đoạn hiệu chỉnh và kiểm định cho thấy giả thiết về điều kiện lƣu vực không thay đổi đáng kể trong tƣơng lai là hoàn toàn chấp nhận đƣợc.
51
RunOff, BASIN (1640.31), 1640.310 RunOff
R2 = 0.811 5-2-1972 23-8-1972 11-3-1973 27-9-1973 15-4-1974 1-11-1974 0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0 450.0 [m^3/s] 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 [] Model Series Data Series Low er Calc Threshold Model Series Data Series Low er Calc Threshold
RunOff, BASIN (1640.31), 1640.310 RunOff
R2 = 0.862 9-2-1975 20-5-1975 28-8-1975 6-12-1975 15-3-1976 23-6-1976 1-10-1976 9-1-1977 19-4-1977 28-7-1977 5-11-1977 0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0 450.0 [m^3/s] 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 [] Model Series Data Series Low er Calc Threshold Model Series Data Series Low er Calc Threshold
RunOff, BASIN (1640.31), 1640.310 RunOff
R2 = 0.811 50.0 100.0 150.0 200.0 250.0 300.0 350.0 [m^3/s] Gauge 0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0 450.0 [m^3/s] Model
Scatter Plot. Model vs Gauge a=0.836 b=11.5 Theoretical
Best Fit Line
Scatter Plot. Model vs Gauge a=0.836 b=11.5 Theoretical
Best Fit Line
RunOff, BASIN (1640.31), 1640.310 RunOff
R2 = 0.862 0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0 [m^3/s] Gauge 0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0 [m^3/s] Model
Scatter Plot. Model vs Gauge a=0.892 b=5.85 Theoretical
Best Fit Line
Scatter Plot. Model vs Gauge a=0.892 b=5.85 Theoretical
Best Fit Line
Hình 3. 3a. Hiệu chỉnh mô hình NAM xây dựng cho tiểu lƣu vực ND1 - Ba Thá (1972-1974). Sai số tổng lƣợng tích lũy = 0.3%, và NSE=0.824.
Hình 3. 3b. Kiểm định mô hình NAM xây dựng cho tiểu lƣu vực ND1 - Ba Thá cho (1976-1978). Sai số cân bằng nƣớc tích lũy = 0.4%, và NSE=0.858.
52
Bảng 3. 3. Bộ thông số tối ƣu cho lƣu vực sông Nhuệ Đáy
Thông số
Umax Lmax CQOF CKIF CK1.2 TOF TIF TG CKBF CK2
14 292 0.568 282.5 88.9 0.0374 0.141 0.0541 1318 10
U L QOF QIF BF Điều kiện
ban đầu
0.6 0.4 5 1 8