Nghiên cứu ứng dụng mô hình marine dự báo lũ sông kỳ lộ, tỉnh phú yên

BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI LUẬN VĂN THẠC SĨ NGHIÊN CỨU ỨNG DỤNG MƠ HÌNH MARINE DỰ BÁO LŨ SÔNG KỲ LỘ TỈNH PHÚ YÊN CHUYÊN NGÀNH: THỦY VĂN HỌC VÕ ANH KIỆT HÀ NỘI, NĂM 2018 BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI LUẬN VĂN THẠC SĨ NGHIÊN CỨU ỨNG DỤNG MƠ HÌNH MARINE DỰ BÁO LŨ SƠNG KỲ LỘ TỈNH PHÚ YÊN VÕ ANH KIỆT CHUYÊN NGÀNH: THỦY VĂN HỌC Mà SỐ: 60440224 NGƯỜI HƯỚNG DẪN KHOA HỌC: TS ĐẶNG THANH MAI TS TRƯƠNG VÂN ANH HÀ NỘI, NĂM 2018 CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI Cán hướng dẫn chính: TS Đặng Thanh Mai Cán chấm phản biện 1: TS Nguyễn Viết Thi Cán chấm phản biện 2: TS Nguyễn Lan Châu Luận văn thạc sĩ bảo vệ tại: HỘI ĐỒNG CHẤM LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI Ngày 27 tháng 12 năm 2017 LỜI CẢM ƠN Học viên xin chân thành cảm ơn Khoa Khí tượng Thủy văn, Phòng Đào tạo, Trường Đại học Tài nguyên Môi trường Hà Nội tạo điều kiện thuận lợi để tác giả học tập nghiên cứu thực luận văn Học viên bày tỏ lòng biết ơn sâu sắc đến TS Đặng Thanh Mai, Trung tâm dự báo Khí tượng Thủy văn Trung ương TS Trương Vân Anh, Trường Đại học Tài nguyên Môi trường Hà Nội, tận tnh hướng dẫn tác giả suốt trình tìm hiểu, nghiên cứu hoàn thành luận văn Học viên chân thành cám ơn quý quan: Trung tâm Dự báo Khí tượng Thủy văn Trung ương, Đài Khí tượng Thủy văn khu vực Nam Trung Bộ, Đài Khí tượng Thủy văn tỉnh Phú Yên, Ban huy PCTT & TKCN tỉnh Phú Yên, giúp đỡ tác giả trình tìm hiểu thực tế, thu thập liệu thực hành tác nghiệp liên quan đến luận văn Cuối Học viên xin gửi lời cám ơn đến thầy cô, đồng nghiệp, bạn bè học viên lớp cao học CH2AT tạo điều kiện thuận lợi, giúp đỡ tác giả suốt trình thực luận văn Người viết Võ Anh Kiệt LỜI CAM ĐOAN Tôi Võ Anh Kiệt, xin cam đoan luận văn cơng trình nghiên cứu riêng Những nội dung luận văn trung thực chưa công bố cơng trình khoa học Người cam đoan Võ Anh Kiệt DANH MỤC CHỮ VIẾT TẮT TANK : Mơ hình thủy văn thông số tập trung Nhật Bản MIKE 11 : Mơ hình thủy lực chiều Viện Thủy lực Đan Mạch MIKE 21 : Mơ hình thủy lực chiều Viện Thủy lực Đan Mạch DHI : Danish Hydraulic Institute - Viện Thủy lực Đan Mạch WMO : World Meteorological Organizaton - Tổ chức Khí tượng Thế giới HR : Ký hiệu trạm Thủy văn RR : Ký hiệu điểm đo mưa nhân dân UCAR : University Corporaton for Atmospheric Research GFS : Global Forecast System - Hệ thống dự báo toàn cầu IFS : Integnated Forecastng system – Hệ thống dự báo tích hợp NCEP : Natonal Centers for Environmental Prediction - Trung tâm dự báo môi trường Quốc gia Mỹ MRF : Medium Range Forecast - Mơ hình dự báo hạn vừa Pbl : Planetary boundary layer - Lớp biên hành tnh IR : Infra Red - Bức xạ hồng ngoại UV : Ultraviolet - Bức xạ tử ngoại K : Độ Kelvin CAPE : Convective available potental energy - Chỉ số bất ổn định SAS : Sơ đồ đối lưu Schubert - Arakawa-Schubert UTC : Giờ quốc tế ATNĐ : Áp thấp nhiệt đới KKL Khơng khí lạnh XTNĐ : : Xoáy thuận nhiệt đới ITCZ : Intertropical Convergence Zone - Dải hội tụ nhiệt đới KTTV : Khí tượng thủy văn CSDL : Cơ sở liệu MỤC LỤC MỞ ĐẦU .1 CHƯƠNG TỔNG QUAN VỀ MƠ HÌNH DỰ BÁO LŨ VÀ LƯU VỰC NGHIÊN CỨU 1.1 KHÁI QT MƠ HÌNH DỰ BÁO 1.1.1 Quá trình phát triển mơ hình dự báo thủy văn 1.1.2 Phân loại mơ hình thủy văn .4 1.2 TÌNH HÌNH NGHIÊN CỨU TRONG VÀ NGỒI NƯỚC 1.2.1 Đánh giá tổng quan tnh hình nghiên cứu ngồi nước 1.2.2 Đánh giá tổng quan tnh hình nghiên cứu nước 10 1.2.3 Đánh giá tổng quan tnh hình nghiên cứu lưu vực 12 1.3 ĐẶC ĐIỂM ĐỊA LÝ TỰ NHIÊN 14 1.3.1 Đặc điểm hình thái lưu vực sông 14 1.3.2 Đặc điểm địa hình, địa chất, thổ nhưỡng, thảm phủ thực vật 15 1.3.3 Khái quát đặc điểm khí hậu thủy văn 20 1.4 ĐẶC ĐIỂM MƯA LŨ LƯU VỰC SÔNG KỲ LỘ 22 1.4.1 Đặc điểm chung .22 1.4.2 Các đợt mưa lũ điển hình 23 1.5 ĐIỀU KIỆN KINH TẾ Xà HỘI .25 1.6 MẠNG LƯỚI TRẠM KHÍ TƯỢNG THỦY VĂN .26 CHƯƠNG 28 CƠ SỞ MƠ HÌNH VÀ BỘ SỐ LIỆU NGHIÊN CỨU 28 2.1 CƠ SỞ LÝ THUYẾT MƠ HÌNH MARINE 28 2.1.1 Các giả thiết nguyên lý mô 28 2.1.2 Lý thuyết sinh dòng chảy 29 2.1.3 Lý thuyết thấm 30 2.2 CƠ SỞ LÝ THUYẾT MƠ HÌNH MUSKINGUM 32 2.2.1 Các phương trình 32 2.2.2 Xây dựng mơ hình Muskingum 34 2.3 TÍNH TỐN VÀ XỬ LÝ DỮ LIỆU 34 2.3.1 Số liệu địa hình 34 2.3.2 Số liệu thảm phủ 37 2.3.3 Số liệu đất 38 2.3.4 Xử lý mưa phân bố theo không gian 39 2.3.5 Số liệu sông .41 2.4 DẪN SUẤT ĐẦU VÀO MƠ HÌNH MARINE 42 2.4.1 Xây dựng đồ dẫn suất từ đồ DEM phục vụ đầu vào cho mơ hình MARINE .42 2.4.2 Kết xây dựng đồ dẫn suất đầu vào mơ hình MARINE .45 CHƯƠNG 49 ỨNG DỤNG MƠ HÌNH MARINE VÀ MUSKINGUM XÂY DỰNG PHƯƠNG ÁN CẢNH BÁO, DỰ BÁO LŨ VÀ ĐỀ XUẤT CÁC GIẢI PHÁP PHÒNG CHỐNG LŨ 49 3.1 HIỆU CHỈNH BỘ THƠNG SỐ MƠ HÌNH 49 3.2 KIỂM ĐỊNH BỘ THƠNG SỐ MƠ HÌNH 52 3.3 DỰ BÁO THỬ TRẬN LŨ NĂM 2017 54 3.4 XÂY DỰNG PHƯƠNG ÁN CẢNH BÁO, DỰ BÁO LŨ SƠNG KỲ LỘ 56 3.5 HƯỚNG DẪN QUY TRÌNH DỰ BÁO LŨ BẰNG MƠ HÌNH MARINE VÀ MUSKINGUM .56 3.6 CÁC GIẢI PHÁP PHÒNG CHỐNG LŨ TRÊN LƯU VỰC SÔNG KỲ LỘ .58 3.6.1 Các giải pháp cơng trình 58 3.6.2 Các giải pháp phi công trình 60 KẾT LUẬN VÀ KIẾN NGHỊ 63 TÀI LIỆU THAM KHẢO .65 PHẦN PHỤ LỤC .67 DANH MỤC HÌNH Hình 1.1 Cây phân loại mơ hình tốn thủy văn Hình 1.2 Bản đồ địa hình lưu vực sông Kỳ Lộ 17 Hình 1.3 Mức nước lũ ngập sâu trạm thủy văn Hà Bằng trạm thủy văn Hà Bằng, trận lũ lịch sử tháng 11/2009 25 Hình 1.4 Lũ trơi cơng trình đo mực nước tự ghi trạm thủy văn Hà Bằng, trận lũ lịch sử tháng 11/2009 25 Hình 1.5 Bản đồ sơng suối lưu vực sơng Kỳ Lộ 27 Hình 2.1 Sơ đồ tổng qt mơ hình MARINE 28 Hình 2.2 Mơ tả vùng thấm lý thuyết Green Ampt 31 Hình 2.3 Code mơ hình Muskingum 34 Hình 2.4 Chi tết địa hình khu vực trung lưu 35 Hình 2.5 Kiểm tra đồ mơ hình số độ cao (DEM) đường bình đồ địa hình 35 Hình 2.6 Kiểm tra đồ mơ hình số độ cao (DEM) hệ thống sông suối 36 Hình 2.7 Bản đồ mơ hình số độ cao (DEM) lưu vực sông Kỳ Lộ 36 Hình 2.8 Bản đồ lớp phủ thực vật lưu vực sông Kỳ Lộ 38 Hình 2.9 Bản đồ đất lưu vực sơng Kỳ Lộ 39 Hình 2.10 Bản đồ phân bố mưa lưới theo đa giác Theissen 41 Hình 2.11 Bản đồ mạng lưới sơng tính tốn liệu cho MARINE 42 Hình 2.12 Sơ đồ mơ tả q trình tạo đồ dẫn suất từ DEM 43 Hình 2.13 Sơ đồ mơ tả q trình tạo hướng chảy từ đồ DEM 44 Hình 2.14 Sơ đồ mô số ô hội tụ nước hình thành mạng sơng 44 Hình 2.15 Tạo đồ dẫn suất từ AML fles 45 Hình 2.16 Bản đồ hướng dòng chảy lưu vực sơng Kỳ Lộ .46 Hình 2.17 Bản đồ độ rỗng đất lưu vực sông Kỳ Lộ 46 Hình 2.18 Bản đồ độ dẫn thủy lực bão hòa lưu vực sơng Kỳ Lộ 46 Hình 2.19 Bản đồ lượng trữ lưu vực sông Kỳ Lộ 46 Hình 2.20 Bản đồ độ dốc lưu vực sông Kỳ Lộ 47 Hình 2.21 Bản đồ lực mao dẫn lưu vực sông Kỳ Lộ .47 Hình 2.22 Hàm lượng nước ban đầu lưu vực sông Kỳ Lộ 47 Hình 2.23 Bản đồ hệ số nhám Manning lưu vực sông Kỳ Lộ 47 Hình 2.24 Bản đồ miền tính lưu vực sơng Kỳ Lộ 48 - Mạng lưới trạm khí tượng, thủy văn đo mưa thiếu so với yêu cầu, cùng thượng nguồn, gây nhiều khó khăn công tác dự báo lũ vùng hạ lưu Vì cần bổ sung thêm mạng lưới trạm khí tượng thủy văn lưu vực sông để đảm bảo quan trắc dự báo lũ chủ động hiệu (hiện bổ sung thêm trạm, xây dựng vận hành thử nghiệm từ năm 2017) - Kiến thức bão lũ nói riêng, loại hình thiên tai khác nói chung cộng đồng xã hội nhiều hạn chế, cần phải thường xuyên phổ cập, thơng tn tun truyền nhiều hình thức để nâng cao nửa ý thức cộng đồng việc phòng tránh thiên tai ứng phó với biến đổi khí hậu toàn cầu./ TÀI LIỆU THAM KHẢO Nguyễn Hữu Khải, Nguyễn Thanh Sơn (2003), Mơ hình tốn thủy văn, Nxb Đại học Quốc gia Hà Nội IFAS System Instruction GuideBook (2011), Internatonal Centre for Water Hazard and Risk Management (ICHARM) Ven Techow, David R.Maidment, Larry W.Mays (1988), Applied Hydrology, New York : McGraw-Hill, c1988 Nguyễn Văn Tuần, Đoàn Quyết Trung, Bùi Văn Đức (2003), Dự báo thủy văn, Nxb Đại học Quốc gia Hà Nội Bùi Văn Chanh (2017), Tiểu luận tổng quan, Nghiên cứu tch hợp mơ hình để khơi phục dòng chảy dự báo thủy văn cho lưu vực thiếu số liệu quan trắc, Viện Khoa học Khí tượng Thủy văn Biến đổi khí hậu Bùi Văn Chanh (2017), Báo cáo khoa học, Nghiên cứu ứng dụng mơ hình mưa dòng chảy thơng số phân bố để xây dựng công cụ dự báo thủy văn cho lưu vực sơng Dinh Ninh Hòa, Trung tâm Khí tượng Thủy văn Quốc gia Nguyễn Văn Lý (2015), báo cáo tổng kết đề tài NCKH cấp sở, Nghiên cứu xây dựng công cụ cảnh báo, dự báo lũ sơng Kỳ Lộ, tỉnh Phú n, Trung tâm Khí tượng Thủy văn quốc gia Nguyễn Lan Châu, Đặng Thanh Mai, Trịnh Thu Phương (2005) Các toán việc ứng dụng mơ hình thủy văn Marine để mơ dự báo lũ sông Đà Hội nghị Khoa học Cơng nghệ phục vụ dự báo Khí tượng thủy văn lần thứ VI Nguyễn Lan Châu (2006) Đề tài cấp Bộ Tài nguyên Môi trường “Nghiên cứu xây dựng công nghệ dự báo lũ phục vụ điều tiết hồ Hòa Bình cơng tác phòng chống lũ lụt” 10 Viện Kỹ thuật Biển (2012) Đánh giá trạng môi trường lưu vực sông Cái (Nha Trang) sơng Dinh (Ninh Hòa), phân tch ngun nhân đề xuất giải pháp bảo vệ, cải thiện chất lượng mơi trường 11 Đỗ Thị Chính (2010) Báo cáo tổng kết đề tài “Nghiên cứu đánh giá hiểm họa ngập lụt sử dụng đất vùng hạ du lưu vực sơng Sài Gòn – Đồng Nai điều kiện nước biển dâng biến đổi khí hậu tồn cầu”, Đề tài cấp sở Viện Khoa học thủy lợi miền Nam 12 Nguyễn Quốc Hiệp (2012) Bài báo: “Hệ thống giám sát điều khiển, dự báo lũ, dự báo ngập lụt vùng hạ du điều hành hồ chứa theo thời gian thực phục vụ cơng tác phòng chống lụt bão, giảm nhẹ thiên tai” 13 Võ Anh Kiệt (2010) Báo cáo Hội thảo “Tình hình bão lũ tỉnh Phú Yên năm gần đây, xác định nguyên nhân trận lũ lụt lịch sử năm 2009 địa bàn tỉnh Phú Yên” Hội thảo khoa học BĐKH Phú Yên tháng 4.2010 14 Trần Công Danh (2012) Đề tài NCKH cấp tỉnh “nghiên cứu, bổ sung đặc điểm Khí hậu thủy văn Phú Yên xây dựng đồ nguy ngập lụt lưu vực sông Kỳ Lộ đoạn từ Phú Mỡ đến hạ lưu” 15 Võ Anh Kiệt, Nguyễn Thị Thu Loan (2014), Đặc điểm khí hậu thủy văn khu vực Nam Trung Bộ năm gần Tập san Quy hoạch Tài nguyên nước Miền Trung 16 Võ Anh Kiệt, Lại Thị Lương (2015), Đặc trưng hình thái sơng ngòi Nam Trung Bộ Tập san Quy hoạch Tài nguyên nước Miền Trung PHẦN PHỤ LỤC Biểu quan hệ Q = f(H) Trạm Hà Bằng Sông: Kỳ Lộ Q H H 300 310 320 330 340 350 360 370 380 390 400 410 420 430 440 450 460 470 480 490 500 510 520 530 540 550 560 570 580 590 600 610 620 630 640 650 660 670 Trạm: Hà Bằng Q 0.156 0.176 0.196 0.216 0.236 0.256 0.276 0.296 0.578 1.36 2.58 4.31 6.58 9.47 13.0 17.3 22.2 28.0 34.6 42.1 50.4 59.7 70.0 81.2 93.5 107 121 137 154 171 191 211 232 255 279 305 332 360 0.158 0.178 0.198 0.218 0.238 0.258 0.278 0.321 0.639 1.46 2.73 4.51 6.84 9.79 13.4 17.7 22.8 28.6 35.3 42.9 51.3 60.7 71.0 82.4 94.8 108 123 138 155 173 192 213 235 258 282 308 335 363 0.160 0.180 0.200 0.220 0.240 0.260 0.280 0.346 0.703 1.57 2.89 4.71 7.11 10.1 13.8 18.2 23.3 29.3 36.0 43.7 52.2 61.7 72.1 83.6 96.1 110 124 140 157 175 194 215 237 260 284 310 337 366 0.162 0.182 0.202 0.222 0.242 0.262 0.282 0.371 0.772 1.68 3.05 4.93 7.38 10.5 14.2 18.7 23.9 29.9 36.8 44.5 53.1 62.7 73.2 84.8 97.4 111 126 142 159 177 196 217 239 262 287 313 340 369 0.164 0.184 0.204 0.224 0.244 0.264 0.284 0.396 0.844 1.80 3.21 5.15 7.66 10.8 14.6 19.2 24.5 30.5 37.5 45.3 54.0 63.7 74.3 86.0 98.7 112 127 143 161 179 198 219 241 265 290 316 343 372 0.166 0.186 0.206 0.226 0.246 0.266 0.286 0.421 0.920 1.92 3.38 5.37 7.95 11.2 15.0 19.6 25.0 31.2 38.2 46.1 54.9 64.7 75.5 87.2 100 114 129 145 162 181 201 221 244 267 292 318 346 375 0.168 0.188 0.208 0.228 0.248 0.268 0.288 0.446 1.00 2.04 3.55 5.60 8.24 11.5 15.5 20.2 25.6 31.9 39.0 47.0 55.9 65.7 76.6 88.5 101 115 130 147 164 183 203 224 246 270 295 321 349 378 0.170 0.190 0.210 0.230 0.250 0.270 0.290 0.471 1.08 2.17 3.73 5.84 8.54 11.9 15.9 20.7 26.2 32.5 39.7 47.8 56.8 66.8 77.7 89.7 103 117 132 148 166 185 205 226 248 272 297 324 351 381 0.172 0.192 0.212 0.232 0.252 0.272 0.292 0.496 1.17 2.30 3.92 6.08 8.84 12.2 16.3 21.2 26.8 33.2 40.5 48.7 57.8 67.8 78.9 91.0 104 118 134 150 168 187 207 228 251 275 300 326 354 384 0.174 0.194 0.214 0.234 0.254 0.274 0.294 0.521 1.27 2.44 4.11 6.33 9.15 12.6 16.8 21.7 27.4 33.9 41.3 49.5 58.7 68.9 80.0 92.2 105 120 135 152 170 189 209 230 253 277 302 329 357 387 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.282 0.782 1.22 1.73 2.27 2.89 3.53 4.30 4.90 5.80 6.60 7.50 8.30 9.30 10.3 11.2 12.3 13.5 14.0 16.0 17.0 17.0 20.0 20.0 21.0 23.0 24.0 26.0 27.0 28.0 30.0 Q H H 680 690 700 710 720 730 740 750 760 770 780 790 800 810 820 830 840 850 860 870 880 890 900 910 920 930 940 950 960 970 980 990 1000 1010 1020 1030 1040 1050 1060 1070 1080 1090 1100 1110 Q 390 421 453 487 523 560 599 639 681 724 769 816 864 914 966 1019 1074 1131 1190 1251 1313 1377 1444 1512 1582 1653 1727 1803 1881 1960 2042 2126 2212 2299 2389 2481 2576 2672 2770 2871 2974 3079 3186 3295 393 424 457 491 527 564 603 643 685 728 774 820 869 919 971 1025 1080 1137 1196 1257 1320 1384 1450 1519 1589 1661 1735 1811 1889 1968 2050 2134 2220 2308 2399 2491 2585 2682 2780 2881 2984 3089 3197 3306 396 427 460 494 530 568 607 647 689 733 778 825 874 924 976 1030 1086 1143 1202 1263 1326 1391 1457 1525 1596 1668 1742 1818 1896 1977 2059 2143 2229 2317 2408 2500 2595 2691 2790 2891 2994 3100 3207 3317 399 430 463 498 534 571 611 651 693 737 783 830 879 929 982 1036 1091 1149 1208 1269 1332 1397 1464 1532 1603 1675 1750 1826 1904 1985 2067 2151 2238 2326 2417 2509 2604 2701 2800 2901 3005 3111 3218 3328 402 434 467 502 538 575 615 655 698 742 787 835 884 934 987 1041 1097 1155 1214 1276 1339 1404 1471 1539 1610 1683 1757 1834 1912 1993 2075 2160 2247 2335 2426 2519 2614 2711 2810 2912 3015 3121 3229 3340 405 437 470 505 541 579 619 660 702 746 792 840 889 940 992 1047 1103 1161 1220 1282 1345 1410 1477 1546 1617 1690 1765 1842 1920 2001 2084 2168 2255 2344 2435 2528 2623 2721 2820 2922 3026 3132 3240 3351 408 440 474 509 545 583 623 664 706 751 797 844 894 945 998 1052 1108 1166 1226 1288 1352 1417 1484 1553 1624 1697 1772 1849 1928 2009 2092 2177 2264 2353 2444 2538 2633 2731 2830 2932 3036 3143 3251 3362 411 444 477 512 549 587 627 668 711 755 801 849 899 950 1003 1058 1114 1172 1232 1294 1358 1424 1491 1560 1632 1705 1780 1857 1936 2017 2100 2186 2273 2362 2454 2547 2643 2741 2840 2943 3047 3153 3262 3373 415 447 481 516 553 591 631 672 715 760 806 854 904 955 1008 1063 1120 1178 1238 1301 1364 1430 1498 1567 1639 1712 1788 1865 1944 2026 2109 2194 2282 2371 2463 2557 2652 2750 2851 2953 3057 3164 3273 3384 418 450 484 519 556 595 635 676 720 764 811 859 909 960 1014 1069 1126 1184 1245 1307 1371 1437 1505 1574 1646 1720 1795 1873 1952 2034 2117 2203 2291 2380 2472 2566 2662 2760 2861 2963 3068 3175 3284 3396 31.0 32.0 34.0 36.0 37.0 39.0 40.0 42.0 43.0 45.0 47.0 48.0 50.0 52.0 53.0 55.0 57.0 59.0 61.0 62.0 64.0 67.0 68.0 70.0 71.0 74.0 76.0 78.0 79.0 82.0 84.0 86.0 87.0 90.0 92.0 95.0 96.0 98.0 101 103 105 107 109 112 68 Q H H 1120 1130 1140 1150 1160 1170 1180 1190 1200 1210 1220 1230 1240 1250 1260 1270 1280 1290 1300 1310 1320 1330 1340 1350 Q 3407 3521 3637 3755 3876 3999 4125 4252 4382 4515 4650 4787 4927 5069 5214 5361 5510 5662 5817 5974 6134 6296 6461 6628 3418 3532 3649 3767 3888 4012 4137 4265 4396 4528 4663 4801 4941 5083 5228 5376 5526 5678 5833 5990 6150 6312 6477 6645 3429 3544 3660 3779 3900 4024 4150 4278 4409 4542 4677 4815 4955 5098 5243 5391 5541 5693 5848 6006 6166 6329 6494 6662 3441 3555 3672 3791 3913 4037 4163 4291 4422 4555 4691 4829 4969 5112 5258 5405 5556 5709 5864 6022 6182 6345 6511 6679 3452 3567 3684 3803 3925 4049 4175 4304 4435 4569 4704 4843 4983 5127 5272 5420 5571 5724 5880 6038 6198 6362 6528 6696 3463 3578 3696 3815 3937 4062 4188 4317 4448 4582 4718 4857 4998 5141 5287 5435 5586 5739 5895 6054 6215 6378 6544 6713 3475 3590 3708 3827 3950 4074 4201 4330 4462 4596 4732 4871 5012 5156 5302 5450 5601 5755 5911 6070 6231 6395 6561 6730 3486 3602 3720 3840 3962 4087 4214 4343 4475 4609 4746 4885 5026 5170 5316 5465 5617 5770 5927 6086 6247 6411 6578 6747 3498 3613 3731 3852 3974 4099 4227 4356 4488 4623 4760 4899 5040 5185 5331 5480 5632 5786 5943 6102 6263 6428 6595 6764 3509 3625 3743 3864 3987 4112 4239 4369 4502 4636 4773 4913 5055 5199 5346 5495 5647 5802 5958 6118 6280 6444 6611 6781 114 116 118 121 123 126 127 130 133 135 137 140 142 145 147 149 152 155 157 160 162 165 167 Kết tối ưu thông số * Trận mưa lũ ngày - 12/11/2010 (hiệu chỉnh) STT 10 11 12 13 14 15 16 17 18 H đo 599 598 598 598 598 602 605 609 618 626 635 646 657 668 679 689 700 707 Qd đo 189 187 187 187 187 194 201 209 228 246 267 295 324 354 387 418 453 477 Qt tính 276 272 269 266 263 264 266 268 278 287 298 314 330 347 366 384 404 417 (Qdi-�̅� )2 152362 153927 153927 153927 153927 148484 143138 137149 123437 111113 97554 80847 65196 50776 36993 26029 15961 10473 (Qd-Qt)2 7569 7168 6669 6188 5725 4947 4225 3520 2500 1681 961 348 36 49 441 1179 2401 3600 STT H đo 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 770 762 755 748 741 734 728 722 717 712 707 702 698 693 690 686 683 679 69 Qd Qt đo tính 724 530 689 506 660 486 631 466 603 446 575 426 553 410 530 394 512 380 494 367 477 354 460 342 447 332 430 319 421 312 408 302 399 294 387 285 (Qdi-̅ )2 20928 12026 6507 2669 560 19 694 2434 4534 7282 10473 14241 17513 22301 25070 29356 32521 36993 (Qd-Qt) 37507 33489 30276 27335 24649 22102 20354 18587 17336 16129 15047 14003 13302 12299 11946 11307 10955 10404 STT 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 H đo 714 721 720 719 718 712 708 703 699 696 692 694 695 697 697 698 698 716 733 751 794 836 875 904 930 949 964 976 981 984 984 980 972 959 945 929 913 901 886 871 860 851 844 839 832 825 819 810 802 795 787 779 Qd đo 502 527 523 519 516 494 481 463 450 440 427 434 437 444 444 447 447 509 571 643 835 1052 1282 1471 1653 1795 1912 2009 2050 2075 2075 2042 1977 1873 1765 1646 1532 1450 1352 1257 1190 1137 1097 1069 1030 992 960 914 874 840 801 764 Qt tính 431 444 439 433 428 410 399 384 372 362 351 352 351 353 324 307 301 340 381 429 557 702 869 1058 1298 1503 1632 1738 1860 1962 1961 1896 1804 1686 1554 1477 1470 1412 1294 1171 1073 991 925 873 818 770 730 683 642 612 583 558 (Qdi-̅ )2 (Qd-Qt) STT 5981 2739 3174 3640 4011 7282 9670 13534 16728 19414 23206 21122 20259 18316 18316 17513 17513 4947 69 4053 65364 223412 493737 795065 1152755 1477840 1775994 2043940 2162853 2237012 2237012 2139387 1953465 1673567 1405800 1137773 907569 758056 597010 459229 372911 310989 267976 239771 203098 170292 144905 112000 86827 67946 49135 34101 5088 6834 7112 7396 7744 7000 6779 6294 6084 6032 5827 6669 7339 8281 14296 19497 21180 28459 36125 45906 77458 122722 170459 170569 126025 85459 78400 73260 36227 12694 12996 21413 29814 35094 44521 28561 3885 1469 3403 7339 13689 21219 29584 38547 44803 49432 53053 53515 53669 52136 47379 42573 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 70 Qd Qt ̅ (Qd-Qt) đo tính (Qdi- ) 676 378 278 40536 10060 673 369 265 44241 10913 670 360 248 48108 12604 667 351 240 52137 12417 663 340 232 57282 11578 659 329 226 62668 10650 656 321 221 66737 9967 653 313 216 70935 9319 651 308 214 73623 8911 648 300 209 78028 8348 646 295 206 80847 8004 643 287 200 85460 7511 641 282 197 88409 7219 639 277 194 91407 6933 636 270 189 95689 6545 634 265 186 98807 6273 631 258 181 103257 5903 629 253 178 106495 5645 626 246 173 111113 5295 624 241 170 114471 5050 623 239 169 115828 4951 621 235 166 118567 4761 619 230 163 122035 4529 617 226 160 124846 4347 615 221 157 128404 4126 613 217 154 131287 3952 612 215 153 132741 3865 610 211 150 135671 3697 608 207 148 138634 3532 607 205 146 140127 3450 605 201 144 143138 3291 604 198 142 145417 3173 602 194 139 148484 3021 601 192 138 150029 2945 600 191 137 150805 2905 598 187 134 153927 2760 597 185 133 155501 2687 596 183 132 157082 2615 595 181 131 158671 2544 594 179 129 160269 2473 593 177 128 161874 2404 593 177 128 161874 2401 592 175 127 163487 2333 592 175 127 163487 2330 591 173 125 165109 2263 591 173 125 165109 2259 591 173 126 165109 2256 591 173 126 165109 2253 591 173 126 165109 2250 591 173 126 165109 2247 591 173 126 165109 2244 591 173 126 165109 2240 T.bình 579 ∑=260021 ∑=18384 H đo STT H đo Qd đo Qt tính (Qdi-̅ )2 (Qd-Qt) STT H đo Qd đo Qt tính (Qdi-̅ )2 (Qd-Qt) Qd Qt đo tính 742 774 681 687 631 614 579 546 487 456 467 420 450 390 430 360 430 346 405 322 381 303 357 286 335 271 313 255 292 241 262 220 251 211 239 202 230 195 221 188 213 181 203 173 203 172 192 163 185 157 175 149 170 144 164 139 159 134 155 130 150 126 147 117 142 102 138 98 134 95 130 93 126 91 123 90 120 88 117 87 114 85 111 83 108 81 105 79 103 78 101 76 (Qdi-̅ )2 (Qd-Qt) 143807 101263 71941 46751 15430 10862 7607 4518 4518 1782 332 33 772 2478 5010 10157 12495 15322 17631 20102 22434 25530 25530 29166 31606 35262 37165 39514 41527 43173 45276 46562 48744 50527 52341 54187 56065 57495 58943 60408 61892 63394 64913 66451 67486 68529 1045 40 278 1089 940 2240 3640 4900 7000 6944 6032 4994 4096 3325 2635 1792 1573 1344 1225 1111 1024 900 982 841 784 671 657 625 608 608 590 928 1568 1635 1516 1359 1213 1102 1007 922 849 784 727 673 637 602 Nash = 92.93 * Trận mưa lũ ngày 19 - 22/10/2011 (kiểm định) STT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 H đo 446 446 447 449 450 452 453 453 459 464 470 475 481 481 488 495 502 509 516 516 523 531 538 557 575 594 619 645 670 705 740 775 818 842 864 883 891 894 891 886 878 868 859 848 835 818 Qd đo 15 15 16 17 17 18 19 19 22 24 28 31 35 35 41 46 52 59 66 66 73 82 91 117 145 179 230 292 360 470 599 746 955 1086 1214 1332 1384 1404 1384 1352 1301 1238 1184 1120 1047 955 Qt tính 130 127 125 122 119 117 115 112 111 109 109 108 107 104 105 106 107 108 110 107 109 112 89 87 100 120 154 195 240 314 414 575 833 1030 1167 1287 1416 1515 1500 1436 1354 1262 1167 1126 1146 1082 (Qdi-̅ )2 (Qd-Qt)2 120952 120952 120257 119565 119565 118874 118186 118186 116132 114773 112078 110079 107441 107441 103543 100350 96585 92283 88079 88079 83973 78838 73865 60408 47429 33776 17631 5010 11496 55799 146857 350723 523045 724573 939385 1042888 1084136 1042888 978554 880254 766008 674400 573380 468155 350723 13225 12544 11808 11095 10472 9801 9152 8587 7862 7225 6507 5878 5232 4807 4139 3560 2988 2434 1936 1681 1272 880 878 2013 3445 5837 9459 14395 24435 34299 29355 14965 3136 2240 2025 1003 12321 13533 7000 2774 592 300 40 9867 16044 STT 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 71 H đo 774 760 748 735 710 704 699 693 693 685 677 669 661 653 645 633 628 623 619 615 611 606 606 601 597 592 589 586 583 581 578 576 573 571 568 566 563 561 559 557 555 553 551 549 547 546 STT 47 48 H đo 802 788 Qd đo 874 806 Qt tính 975 871 (Qdi-̅ )2 (Qd-Qt) 261345 196443 10201 4182 STT 95 96 T.bình Nash = 97.12 H đo 544 542 Qd Qt ̅ đo tính (Qdi- ) 99 75 69581 96 73 71172 363 ∑=162011 (Qd-Qt) 568 520 ∑=4673 Kết dự báo thử trận lũ từ ngày – 8/11/2017 STT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 Thời gian 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/4/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 Thực đo 592 598 606 613 636 653 709 783 858 932 1006 1068 1095 1117 1132 1138 1135 1125 1115 1098 1079 1059 1040 1015 989 967 945 923 901 886 868 855 837 822 813 794 783 777 767 Dự báo 675 677 681 684 694 703 733 776 825 875 926 1014 1015 1093 1160 1187 1161 1104 1053 1018 1020 1000 985 972 955 944 927 916 894 870 841 815 789 773 762 753 740 732 724 Sai số 83 79 76 71 58 50 24 33 57 80 54 80 24 28 49 26 21 62 80 59 59 55 43 34 23 18 7 16 27 40 48 49 51 41 43 45 43 Đánh STT giá sai 60 sai 61 sai 62 sai 63 sai 64 sai 65 66 67 68 sai 69 sai 70 sai 71 sai 72 73 74 sai 75 76 77 sai 78 sai 79 sai 80 sai 81 sai 82 83 84 85 86 87 88 89 90 91 sai 92 sai 93 sai 94 95 96 sai 97 98 Thời Thực Dự gian đo báo 11/6/2017 743 735 11/6/2017 740 715 11/6/2017 735 702 11/6/2017 730 701 11/6/2017 723 696 11/6/2017 715 693 11/6/2017 709 696 11/6/2017 703 705 11/6/2017 700 687 11/6/2017 697 675 11/6/2017 697 673 11/6/2017 697 672 11/7/2017 697 672 11/7/2017 697 673 11/7/2017 695 671 11/7/2017 692 670 11/7/2017 689 670 11/7/2017 685 668 11/7/2017 680 664 11/7/2017 675 661 11/7/2017 671 657 11/7/2017 666 655 11/7/2017 662 651 11/7/2017 657 648 11/7/2017 654 646 11/7/2017 651 644 11/7/2017 648 642 11/7/2017 645 640 11/7/2017 643 638 11/7/2017 640 636 11/7/2017 638 635 11/7/2017 635 633 11/7/2017 632 631 11/7/2017 629 629 11/7/2017 627 627 11/7/2017 624 626 11/8/2017 623 625 11/8/2017 622 624 11/8/2017 620 623 Sai số 25 33 29 27 22 13 13 22 24 25 25 24 24 22 19 17 16 14 14 11 11 5 1 2 Đánh giá đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/5/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 11/6/2017 753 746 735 727 720 715 708 702 702 702 707 711 717 722 726 730 736 742 744 745 713 708 700 693 689 685 680 676 679 690 692 690 696 698 699 705 709 715 723 740 40 99 11/8/2017 38 100 11/8/2017 35 101 11/8/2017 34 102 11/8/2017 31 103 11/8/2017 30 104 11/8/2017 28 105 11/8/2017 26 106 11/8/2017 23 107 11/8/2017 12 108 11/8/2017 15 109 11/8/2017 21 110 11/8/2017 21 111 11/8/2017 24 112 11/8/2017 27 113 11/8/2017 25 114 11/8/2017 27 115 11/8/2017 27 116 11/8/2017 21 117 11/8/2017 118 11/8/2017 Tống số lần dự báo: 118 Số lần dự báo đúng: 98 Chất lượng dự báo: 83,1% 618 616 614 612 610 608 607 605 604 602 601 600 598 597 596 595 594 593 592 591 622 621 619 617 616 615 614 613 612 611 610 609 608 607 606 605 604 603 602 601 5 7 8 9 10 10 10 10 10 10 10 10 đúng đúng đúng đúng đúng đúng đúng đúng đúng đúng Tính sai số dự báo lũ (Phương pháp tnh theo Thông tư 42) Tháng 11 Ngày 12 13 14 15 16 17 Năm 2003 Giờ H 439 469 13 534 19 800 905 977 13 955 19 901 919 916 13 820 19 734 685 657 13 635 19 616 601 589 13 579 19 568 560 ΔH 30 65 266 105 72 -22 -54 18 -3 -96 -86 -49 -28 -22 -19 -15 -12 -10 -11 -8 -9 (∆� − ∆̅�̅)� 663 3689 68507 10148 4589 690 3394 189 53 10052 8147 2837 1041 690 541 371 264 203 233 150 176 Tháng Ngày 10 11 12 13 14 15 Năm 2004 Giờ H ΔH 405 23 428 13 428 -1 19 427 436 36 472 -10 13 462 -16 19 446 -7 439 -2 437 13 438 19 442 14 456 319 775 -101 13 674 -93 19 581 -44 537 -25 512 -20 13 492 -11 19 481 -17 464 -8 (∆� − ∆̅�̅)� 458 55 1183 135 310 74 13 154 100737 10529 8951 2080 708 467 159 346 92 13 19 Tháng 10 Ngày 19 20 21 22 23 24 25 26 27 28 551 545 537 TB -6 -8 105 150 DLC SS 5082 71 48 Năm 2005 Giờ H ΔH 450 -4 446 13 449 19 454 457 18 475 295 13 770 -6 19 764 -102 662 -55 607 -33 13 574 -17 19 557 564 62 626 83 13 709 -7 19 702 68 770 95 865 103 13 968 -17 19 951 -107 844 -78 766 -4 13 762 18 19 780 -10 770 10 780 15 13 795 37 19 832 -15 817 -63 754 -51 13 703 -34 19 669 -12 657 -1 656 -15 13 641 -19 19 622 -18 604 -12 592 -11 13 581 -11 19 570 13 19 (∆� − ∆̅�̅)� 50 223 85219 82 11041 3373 1302 403 15 3472 6388 102 4215 8450 9985 403 12117 6573 50 223 171 48 142 1151 327 4366 2924 1375 227 17 327 487 444 227 198 198 Tháng 12 Ngày 74 456 448 442 TB -8 -6 92 58 DLC SS 5505 74 50 Năm 2006 (∆� − ∆̅�̅)� Giờ H ΔH 408 1 409 13 409 409 19 409 409 1 13 410 47 19 413 422 68 4336 490 74 5163 25 522 13 564 10 62 19 589 599 -8 103 591 -10 148 -25 737 13 581 -15 294 19 556 541 -9 124 532 -12 200 -11 173 13 520 -10 148 19 509 499 -9 124 490 -6 66 -7 84 13 484 -6 66 19 477 471 -4 38 467 -2 17 1 13 465 19 466 TB 2.148148 462 DLC 21 SS 14 TB Tháng 11 Ngày DLC SS Năm 2007 Giờ H ΔH 542 18 560 44 13 604 80 19 684 -31 653 73 726 97 13 823 118 19 941 79 1020 -47 973 31 13 1004 87 19 1091 1 1092 -82 1010 -132 13 878 -92 19 786 -55 731 -31 700 -23 13 677 -18 19 659 -20 639 -9 630 -11 13 619 -9 19 610 -9 601 -7 594 -5 13 589 -6 19 583 TB DLC SS 4265 65 44 (∆� − ∆̅�̅)� Tháng Ngày 11 17 272 1805 6159 1057 5110 9117 13568 6003 2354 869 7307 6975 17827 8746 3194 1057 601 381 463 111 157 111 111 73 42 57 18 19 20 21 22 23 3464 59 40 24 25 26 27 28 75 Năm 2008 Giờ H ΔH 591 -15 576 -2 13 574 120 19 694 156 850 -12 838 -80 13 758 -63 19 695 -34 661 -14 647 120 13 767 121 19 888 -43 845 -3 842 -11 13 831 -34 19 797 -22 775 -25 750 -24 13 726 -30 19 696 699 42 741 13 750 43 19 793 -6 787 -30 757 -49 13 708 -23 19 685 -17 668 33 701 13 707 31 19 738 56 794 55 849 111 13 960 54 19 1014 -83 931 -78 853 13 861 19 862 12 874 -74 800 -41 13 759 -4 19 755 -21 734 -27 707 -24 13 683 -20 19 663 TB DLC (∆� − ∆̅�̅)� 273 14400 24336 144 6400 3969 1156 196 14400 14641 1849 121 1156 484 625 576 900 1764 81 1849 36 900 2401 529 289 1089 36 961 3136 3025 12321 2916 6889 6084 64 144 5476 1681 16 441 729 576 400 2968 54 SS Tháng 11 Ngày Năm 2009 Giờ H 438 438 13 437 19 436 435 434 13 445 19 765 1261 1332 13 1227 19 1105 961 845 13 768 19 739 725 727 13 719 19 706 690 677 13 668 19 659 TB Năm 2011 Ngày 10/18/2011 10/18/2011 10/18/2011 10/18/2011 10/19/2011 10/19/2011 10/19/2011 10/20/2011 10/20/2011 10/20/2011 10/20/2011 10/21/2011 10/21/2011 10/21/2011 ΔH -1 -1 -1 -1 11 320 496 71 -105 -122 -144 -116 -77 -29 -14 -8 -13 -16 -13 -9 -9 10 DLC SS (∆� − ∆̅�̅)� Tháng 11 92 113 113 113 113 96343 236577 3769 13135 17321 23596 15778 7501 1491 557 58 310 511 656 511 346 346 Ngày 10 11 18233 135 91 Giờ 13 19 13 19 13 19 13 19 13 19 13 19 37 Năm 2010 H ΔH 652 -18 634 -15 619 -12 607 -9 598 11 609 59 668 53 721 -18 703 -6 697 54 751 198 949 31 980 -79 901 -62 839 -44 795 -47 748 -36 712 -26 686 -19 667 -19 648 -14 634 -13 621 -11 610 TB -2 DLC SS Giờ H ΔH (∆� − ∆̅�̅)� 13 19 13 13 19 13 453 451 448 446 446 453 481 594 775 894 848 760 693 645 -2 -3 -2 28 113 181 119 -46 -88 -67 -48 -30 13 22 13 29 694 12398 32165 13770 2271 8038 4713 2466 1002 76 (∆� − ∆̅�̅)� 262 174 104 51 165 3700 3006 262 17 3117 39930 1078 5956 3621 1779 2041 1168 584 295 295 148 125 84 2955 54 37 10/21/2011 10/22/2011 10/22/2011 10/22/2011 10/22/2011 10/23/2011 10/23/2011 10/23/2011 10/23/2011 10/24/2011 10/24/2011 10/24/2011 10/24/2011 19 13 19 13 19 13 19 615 592 576 561 549 539 530 523 517 511 506 501 496 TB 77 -23 -16 -15 -12 -10 -9 -7 -6 -6 -5 -5 -5 608 312 277 186 136 114 75 59 59 44 44 44 DLC SS 3060 55 37 ... luận văn: Nghiên cứu ứng dụng mô hình Marine dự báo lũ sơng Kỳ Lộ tỉnh Phú Yên + Tóm tắt luận văn: Kỳ Lộ sông lớn thứ hai tỉnh Phú Yên, trải rộng hai huyện Đồng Xuân Tuy An Diễn biến lũ ngập lụt... mơ hình dự báo lũ lưu vực nghiên cứu Chương Cơ sở mơ hình số liệu nghiên cứu Chương Ứng dụng mơ hình MARINE MUSKINGUM xây dựng phương án cảnh báo, dự báo lũ giải pháp phòng chống lũ lưu vực sông. ..BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI LUẬN VĂN THẠC SĨ NGHIÊN CỨU ỨNG DỤNG MƠ HÌNH MARINE DỰ BÁO LŨ SÔNG KỲ LỘ TỈNH PHÚ YÊN VÕ ANH KIỆT CHUYÊN NGÀNH: THỦY