CHƯƠNG 3 KẾT QUẢ THỰC NGHIỆM VÀ THẢO LUẬN
3.3. Thực nghiệm xác định hệ số a,b của mô hình Priestley – Taylor với điều kiện
3.3.3 Kết quả tính hệ số a,b của mô hình Priestley – Taylor với điều kiện địa hình,
khí hậu tỉnh Hòa Bình, Sơn La thuộc vùng Tây Bắc Việt Nam từ dữ liệu quan trắc khí tượng thủy văn
3.3.3.1. Kết quả xác định hệ số tuyến tính a, b của mô hình Priestley-Taylor
Sau khi tính được các tham số Rnd, λ, γ, Δ từ dữ liệu khí tượng quan trắc tại 6 trạm khí tượng thủy văn Hòa Bình và 2 trạm khí tượng Sơn La, kết hợp với lượng bốc thoát hơi nước đo được tại thời điểm ngày 01/7/2015, ngày 15/6/2016, ngày
04/6/2017, ngày 25/7/2018, ngày 09/7/2019, ngày 07/7/2020 và ngày 18/8/2021. Thay các giá trị Rnd, λ, γ, Δ và ETa_Đo vào mô hình Priestley – Taylor, tại mỗi thời điểm lập được hệ 8 phương trình dạng (2.40) với 2 ẩn số a, b giải bài toán theo nguyên lý số bình phương nhỏ nhất xác định được hệ số a, b cho các thời điểm ngày 01/7/2015, ngày 15/6/2016, ngày 04/6/2017, ngày 25/7/2018, ngày 09/7/2019, ngày 07/7/2020 và ngày 18/8/2021. Đề xuất hệ số a, b phù hợp với điều kiện địa hình khí hậu tỉnh Hòa Bình thuộc vùng Tây Bắc Việt Nam là giá trị trung bình của 7 thời điểm ngày 01/7/2015, ngày 15/6/2016, ngày 04/6/2017, ngày 25/7/2018, ngày 09/7/2019, ngày 07/7/2020 và ngày 18/8/2021. Kết quả tính toán được thể hiện trong bảng sau:
Bảng 3.14. Kết quả tính hệ số a, b của mô hình Priestley - Taylor phù hợp với điều kiện địa hình, khí hậu tỉnh Hòa Bình thuộc vùng Tây Bắc Việt Nam
Kết quả tính hệ số a, b của mô hình Priestley - Taylor theo chuỗi thời điểm
01/7/2015 15/6/2016 4/6/2017 25/7/2018 9/7/2019 07/7/2020 18/8/2021 Trung bình a b a b a b a b a b a b a b a b 0,795 - 0,087 0,771 - 0,118 0,862 0,123 0,767 - 0,081 0,835 - 0,171 0,782 0,123 0,730 0,029 0,792 - 0,026
Qua kết quả tính hệ số a, b của mô hình Priestley - Taylor với điều kiện địa hình, khí hậu tại tỉnh Hòa Bình thuộc vùng Tây Bắc Việt Nam từ dữ liệu khí tượng đo trực tiếp tại 6 trạm khí tượng thủy văn tỉnh Hòa Bình và 2 trạm khí tượng Sơn La tại thời điểm ngày 01/7/2015, ngày 15/6/2016, ngày 04/6/2017, ngày 25/7/2018, ngày 09/7/2019, ngày 07/7/2020 và ngày 18/8/2021 cho thấy giá trị a của mô hình Priestley - Taylor tính từ dữ liệu khí tượng thủy văn dao động từ 0,730 tại thời điểm ngày 18/8/2021 đến 0,862 tại thời điểm ngày 04/6/2017. Giá trị b của mô hình dao động từ -0,171 tại thời điểm ngày 09/7/2019 đến 0,123 tại thời điểm ngày 04/6/2017 và ngày 07/7/2020. Giá trị a, b trung bình tại các thời điểm là a = 0,792 và b = -0,026. Giá trị a, b trung bình sẽ được sử dụng để ước tính, giám sát lượng bốc thoát hơi nước trên địa bàn tỉnh Hòa Bình với các tham số Rnd, λ, γ, Δ được chiết xuất, tính toán từ ảnh vệ tinh và giá trị độ cao tính từ DEM.
3.3.3.2. Kết quả kiểm chứng xác định hệ số a, b của mô hình Priestley-Taylor
Sử dụng hệ số a, b xác định được từ các thời điểm ngày 01/7/2015, ngày 15/6/2016, ngày 04/6/2017, ngày 25/7/2018, ngày 09/7/2019, ngày 07/7/2020 và ngày 18/8/2021 để kiểm chứng kết quả tính tại 3 trạm (Thủy văn Hưng Thi, Thủy văn Lâm Sơn, Khí tượng Bắc Yên). Kết quả kiểm chứng được thể hiện trong bảng sau:
Bảng 3.15. Kết quả kiểm chứng xác định hệ số a, b của mô hình Priestley – Taylor
Tên trạm
Kết quả kiểm chứng tại các trạm (mm/ngày)
01/7/2015 15/6/2016 4/6/2017 25/7/2018 9/7/2019 07/7/2020 18/8/2021 ETa_ Đo ETa_ Ktra ETa_ Đo ETa_ Ktra ETa_ Đo ETa_ Ktra ETa_ Đo ETa_ Ktra ETa_ Đo ETa_ Ktra ETa_ Đo ETa_ Ktra ETa_ Đo ETa_ Ktra Thủy văn Hưng Thi 7,0 7,0 6,3 6,3 8,0 8,3 4,9 4,5 3,5 4,2 7,1 6,5 7,1 6,6 Thủy văn Lâm Sơn 6,8 7,2 6,5 6,3 9,2 8,6 4,8 4,4 5,1 4,4 6,8 6,4 6,2 6,6 Khí tượng Bắc Yên 7,2 7,3 7,1 6,8 7,1 8,1 4,7 4,8 3,1 3,9 7,8 6,7 6,5 6,5
Kết quả kiểm chứng xác định hệ số a,b của mô hình Priestley-Taylor tại 3 trạm (Thủy Văn Hưng Thi, Thủy Văn Lâm Sơn và Khí tượng Bắc Yên) cho thấy chênh lệch lượng bốc thoát hơi nước đo trực tiếp (ETa_Đo) và lượng bốc thoát hơi nước tính từ hệ số a,b theo mô hình Priestley-Taylor có giá trị chênh lệch lớn nhất là 1,1mm/ngày tại thời điểm ngày 07/7/2020 và chênh lệch trung bình tại 7 thời điểm ngày 01/7/2015, ngày 15/6/2016, ngày 04/6/2017, ngày 25/7/2018, ngày 09/7/2019, ngày 07/7/2020 và ngày 18/8/2021 là 0,4 mm/ngàỵ
3.4. Thực nghiệm tính lượng bốc thoát hơi nước thực tế từ bề mặt lớp phủ tại tỉnh Hòa Bình khu vực Tây Bắc Việt Nam sử dụng kết hợp mô hình SEBAL với