Trong nghiên cứu này đề xuất một quy trình các bước tính ứng dụng hệ các mô hình diễn toán SWAT, NAM và MIKE, kết hợp với các dữ liệu đo đạc thực tế, cùng các bước hiệu chỉnh và kiểm định xây dựng bộ hệ số thủy văn và thủy lực phục vụ cho tính toán dòng chảy sông Vệ, Quảng Ngãi.
BÀI BÁO KHOA HỌC ỨNG DỤNG MƠ HÌNH DIỄN TỐN SWAT/NAM/MIKE XÂY DỰNG BỘ THÔNG SỐ THỦY VĂN VÀ THỦY LỰC PHỤC VỤ CHO TÍNH TỐN DỊNG CHẢY TRƯỜNG HỢP SÔNG VỆ, QUẢNG NGÃI Lê Thị Mỹ Diệp1, Bùi Huỳnh Anh2, Bùi Tá Long2* Tóm tắt: Khu vực duyên hải tỉnh Quảng Ngãi có tốc độ tăng trưởng kinh tế nhanh Trước ảnh hưởng biến đổi khí hậu (BĐKH), tình trạng ngập lụt xâm nhập mặn diễn biến phức tạp Cả hai loại hình u cầu phải tính tốn dự báo diễn biến dịng chảy, vậy đới tượng nghiên cứu nhiều tác giả Trong nghiên cứu đềxuất qui trình bước tính ứng dụng hệ mơ hình diễn tốn SWAT, NAM MIKE, kết hợp với liệu đo đạc thực tế, bước hiệu chỉnh kiểm định xây dựng bộhệsốthủy văn thủy lực phục vụcho tính tốn dịng chảy sơng Vệ, Quảng Ngãi Kết hiệu chỉnh kiểm định được thực dự sốthống kê được sử dụng rộng rãi thủy văn cho phép khẳng định khả áp dụng thực tếlà chấp nhận được Điểm nghiên cứu so với nghiên cứu trước được thểhiện sơ đờtích hợp, thể rõ kết hợp số liệu, mơ hình bước hiệu chỉnh thông số thủy văn thủy lực Từ khóa: Mơ hình mưa - dịng chảy, SWAT, NAM, MIKE, sông Vệ Ban Biên tập nhận bài: 09/04/2019 Ngày phản biện xong: 20/6/2019 Đặt vấn đề Vùng duyên hải tỉnh Quảng Ngãi bao gồm huyện Bình Sơn, Sơn Tịnh, Tư Nghĩa, Mộ Đức, Đức Phổ Tp.Quảng Ngãi Nhờ thiên nhiên ưu đãi với vị trí địa lý thuận lợi mà huyện có tốc độ tăng trưởng kinh tế nhanh Tuy nhiên, năm qua trước ảnh hưởng BĐKH, tình trạng ngập lụt xâm nhập mặn khu vực diễn biến phức tạp, đối tượng nghiên cứu nhiều đề tài, dựán cấp Trong xây dựng công nghệ dự báo lũ lớn cảnh báo ngập lụt hệ thống sông Vệ - Trà Khúc khuôn khổ đề tài cấp nhà nước [1] Quá trình động lực diễn biến hình thái cửa Lở thực [3-4] Nghiên cứu dịng chảy ảnh hưởng tiêu lũ tỉnh Quảng Ngãi xây dựng tuyến đường cao tốc Đà Nẵng - Quảng Ngãi được thực Viện Khoa học Thủy lợi Miền Nam Trường Đại học Bách Khoa Tp HCM Email: longbt62@hcmut.edu.vn TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 Ngày đăng bài: 25/06/2019 [5] Nghiên cứu mơ hình dự báo thủy văn sông Trà Khúc thực [6] Đặc biệt, nghiên cứu [2] nhóm tác giả thực xác định thông số thủy văn thủy lực cho lưu vực sông Trà Khúc sông Vệ Một điểm chung nghiên cứu ứng dụng nhiều mơ hình thủy văn khác NAM, MARINE, SWAT, Trong tính tốn thủy văn, mơ hình “mưa rào - dịng chảy” thường sử dụng rộng rãi dù rằng, mơ hình hóa ln nhiệm vụ phức tạp trình diễn hệ thống thủy văn phức tạp, chưa biết đầy đủ [7] Mợt sớmơ hình dịng chảy “mưa rào - dòng chảy” được nghiên cứu nghiên cứu tác giả nước [8-14] Tại Việt Nam, có nhiều nghiên cứu áp dụng NAM dự báo dịng chảy từ mưa cho lưu vực sơng khác Nghiên cứu [15] ứng dụng NAM dự báo dịng chảy cho sơng Ba, Phú n phục vụ đánh giá tài nguyên nước cho lưu vực Nghiên cứu nhóm tác giả [16-17] ứng dụng BÀI BÁO KHOA HỌC NAM dự báo dịng chảy lưu vực sơng Cả, Bình Định khu vực hạ lưu sơng Sài Gịn Theo nghiên cứu [18], thập kỷ qua, hàng trăm mơ hình “mưa rào - dịng chảy” đề xuất với mục đích ứng dụng thời gian thực Các mơ hình phân loại theo quy mơ thời gian - khơng gian, theo phương pháp giải phương trình sử dụng Các tính để phân biệt cách tiếp cận là: thuật toán (dựa kinh nghiệm, khái niệm dựa trình); cách tiếp cận: ngẫu nhiên tiền định [7] Trong nghiên cứu này, ứng dụng hệ thống mơ hình diễn tốn SWAT/NAM/MIKE21 HD nhằm xác định thơng sớthủy văn thủy lực đạt yêu cầu khu vực cửa sông Vệ - bước cho nghiên cứu mô xâm nhập mặn khu vực cửa sông Vệ, khuôn khổ đề tài nghiên cứu sinh Phương pháp nghiên cứu thu thập tài liệu 2.1 Giới thiệu khu vực nghiên cứu Sông Vệ bắt nguồn từ rừng núi phía Tây, chảy theo hướng Tây Nam - Đông Bắc, huyện Tư Nghĩa, đổ biển Đông cửa Cổ Lũy cửa Đức Lợi Phía Bắc phía Tây giáp sơng Trà Khúc, phía Nam giáp tỉnh Bình Định phía Đơng giáp biển Lưu vực sơng Vệ có diện tích khoảng 1260km2, sơng có chiều dài khoảng 90km có 2/3 chiều dài chảy vùng núi có độ cao 100 - 1000m, mật độ sông suối lưu vực đạt 0,79km/km2, độ dốc bình quân lưu vực khoảng 19,9% [19] Thực vật che phủ bề mặt lưu vực vùng thượng lưu phần lớn rừng già, bụi rậm, vùng hạ lưu chủ yếu vùng đất canh tác nông nghiệp Hình thểhiện vị trí phạm vi nghiên cứu 2.2 Mơ hình SWAT, MIKE NAM, MIKE 21 HD 2.2.1 Mơ hình SWAT SWAT (Soil and Water Assessment Tool) công cụ đánh giá nước đất, xây dựng Trung tâm phục vụ nghiên cứu nông nghiệp (ARS - Agricultural Research Service) Mơ hình giúp đánh giá dự đoán tác động sử dụng đất đai tác động đến nguồn nước, phù sa, lượng hóa chất nơng nghiệp sinh lưu vực rộng lớn phức tạp Mơ hình mơ nhiều trình vật lý lúc lưu vực cho phép mô với mức độ chi tiết hóa cao cách chia lưu vực thành tiểu lưu vực theo địa hình mặng lưới thủy văn, sau tiểu lưu vực phân chia thành đơn vị thủy văn (Hydrological response units - HRUs) có đặc trưng riêng đất sử dụng đất dựa vào loại đất lớp phủ thực vật bên tiểu lưu vực Sự phân chia giúp người sử dụng áp dụng kết nghiên cứu vùng vào vùng khác chúng có tương đồng định Trong nghiên cứu này, SWAT kết hợp với ArcGIS được sử dụng đểphân lưu vực, phân tích đơn vị thủy văn, biên tập liệu khí tượng (vận tốc gió, lượng mưa, nhiệt độ max/min, đợẩm, cường độ xạ) Kết chạy SWAT nghiên cứu thể [19] Một số kết tương tự ứng dụng SWAT nhóm tác giả thể Hình Vị trí, giới hạn khu vực nghiên cứu với 12 tiểu lưu vực [19] (1) TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 BÀI BÁO KHOA HỌC công bố khoa học [20] 2.2.2 Mô hình thủy văn NAM Trong nghiên cứu sử dụng hệ thớng mơ hình thủy lực thủy văn nên trước tiên tổng quan công cụ được sử dụng Hệ thống dự báo thủy văn - thủy lực Mike11 gồm modules: (1) Mưa rào - dòng chảy (Rain - Runnoff hay ký hiệu RR), (2) Thủy lực (Hydrodynamic - HD), (3) Quy trình Flood Forecasting (FF) Mơ hình NAM (viết tắt cụm từ Nedbør-Afstrømnings-Model) phần hệ thống mơ hình mơ kênh sơng Mike11 Mơ hình có tên gọi “mưa rào dịng chảy” mơ hình tập trung gồm nhiều module con, đóng gói, mơ dòng chảy mặt đất, dòng kết nối với nhau, với dịng chảy Module RR áp dụng độc lập sử dụng để thể nhiều lưu vực gia nhập vào mạng lưới sông Module FF thiết kế để dự báo thay đổi mực nước lưu lượng hệ thống kênh sơng q trình “mưa rào - dịng chảy“ Mơ hình NAM mơ hình cải tiến từ mơ hình Nielsen-Hansen, cơng bố tạp chí “Nordic Hydrology” năm 1973, sau Viện Thủy lực Đan Mạch phát triển đổi thành NAM (từ từ viết tắt tiếng Đan Mạch: NedbørAfstrømnings-Model) Mơ hình gồm bể chứa, ngun lý tính tốn bể chứa giải phương trình cân theo quy luật phi tuyến (dạng đường cong nước rút) Mơ hình NAM mơ q trình mưa - dịng chảy cách liên tục thơng qua việc tính tốn cân nước bốn bể chứa thẳng đứng, có tác dụng qua lại lẫn để diễn tả tính chất vật lý lưu vực Các bể chứa gồm: Bể tuyết (chỉ áp dụng cho vùng có tuyết); Bể mặt; Bể sát mặt hay bể tầng rễ cây; Bể ngầm Dữ liệu đầu vào mơ hình mưa, bốc tiềm năng, nhiệt độ Kết đầu mơ hình dịng chảy lưu vực, mực nước ngầm, thơng tin khác chu trình thuỷ văn, thay đổi tạm thời độ ẩm đất khả bổ xung nước ngầm Dòng chảy lưu vực phân cách gần thành dòng chảy mặt, dòng chảy sát mặt, dòng chảy ngầm Mơ hình NAM thuộc loại mơ hình tất định, thơng số tập trung, mơ hình mơ liên tục Trong nghiên cứu này, NAM được sử dụng để tạo biên thủy lực cho mơ hình thủy lực MIKE21 HD Một số kết thuộc mực thểhiện [19] 2.2.3 Mơ hình thủy lực Mike 21 HD MIKE 21 HD DHI [27] thực phần mềm kỹ thuật để tính tốn dịng chảy, sông, hồ, cửa sông, vịnh, vùng biển ven bờ ngồi khơi Phương trình tốn MIKE21 HD được viết cho dịng hai chiều khơng đờng nhất lớp chất lỏng đồng theo phương thẳng đứng dựa phương trình Navier-Stokes, Reynolds khơng nén chiều với giả định Boussinesq áp suất thủy tĩnh Phương trình liên tục: ∂h ∂hU ∂hV + + = hS ∂t ∂x ∂y (1) Các phương trình động lượng theo phương ngang gồm ∂h ∂t + ∂hU ∂x + ∂hUV ∂y = fVh − gh + ∂ ∂x ∂η ∂x − h ∂p a ρ0 ∂x ( hTxx ) + ∂ ∂y − gh ∂ρ 2ρ0 ∂x + (1) τ xx ρ0 − τbx ρ0 − ∂Sxy + ρ0 ∂x ∂x ∂Sxx ( hT ) + hU S xy (2) s ∂hV ∂hV ∂hUV ∂η h ∂p a gh ∂ρ τxx τbx ∂Sxx ∂Sxy + + = −fUh − gh − − + − − + ∂t ∂x ∂y ∂y ρ0 ∂y 2ρ0 ∂y ρ0 ρ0 ρ0 ∂x ∂x ∂ ∂ + ( hTxy ) + ( hTyy ) + hVsS ∂x ∂y TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 (1) (3) BÀI BÁO KHOA HỌC Trong U, V thành phần vận tốc trung bình theo độ sâu thành phần vận tốc u, v theo phương tọa độ x, y; t thời gian; h = η + d: với mực nước, d độsâu Với Txx, Txy, Tyy thành phần ứng suất nhớt tổng cộng nhớt túy, nhớt rối khuếch tán S thành phần lưu lượng từ nguồn thải điểm; Us, Vs thành phần tốc độ nguồn thải nguồn điểm; g gia tốc trọng trường; f = 2Ω sin φ tham số Coriolis; ρ mật độ nước, ρ0 mật độ tiêu chuẩn; pa áp suất khí Hệthớng các phương trình được giải theo thuật tốn thể tích hữu hạn với lưới mềm dẻo 2.3 Thu thập tài liệu thiết lập mơ hình 2.3.1 Dữ liệu lưu vực Dữ liệu DEM lấy từ trang web http://gdex.cr.usgs.gov/gdex/ sử dụng phần mềm SWAT phân chia khu vực nghiên cứu thành lưu vực, sử dụng ArcGIS số hóa đoạn sơng Sau phân chia lưu vực nhận số liệu dạng shapefile tiểu lưu vực Các kết được trình bày [19] 2.3.2 Dữ liệu địa hình chạy MIKE Dữ liệu chạy MIKE 21 HD nghiên cứu chia làm nhóm Nhóm thứ liên quan tới vùng ven biển thu thập, xử lý chuyển vào phần mềm MIKE21 [22], nhóm thứ hai liên quan tới đất liền cửa sông gồm: số liệu thực đo 28 mặt cắt thực đo được kếthừa từ dự án trước [23-24] Đoạn sông Vệ được xem xét nghiên cứu giới hạn từ thượng nguồn sông đến cửa Lở dài 21,47km (Hình 3) Hình Địa hình biển được lựa chọn TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 BÀI BÁO KHOA HỌC Hình Địa hình vùng cửa sơng mặt cắt sông Vệ Trong nội dung chạy MIKE21 HD thuộc nghiên cứu sử dụng liệu biên mơ hình thủy lực vùng ven bờ biển lấy từ công cụ Tide Prediction of Height MIKE 21 Toolbox (.21t) cho năm 2015, 2017, 2018, nguồn [27] Các liệu sử dụng để chạy MIKE21 HD cho vùng ven biển nằm phạm vi nghiên cứu Hình 2.3.5 Chỉ tiêu đánh giá: Chỉ số đánh giá mức độ tương quan kết tính tốn kết đo đạc, xác định theo cơng thức sau: 2.3.3 Dữ liệu khí tượng - thủy văn Thơng tin số liệu khí tượng từ 02 trạm đo khí tượng thuộc tỉnh Quảng Ngãi Ba Tơ Quảng Ngãi năm 2013 - 2015, thông tin số liệu thủy văn thuộc trạm An Chỉtrong năm 2013 - 2015 được sử dụng từ nguồn [25] Các số liệu có dạng file thống kê thơng dụng xử lý u cầu mơ hình sử dụng Thông tin số liệu hải văn từ trạm Tam Quan, Bình Định được ssử dụng để thực hiện nghiên cứu lấy từ tài liệu [26] 2.3.4 Dữ liệu biên thủy lực n NSE = ∑ (Q sim i i=1 n ∑ (Q i=1 - Qiobs ) obs i -Q ) 2 , n PBIAS = RMSE = STDEVobs Tiêu chuẩn đánh giá theo số được thể hiện bảng 1, nghiên cứu [28] được sử dụng nghiên cứu 2.3.6 Thiết lập mơ hình tạo bộ̣thơng sớthủy văn Bộ thông số NAM sử dụng nghiên cứu gồm: Lmax: Lượng ẩm lớn bể chứa tầng rễ cây; Umax: Lượng nước tối đa bể trữ mặt, lượng nước để điền trũng đọng TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 obs i i=1 n - Qisim ) ×100 ∑ (Q ) obs i i=1 n RSR = ∑ (Q ∑ (Q obs i i=1 ∑ (Q n - Qisim ) obs i −Q ) 2 mặt thực vật chứa vài cm bề mặt đất; CQOF: Hệ số dòng chảy mặt (0≤ CQOF≤1), định phân phối mưa hiệu cho dòng chảy mặt thấm; CKIF: Hằng số thời gian dòng chảy sát mặt với Umax định dòng chảy sát mặt Nó chi phối thơng số diễn tốn dịng chảy trao đổi CKIF>>CK12, CK1,2: số thời gian cho diễn tốn dịng chảy mặt dịng chảy i=1 sát mặt dọc theo sườn dốc lưu vực qua lòng dẫn đến cửa lưu vực Ký hiệu CKBF: số thời gian dòng chảy ngầm, để diễn tốn hồn trả nước thơng qua lượng trữ nước ngầm, thường CKBF>>CK12 Ký hiệu TOF: giá trị ngưỡng dòng chảy mặt (0≤TOF≤1) Ký hiệu TIF: giá trị ngưỡng dòng chảy sát mặt (0≤TIF≤1) Dòng chảy sát mặt hình thành số ẩm tương đối đất tầng rễ lớn TIF Ký hiệu TG: giá trị ngưỡng lượng nước bổ sung cho dòng chảy ngầm (0≤TG≤1) Lượng nước bổ sung cho bể chứa ngầm hình thành số ẩm tương đối đất tầng rễ lớn TG [27] Các thông sốnày dạng bảng thể Bảng Hình Các bước thiết lập thông sốthủy văn cho lưu vực Các bước thiết lập mơ hình lập bợthơng sớ thủy văn được thực sau: liệu đầu vào gồm đồ địa hình (DEM), số liệu khí tượng - thủy văn, đồ sử dụng đất, đồ thổ nhưỡng đưa vào SWAT Kết triết xuất sau chạy SWAT gồm diện tích tiểu lưu vực kết hợp với số liệu mưa trạm khu vực đưa vào NAM (Hình 4) để thực bước hiệu chỉnh tham số mơ hình NAM Sử dụng tiêu mục 2.3.5 để đánh giá độ tin cậy thông số Để thực hiệu chỉnh, trước tiên, từ thông sốban đầu, NAM thực hiệu chỉnh tựđộng phương pháp thử dần để tăng tốc độ xác tới mức ổn định với sai số cho phép Bộ thông số sau hiệu chỉnh NAM sử dụng để tính tốn biên lưu lượng cho mơ hình thủy lực Một sớkết nội dung trình bày [19] 2.3.7 Thiết lập mơ hình tạo biên cho mơ hình thủy lực MIKE21 HD Biên xác định tiểu lưu vực số có lưu ý tới nhánh sông phụ chảy vào nhánh sông tiểu lưu vực 8, 10, 11, 12 (Hình 1) Từ mơ hình SWAT tính tốn tổng diện tích lưu vực Số liệu mưa bốc ngày lấy trạm khí tượng Ba Tơ Do điểm hợp lưu nằm vùng hạ lưu sơng Vệ nên mưa lấy từ trạm khí tượng Quảng Ngãi bốc tính tương ứng với tiểu lưu vực bằng sự hỗ trợ từ SWAT Từ thông số NAM hiệu chỉnh kiểm định mục 2.3.6, với dạng chuỗi theo thời gian lượng mưa, bốc vào mô hình NAM Các bước tạo file biên lưu lượng thượng lưu được thểhiện Hình TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 BÀI BÁO KHOA HỌC Hình Các bước thiết lập hệ thống mơ hình diễn tốn SWAT, NAM tạo biên lưu lượng cho mơ hình thủy lực Sử dụng liệu biên cho mơ hình thủy lực vùng ven bờ biển (Hình 2) lấy từ cơng cụ Tide Prediction of Height MIKE 21 Toolbox(.21t) cho năm 2015, nguồn [27] Thực hiện bước hiệu chỉnh kiểm mơ hình thủy lực MIKE21 HD Sau sử dụng liệu biên cho mơ hình thủy lực vùng ven bờ biển từ công cụ Tide Prediction of Height MIKE 21Toolbox(.21t) cho năm cần tính để xuất kết mực nước vị trí cửa Lở, sơng Vệ Các bước thiết lập thểhiện Hình Hình Các bước thiết lập hệ thớng mơ hình diễn tốn MIKE21 HD tạo mực nước cho mơ hình thủy lực TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 BÀI BÁO KHOA HỌC 2.3.8 Thiết lập mô hình hiểu chỉnh kiểm định thơng sớthủy lực Trong nghiên cứu này, hệ số nhám manning hệ số nhớt rối lựa chọn làm thông sốthủy lực Từ thông số thủy văn hiệu chỉnh kiểm định mục 2.3.6, với biên lưu lượng mực nước tạo từ 2.3.7 giúp vận hành MIKE 21 HD cho sông Vệ Dựa số liệu thực đo thủy văn tại trạm An Chỉ để hiệu chỉnh kiểm định thông sốthủy lực Kết thảo luận 3.1 Hiệu chỉnh kiểm định mơ hình NAM Bộ thủy văn thực đo trạm An Chỉ, giai đoạn 2013-2015 sử dụng để hiệu chỉnh kiểm định NAM Đểđánh giá mức độ tương quan, nghiên cứu sử dụng số Nash (các nghiên cứu tương tự được thực hiện [6], [15], [16]) Bộ số liệu thực đo lưu lượng năm 2013 sử dụng đểhiệu chỉnh đạt số NASH 92% Bước kiểm định thực theo sốliệu đo lưu lượng năm 2014 2015 đạt số NASH 90% 93% Kết thông số chọn để mô lưu lượng thể Bảng Kết hiệu chỉnh kiểm định thể Hình - Hình Hình Biểu đờq trình lưu lượng thực đo tính tốn theo NAM năm 2013 Hình Biểu đồ q trình lưu lượng thực đo tính tốn theo NAM năm 2014 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 BÀI BÁO KHOA HỌC Hình Biểu đồ q trình lưu lượng thực đo tính toán theo NAM năm 2015 Bảng Bảng kết thơng số mơ hình MIKE NAM Thơng số Umax Lmax CQOF TOF TIF TG CKIF CK12 CKBF Ý nghĩa GiỈ trị Lượng nước tối đa bể chứa mặt (mm) Lượng ẩm lớn bể chứa tầng rễ (mm) Hệ số dịng chảy mặt, khơng thứ nguyŒn, phản ánh điều kiện thấm Ngưỡng dòng chảy tràn Ngưỡng dịng chảy sỈt mặt GiỈ trị ngưỡng tầng rễ Hệ số thời gian dòng chảy sỈt mặt Hằng số thời gian chảy truyền dịng chảy mặt Hằng số thời gian chảy truyền dòng chảy ngầm Kết bảng có tương đồng với kết nhóm tác giả Linh cs (2018) [2], sử dụng cách tiếp cận khác Bộ thông số đảm bảo độ xác tính ổn định để sử dụng mơ dịng chảy thượng lưu cho năm khác thông số dùng để tính tốn dịng chảy phụ lưu khu làm biên đầu vào cho mơ hình dịng chảy lan truyền chất 3.2 Hiệu chỉnh, kiểm định Mike 21 HD 17 172 0,185 0,531 0,114 0,404 655,8 19 3972 vùng ven bờ biển Chuỗi số liệu biên lựa chọn sau: thời gian chọn để hiệu chỉnh: 1/1/2015 28/2/2015; thời gian chọn cho kiểm định: 1/3/2015 - 31/5/2015, mục 2.3.7 Dữ liệu thực đo mực nước trạm Tam Quan sử dụng để hiệu chỉnh kiểm định [26] Kết hiệu chỉnh kiểm định theo mục 2.3.5 được thểhiện Bảng Hình 10 Bảng Kết kiểm định hiệu chỉnh mô hình thủy lực MIKE21 HD cho vùng ven biển Hệ số tương quan Hệ số Nash PBIAS (%) RSR Hiệu Chỉnh 0,916 0,988 7,323 0.460 Kiểm Định 0,919 0,991 7,933 0.490 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2019 Kết quả, thông số lựa chọn sau: hệ số nhớt 0,28 (m2/s), hệ số nhám 30 (m1/3/s) 3.3 Hiệu chỉnh, kiểm định vùng cửa sông Vệ Chuỗi số liệu biên lựa chọn để chạy MIKE21 HD sau: chuỗi thời gian cho hiệu chỉnh: 1/4/2017 - 30/4/2017; thời gian chọn cho kiểm định: 1/4/2018 - 30/4/2018 Biên (lưu lượng) trích xuất từ kết chạy NAM cho phạm vi từ 1-30/4/2017 1-30/4/2018 Biên đước (mực nước) triết xuất từ kết chạy BÀI BÁO KHOA HỌC thủy lực cho khu vực biển mô tả mục 2.3.5 báo Dữ liệu thực đo dùng để hiệu chỉnh kiểm định mơ hình liệu dạng chuỗi thời gian đo mực nước trạm Sông Vệ [25] (mỗi bước giờ) Kết mục chọn thông số hệ số nhớt 0.28 hệ số nhám Manning 32m1/3/s Kết hiệu chỉnh kiểm định theo mục 2.3.5 được thể hiên Bảng Hình 10, dựa [28] có thểkết luận đảm bảo độ tin cậy Bảng Kết kiểm định hiệu chỉnh mơ hình thủy lực MIKE21 HD sông Vệ Hệ số tương quan Hệ số Nash PBIAS < 10% RSR