thử nghiệm dự báo cường độ bão bằng phương pháp lọc kalman tổ hợp đa vật lý trong điều kiện có không khí lạnh

MỤC LỤC DANH MỤC HÌNH .I DANH MỤC BẢNG III MỞ ĐẦU……………………………………………………………………………….4 CHƯƠNG 1: PHƯƠNG PHÁP LỌC KALMAN TỔ HỢP ĐA VẬT LÝ VÀ MƠ HÌNH WRF 1.1 CÁC NGHIÊN CỨU TRONG NƯỚC VÀ NGOÀI NƯỚC 1.1.1 Các nghiên cứu nước 1.1.2 Các nghiên cứu nước 1.2 PHƯƠNG PHÁP LỌC KALMAN TỔ HỢP ĐA VẬT LÝ 1.3 SƠ LƯỢC VỀ MƠ HÌNH WRF 1.3.1 Giới thiệu mơ hình WRF 1.3.2 Dữ liệu đầu vào, đầu bước chạy mơ hình WRF 10 CHƯƠNG 2: HỒN LƯU QUY MƠ LỚN TRONG ĐIỀU KIỆN CĨ KHƠNG KHÍ LẠNH SỰ TƯƠNG TÁC GIỮA KHƠNG KHÍ LẠNH VÀ BÃO 12 2.1 HỒN LƯU QUY MƠ LỚN TRONG THỜI KÌ CĨ KHƠNG KHÍ LẠNH 12 2.1.1 Áp cao Siberia .12 2.1.2 Áp cao Thái Bình Dương 13 2.2 SỰ TƯƠNG TÁC GIỮA KHƠNG KHÍ LẠNH VÀ BÃO 13 2.2.1 Dải hội tụ nhiệt đới (ITCZ): 13 2.2.2 Tín phong rãnh thấp xích đạo 14 2.2.3 Áp cao Thái Bình Dương 14 2.2.4 Áp cao lục địa châu Á 15 2.3 CƠN BÃO ĐƯỢC CHỌN ĐỂ THỬ NGHIỆM 15 2.3.1 Bão Nari (2013) 15 2.3.2 Đường bão 16 2.3.3 Ảnh vệ tinh ngày .18 2.4 THIẾT KẾ THÍ NGHIỆM 18 2.4.1 Cấu trúc miền lưới .18 2.4.2 Cấu trúc tổ hợp 19 2.4.3 Các trường hợp thử nghiệm .21 2.4.4 Các phương pháp đánh giá 21 CHƯƠNG 3: KẾT QUẢ THỬ NGHIỆM 23 3.1 KẾT QUẢ MƠ PHỎNG HỒN LƯU 23 3.2 KẾT QUẢ DỰ BÁO QUỸ ĐẠO BÃO 30 3.3 KẾT QUẢ DỰ BÁO CƯỜNG ĐỘ BÃO 30 3.3.1 Dự báo áp suất mực biển cực tiểu tâm 31 3.3.2 Dự báo tốc độ gió cực đại gần tâm .37 KẾT LUẬN 45 KHUYẾN NGHỊ 46 TÀI LIỆU THAM KHẢO 47 PHỤ LỤC .PL.1 PHỤ LỤC PL.11 DANH MỤC HÌNH Hình 1.1: Sơ đồ hệ thống mơ hình WRF 10 Hình 2.1: Đường bão Nari .17 Hình 2.2: Ảnh vệ tinh ngày 18 Hình 2.3 Cấu trúc miền tính .19 Hình 3.1: Trường đường dòng mực 850 hPa với tốc độ gió mực tương ứng bão Nari dự báo lúc 00 UTC ngày 12 tháng 10 năm 2013 ba trường hợp MI (a), MP (b), PF (c) 23 Hình 3.2: Trường đường dòng mực 850 hPa với tốc độ gió mực tương ứng bão Nari dự báo lúc 12 UTC ngày 12 tháng 10 năm 2013 ba trường hợp MI (a), MP (b), PF (c) .24 Hình 3.3: Trường đường dòng mực 850 hPa với tốc độ gió mực tương ứng bão Nari dự báo lúc 00 UTC ngày 13 tháng 10 năm 2013 ba trường hợp MI (a), MP (b), PF (c) 25 Hình 3.4: Trường đường dòng mực 850 hPa với tốc độ gió mực tương ứng bão Nari dự báo lúc 12 UTC ngày 13 tháng 10 năm 2013 ba trường hợp MI (a), MP (b), PF (c) .25 Hình 3.5: Trường đường dòng mực 850 hPa với tốc độ gió mực tương ứng bão Nari dự báo lúc 00 UTC ngày 14 tháng 10 năm 2013 ba trường hợp MI (a), MP (b), PF (c) .26 Hình 3.6: Trường đường dòng mực 850 hPa với tốc độ gió mực tương ứng bão Nari dự báo lúc 12 UTC ngày 14 tháng 10 năm 2013 ba trường hợp MI (a), MP (b), PF (c) .27 Hình 3.7: Trường đường dòng mực 850 hPa với tốc độ gió mực tương ứng bão Nari dự báo lúc 00 UTC ngày 15 tháng 10 năm 2013 ba trường hợp MI (a), MP (b), PF (c) .28 Hình 3.8: Quỹ đạo quan trắc (màu đen), thành phần tổ hợp (đường mảnh), thử nghiệm PF (hình a); MI (hình b) MP (hình c) bắt đầu dự báo lúc 00 UTC ngày 12 tháng 10 năm 2013 29 Hình 3.9: Biểu đồ khí áp mực biển cực tiểu tâm (UTC) 31 I Hình 3.10: Biến trình PMIN quan trắc PMIN trung bình thử nghiệm MI, MP, PF 33 Hình 3.11: Biến trình áp suất mực biển cực tiểu thành phần tổ hợp (đường màu tím mảnh); trung bình tổ hợp (màu đen); quan trắc (xanh dương) trường hợp MI (hình a); MP (hình b); PF (hình c); bão Nari dự báo lúc 00 UTC ngày 12/10/2013 34 Hình 3.12: Biến trình tốc độ gió cực đại gần tâm quan trắc 38 Hình 3.13: Biến trình VMAX quan trắc VMAX trung bình thử nghiệm MI, MP, PF 40 Hình 3.14: Biến trình tốc độ gió bề mặt cực đại gần tâm thành phần tổ hợp (đường màu tím mảnh); trung bình tổ hợp (màu đen); quan trắc (xanh dương) thử nghiệm MI (hình a); MP (hình b); PF (hình c); bão Nari dự báo lúc 00 UTC ngày 12/10/2013 42 II DANH MỤC BẢNG Bảng 2.1 Thông tin bão Nari .16 Bảng 2.2 Sơ đồ tham số hóa vật lý mơ hình WRF ứng với option cụ thể .19 Bảng 2.3 Thành phần tổ hợp tương ứng với lựa chọn sơ đồ vật lý khác .20 Bảng 2.4 Danh sách trường hợp thử nghiệm 21 Bảng 3.1 Sai số qũy đạo bão Nari với thời điểm dự báo lúc 00 UTC ngày 12/10/2013 (km) 30 Bảng 3.2 Kết dự báo áp suất mực biển cực tiểu tâm quan trắc trung bình thử nghiệm MI, MP, PF .32 Bảng 3.3 Sai số áp suất mực biển cực tiểu tâm bão Nari dự báo lúc 00 UTC ngày 12/10/2013 (km) 34 Bảng 3.4 Kết dự báo áp suất mực biển cực tiểu tâm quan trắc trung bình thử nghiệm MI, MP, PF 39 Bảng 3.5 Sai số tốc độ gió cực đại gần tâm bão Nari dự báo lúc 00 UTC ngày 12/10/2013 (m/s) 41 III MỞ ĐẦU Đặt vấn đề Việt Nam quốc gia ven biển nằm bên bờ Tây Biển Đông, với đường bờ biển trải dài từ Bắc xuống Nam Tuy thiên nhiên ưu đãi mặt khác, Việt Nam lại nằm khu vực ổ bão Tây Bắc Thái Bình Dương, năm chịu ảnh hưởng bão với tần suất lớn Mùa bão năm thường rơi vào mùa hè mùa thu Tuy nhiên, có bão xuất vào mùa đơng, điều kiện có khơng khí lạnh Khơng khí lạnh ảnh hưởng đến cường độ quỹ đạo bão, làm thay đổi cường độ bão khiến quỹ đạo bão trở nên khó dự báo Chính vậy, trường hợp bão kết hợp với khơng khí lạnh gây thiệt hại lớn người tài sản bão đổ vào đất liền Vì thường xuyên chịu ảnh hưởng bão nên việc dự báo bão trọng Việt Nam Theo nghiên cứu NOAA (National Oceanic and Atmospheric Administration Commissioned Corps) năm 2017, công tác dự báo bão khoảng 10 năm trở lại đây, sai số dự báo quỹ đạo bão mơ hình số giảm khoảng 50% nhờ cải tiến mơ hình, thám sát xung quanh tăng cường liên tục kỹ người dự báo cải thiện Tuy nhiên, việc dự báo cường độ bão cải thiện khoảng 10 đến 20% nhờ mơ hình động lực thống kê Vì việc nâng cao độ xác việc dự báo cường độ bão yêu cầu cấp thiết Mặt khác, năm gần tốn cải thiện trường ban đầu mơ hình phát triển mạnh Một phương pháp tiên tiến phương pháp lọc Kalman tổ hợp đa vật lí Xuất phát từ lí trên, đề tài đồ án chọn: “Thử nghiệm dự báo cường độ bão phương pháp lọc Kalman tổ hợp đa vật lý điều kiện có khơng khí lạnh” Mục tiêu Mục tiêu đồ án khảo sát xem phương pháp lọc Kalman tổ hợp đa vật lý có cải thiện khả dự báo cường độ bão điều kiện có khơng khí lạnh so với phương pháp khác Nhiệm vụ Để đạt mục tiêu đồ án thực nhiệm vụ: Tìm hiểu phương pháp lọc Kalman tổ hợp đa vật lý mơ hình WRF, hồn lưu qui mơ lớn trường hợp có bão hoạt động điều kiện khơng khí lạnh, chọn bão để thử nghiệm, áp dụng lọc Kalman tổ hợp đa vật lý để dự báo cường độ bão điều kiện có khơng khí lạnh cuối đánh giá kết sai số dự báo cường độ bão so với quan trắc, so sánh sai số dự báo cường độ bão với kết dự báo phương pháp khác Nội dung nghiên cứu Nội dung nghiên cứu bao gồm nghiên cứu hồn lưu qui mơ lớn thời kì có bão khơng khí lạnh hoạt động, nghiên cứu phương pháp lọc Kalman tổ hợp đa vật lí mơ hình WRF, chọn bão nằm phạm vi nghiên cứu để thử nghiệm, áp dụng phương pháp lọc Kalman tổ hợp đa vật lý để dự báo cường độ bão đánh giá sai số dự báo Phạm vi nghiên cứu Trong khuôn khổ đồ án, phạm vi nghiên cứu giới hạn theo không gian khu vực biển Đông Việt Nam, theo thời gian thời gian có khơng khí lạnh hoạt động Để tiến hành nghiên cứu, đồ án sử dụng ba phương pháp chính: Phương pháp số, phương pháp phân tích đánh giá phương pháp thống kê Ý nghĩa thực tiễn đồ án Trong trình tìm hiểu đồ án, em đọc nhiều tài liệu chun ngành ngồi nước, từ hiểu thêm chun ngành Khí tượng Q trình thực đồ án giúp em rèn luyện kĩ phân tích tổng hợp kiến thức, cách sử dụng từ ngữ nghiên cứu khoa học hết học hỏi phương pháp luận thơng qua cơng trình nghiên cứu liên quan đến đề tài Kết cấu đồ án Với nội dung đồ án có bố cục gồm chương Mở đầu Chương 1: Phương pháp lọc Kalman tổ hợp đa vật lý mơ hình WRF Chương 2: Hồn lưu quy mơ lớn điều kiện có khơng khí lạnh Sự tương tác khơng khí lạnh bão Chương 3: Kết thử nghiệm Kết luận khuyến nghị CHƯƠNG 1: PHƯƠNG PHÁP LỌC KALMAN TỔ HỢP ĐA VẬT LÝ VÀ MƠ HÌNH WRF 1.1 CÁC NGHIÊN CỨU TRONG NƯỚC VÀ NGOÀI NƯỚC Dự báo bão quan tâm nhiều ba khía cạnh dự báo quỹ đạo, dự báo cường độ dự báo bão đổ Trong vài thập kỉ qua, với phát triển vượt bậc khoa học kĩ thuật, việc dự báo quỹ đạo bão có bước tiến đáng kể Tuy nhiên, khơng có cải thiện nhiều kỹ dự báo cường độ bão suốt thời gian Những yếu tố ảnh hưởng đến cường độ bão bao gồm: nhiệt độ mặt nước biển, tương tác với địa chính, địa hình, gió thẳng đứng, độ ẩm môi trường bão, thay đổi cấu trúc tương tác với hệ thống thời tiết khác Trong đó, việc thay đổi cấu trúc mắt bão phức tạp khó dự báo 1.1.1 Các nghiên cứu nước Theo nghiên cứu NOAA (Cơ quan quản trị Khí Đại dương Quốc gia) năm 2017, có cách để dự báo cường độ bão so với dự báo quỹ đạo Dù vòng 20 năm qua, có số cải thiện nhỏ việc dự báo cường độ bão, khoảng từ 10 đến 20% nhờ mơ hình động lực thống kê [3] Một số mơ hình động lực tiêu biểu HWRF (Hurricane Weather Research and Forecast System), HMON (Hurricanes in a Multi-scale Ocean coupled Nonhydrostatic model) Các mô hình hoạt động dựa vào việc giải phương trình điều chỉnh khí (và đại dương) Những mơ hình nhiều hạn chế mạng lưới quan trắc thưa, mô vật lý thay đổi cường độ khơng đầy đủ có vấn đề với biểu diễn cắt Trong nghiên cứu năm 2013, Kiều cộng tiến hành thử nghiệm bão Conson (2010) [2] Để xem xét phương pháp lọc Kalman tổ hợp địa phương biến đổi (Local ensemble transform Kalman filter: LETKF) điều kiện thực, tác giả tiến hành thực thuật toán LETKF kết hợp với mơ hình WRF-ARW phiên phi thủy tĩnh Kết cho thấy, phương pháp tăng cấp nhân (MI) cho dự báo cường độ thành phần dự báo gần tương tự toàn thời gian dự báo, phương pháp đa vật lý (MP) cho biến trình thành phần dự báo có độ mở rộng lớn, nửa số biểu diễn cường độ yếu hơn, nửa lại biểu diễn cường độ mạnh thực tế, nên dự báo xác so với phương pháp tăng cấp nhân 1.1.2 Các nghiên cứu nước Tại Việt Nam có nghiên cứu ứng dụng phương pháp lọc Kalman tổ hợp vào dự báo Trong báo “Ứng dụng phương pháp lọc Kalman tổ hợp vào dự báo quỹ đạo cường độ bão ngày” Trần Tân Tiến, phương pháp lọc Kalman tổ hợp (LETKF) kết hợp với mơ hình WRF nhằm cải thiện việc dự báo quỹ đạo cường độ bão Việt Nam hạn ngày thu kết khả quan [11] Nghiên cứu dự báo cường độ bão phương pháp cho hiệu hẳn phương pháp khác sử dụng Việt Nam Cụ thể dự báo hạn ngày sai số tốc độ gió cực đại bão khoảng 9m/s áp suất cực tiểu khoảng 7mb Trong nghiên cứu khác, GS.TS Trần Tân Tiến tiến hành khảo sát vai trò lọc Kalman tổ hợp đồng hóa số liệu vệ tinh cao khơng mơ hình WRF Cơn bão chọn để thử nghiệm bão Nari, hạn dự báo ngày [12] Trong hầu hết trường hợp dự báo, bão Nari chịu ảnh hưởng mở rộng hoạt động áp cao Tây Thái Bình Dương sang phía tây, mức độ lấn tây hay nhiều tùy trường hợp dự báo, hồn tồn khơng có xuất khơng khí lạnh Và chưa có nghiên cứu dự báo cường độ bão lọc Kalman tổ hợp điều kiện có khơng khí lạnh 1.2 PHƯƠNG PHÁP LỌC KALMAN TỔ HỢP ĐA VẬT LÝ Hiện nay, tin dự báo thời tiết cho kết sai lệch điều kiện ban đầu khơng xác [8], điều đặc thù mơ hình dự báo thời tiết có tính phụ thuộc mạnh vào trường ban đầu (trường đầu vào mơ hình) Một phương pháp làm xác điều kiện ban đầu cho mơ hình đồng hóa số liệu, q trình tạo trường ban đầu tốt cho mơ hình dự báo, dựa mối quan hệ động lực xác suất thống kê Về đồng hóa số liệu chia Bảng 1.2b: Kết dự báo PMIN thử nghiệm MP với thành phần tổ hợp từ 008 đến 014 HẠN DỰ BÁO 008 009 010 011 012 013 014 995,325 995,325 995,325 995,325 995,325 995,325 995,325 987,513 987,985 987,985 988,170 987,045 987,004 987,457 984,566 986,929 986,929 986,493 983,808 983,476 986,675 980,813 983,833 983,833 983,564 979,474 979,552 983,953 12 976,356 983,683 983,683 983,102 976,046 975,721 982,369 15 975,330 982,939 982,939 982,473 974,683 974,513 981,602 18 973,678 981,269 981,269 980,799 971,116 971,337 980,013 21 971,772 981,563 981,563 983,376 970,834 970,545 980,346 24 973,044 980,806 980,806 983,520 972,305 971,717 980,293 27 971,562 981,379 981,379 984,388 971,033 970,594 980,763 30 968,581 981,054 981,054 983,909 969,263 968,940 980,658 33 971,172 981,808 981,808 984,049 969,837 968,175 982,043 36 974,269 985,513 985,513 984,991 972,929 972,067 983,370 39 976,116 985,884 985,884 984,863 975,044 974,694 983,787 42 975,851 985,999 985,999 985,927 975,225 975,413 984,667 45 976,495 988,057 988,057 985,696 975,802 976,703 984,578 48 978,553 986,856 986,856 987,393 975,826 976,599 985,379 51 977,216 988,386 988,386 987,966 973,762 974,808 985,686 54 976,658 988,103 988,103 988,052 972,822 974,480 985,006 57 977,661 988,746 988,746 989,301 973,950 975,408 986,106 60 978,305 991,794 991,794 991,725 975,955 977,581 988,665 63 976,963 990,937 990,937 990,872 976,355 978,163 988,949 66 976,898 991,065 991,065 989,980 990,003 990,760 994,244 69 979,785 993,516 993,516 990,244 993,846 993,846 990,649 72 PL.6 Bảng 1.2c: Kết dự báo PMIN thử nghiệm MP với thành phần tổ hợp từ 015 đến 021 HẠN DỰ BÁO 015 016 017 018 019 020 021 995,325 995,325 995,325 995,325 995,325 995,325 995,325 987,438 987,077 987,037 987,624 987,591 986,936 987,002 986,413 983,962 983,872 985,721 985,726 983,660 983,607 983,663 979,679 978,931 982,990 983,090 978,574 978,497 12 982,322 974,663 974,222 979,142 980,150 973,836 973,750 15 981,534 972,925 973,081 979,705 981,738 971,879 971,879 18 980,385 969,152 968,965 976,418 975,868 966,782 965,361 21 980,789 966,193 967,522 976,746 975,548 964,130 965,597 24 982,245 964,913 964,782 977,983 979,035 963,326 963,813 27 980,258 963,366 962,108 977,232 978,045 962,437 962,004 30 979,940 962,822 961,522 977,538 978,324 963,593 962,914 33 979,481 963,251 959,975 978,680 978,811 964,199 961,987 36 981,462 965,931 961,957 981,058 981,780 967,535 965,203 39 983,324 967,221 963,577 980,968 983,253 970,287 967,738 42 983,614 967,327 963,196 979,584 982,032 970,774 969,341 45 984,714 969,243 965,311 980,207 981,458 971,722 970,489 48 984,934 970,361 965,899 982,998 981,061 972,641 970,819 51 982,820 970,822 965,356 984,455 982,833 972,227 970,773 54 984,493 969,208 965,040 983,184 983,572 972,484 970,616 57 984,017 970,495 965,520 982,847 983,952 973,362 969,705 60 986,064 971,866 966,400 984,677 984,562 973,950 970,655 63 992,601 972,504 967,004 985,530 982,560 973,943 972,430 66 993,531 974,378 990,487 984,566 983,160 977,549 990,045 981,854 984,796 69 72 986,021 PL.7 Bảng 1.3a: Kết dự báo PMIN thử nghiệm PF với thành phần tổ hợp từ 001 đến 007 HẠN DỰ BÁO 001 002 003 004 005 006 007 995,325 995,325 995,325 995,325 995,325 995,325 995,325 986,493 986,491 986,608 986,602 986,693 986,637 986,532 984,578 984,631 984,835 984,765 984,741 984,589 984,683 981,452 981,445 981,568 981,325 981,370 981,214 981,376 12 979,280 979,428 979,525 979,046 979,236 979,245 979,269 15 978,057 978,069 977,872 978,122 977,878 977,990 977,820 18 971,978 972,247 972,122 972,281 971,866 971,906 972,454 21 971,809 972,188 971,948 972,149 972,145 972,124 972,190 24 970,327 970,980 970,313 970,503 970,816 970,539 970,654 27 967,629 968,024 967,562 967,235 966,731 967,395 967,134 30 964,274 964,732 964,392 964,685 964,440 964,680 964,378 33 963,360 963,857 963,066 963,279 962,952 962,969 963,144 36 963,180 963,724 963,212 962,783 963,078 962,818 962,615 39 962,387 963,157 962,310 961,931 961,706 961,784 962,036 42 961,115 960,797 962,129 960,016 960,915 960,661 961,121 45 961,160 959,262 962,484 959,562 961,195 961,367 961,453 48 960,458 959,172 961,503 959,526 961,231 961,018 961,236 51 959,969 959,670 960,708 958,824 960,333 960,107 960,304 54 957,784 957,729 958,997 957,468 959,040 957,163 958,159 57 958,122 957,345 958,073 956,648 958,469 957,179 957,320 60 958,731 957,837 958,072 957,413 959,752 958,432 956,852 63 958,537 957,365 957,513 957,069 959,702 957,812 956,528 66 957,260 956,490 956,664 955,797 958,670 956,487 956,498 69 958,392 957,655 957,302 956,907 959,651 957,224 957,939 72 961,046 960,488 959,094 958,342 961,237 958,225 960,008 PL.8 Bảng 1.3b: Kết dự báo PMIN thử nghiệm PF với thành phần tổ hợp từ 008 đến 014 HẠN DỰ BÁO 008 009 010 011 012 013 014 995,325 995,325 995,325 995,325 995,325 995,325 995,325 986,545 986,526 986,529 986,525 986,679 986,468 986,554 984,600 984,878 984,854 984,655 984,927 984,634 984,762 981,079 981,496 981,477 981,486 981,295 981,334 981,270 12 979,440 979,298 979,236 979,189 979,297 979,102 979,427 15 977,984 977,903 978,062 978,162 978,198 977,957 978,013 18 972,406 972,254 972,196 972,431 972,641 972,220 972,169 21 972,205 972,014 972,172 972,305 972,302 972,190 972,208 24 970,347 970,741 970,667 970,720 970,680 970,375 970,322 27 966,995 967,548 967,827 967,686 967,521 966,986 967,126 30 964,390 964,904 964,410 964,974 964,678 964,965 964,957 33 963,111 963,857 963,458 963,552 963,477 963,636 963,572 36 963,216 963,626 963,220 963,060 963,401 963,009 963,274 39 962,094 962,656 962,289 962,676 962,314 962,011 961,963 42 961,773 961,148 961,474 961,686 961,510 961,137 961,338 45 962,446 961,037 962,098 961,741 961,011 961,305 961,971 48 962,024 960,430 961,582 960,812 960,881 961,392 960,995 51 960,616 959,533 961,177 960,183 960,009 960,541 960,634 54 958,490 957,701 959,702 957,723 958,099 958,568 958,374 57 958,245 956,929 958,478 957,456 957,022 958,113 957,452 60 958,634 957,772 959,021 957,551 957,605 958,888 958,961 63 957,606 957,085 958,286 957,318 956,931 958,751 958,213 66 956,963 956,372 957,691 956,514 956,249 956,453 957,450 69 958,113 957,536 958,847 957,403 957,440 957,079 958,412 72 960,135 959,258 961,037 958,933 959,775 959,04 960,244 PL.9 Bảng 1.3c: Kết dự báo PMIN thử nghiệm PF với thành phần tổ hợp từ 015 đến 021 HẠN DỰ BÁO 015 016 017 018 019 020 021 995,325 995,325 995,325 995,325 995,325 995,325 995,325 986,441 986,400 986,291 986,429 986,530 986,553 986,506 984,647 984,752 984,707 984,649 984,597 984,641 984,586 981,267 981,291 981,403 981,284 981,258 981,252 981,375 12 979,344 979,385 979,266 979,419 979,193 979,197 979,171 15 977,763 978,112 977,958 978,263 977,763 977,782 977,902 18 971,998 972,105 972,290 971,650 972,028 971,945 972,497 21 972,104 972,153 972,066 971,651 972,076 972,114 972,125 24 970,523 970,722 970,795 970,525 970,399 970,375 970,509 27 967,069 967,862 967,389 967,067 967,047 967,400 967,338 30 964,381 964,808 964,353 964,590 964,627 964,456 964,489 33 963,000 963,401 963,131 962,712 962,908 963,080 963,271 36 963,209 963,660 962,603 962,997 963,237 963,068 962,905 39 962,448 963,049 961,370 961,989 962,487 962,924 962,033 42 961,735 961,443 959,949 960,517 961,570 961,625 961,038 45 962,219 962,403 960,169 960,686 961,250 962,082 961,323 48 962,159 961,312 960,615 960,504 960,351 961,435 961,350 51 962,240 959,384 960,020 959,999 959,343 959,398 959,562 54 960,059 957,397 958,027 958,159 956,823 956,334 956,930 57 960,029 955,912 956,821 958,026 955,774 955,290 956,037 60 960,847 957,040 957,440 958,289 956,388 957,246 956,524 63 960,070 956,962 956,861 957,848 955,618 957,115 955,479 66 958,965 956,311 956,792 957,166 954,605 956,200 954,729 69 960,495 957,356 958,704 958,935 955,789 956,929 956,194 72 962,618 959,46 960,346 961,141 957,674 957,864 957,718 PL.10 PHỤ LỤC 2: KẾT QUẢ DỰ BÁO TỐC ĐỘ GIÓ CỰC ĐẠI GẦN TÂM CỦA CÁC THÀNH PHẦN TỔ HỢP TRONG THỬ NGHIỆM MI, MP, PF PL.11 Bảng 2.1a: Kết dự báo VMAX thử nghiệm MI với thành phần tổ hợp từ 001 đến 007 HẠN DỰ BÁO 001 002 003 004 005 006 007 28,393 28,393 28,393 28,393 28,393 28,393 28,393 28,100 28,055 28,152 27,877 28,166 28,097 28,297 30,890 31,035 31,027 31,426 30,968 31,398 31,071 32,359 33,598 32,897 33,537 32,860 33,615 32,970 12 33,476 33,180 33,858 33,455 33,389 33,463 33,069 15 36,508 36,163 36,305 36,827 36,833 36,669 36,795 18 42,797 42,191 42,033 42,959 42,390 42,128 42,334 21 42,823 42,621 42,378 44,263 42,793 42,663 43,357 24 42,769 41,652 41,679 43,226 43,086 42,575 42,193 27 48,295 47,597 48,280 49,031 48,665 48,550 48,746 30 48,057 47,539 50,140 48,557 49,945 48,369 50,651 33 49,127 49,592 50,303 50,213 49,927 50,208 50,651 36 49,224 49,875 49,719 49,366 49,109 49,410 49,758 39 48,139 48,484 48,763 49,712 48,670 49,274 48,184 42 49,233 50,751 50,947 50,382 49,752 49,942 49,076 45 50,134 48,380 46,331 48,695 49,022 48,412 49,540 48 48,683 46,241 49,017 47,009 44,744 49,653 47,191 51 46,807 44,518 43,091 46,078 45,920 43,942 46,222 54 46,317 44,126 44,289 43,460 44,965 47,550 44,566 57 46,579 47,022 46,554 46,267 46,389 48,027 47,335 60 47,145 47,000 46,649 50,847 47,766 48,325 48,985 63 47,345 48,762 47,340 46,584 49,889 48,935 48,884 66 46,131 47,609 47,246 49,588 47,462 49,549 48,197 69 45,156 47,880 46,216 47,999 47,225 49,008 47,726 72 45,169 45,433 46,478 47,696 47,071 46,416 46,339 PL.12 Bảng 2.1b: Kết dự báo VMAX thử nghiệm MI với thành phần tổ hợp từ 008 đến 014 HẠN DỰ BÁO 008 009 010 011 012 013 014 28,393 28,393 28,393 28,393 28,393 28,393 28,393 28,372 27,998 28,287 28,134 28,178 28,258 27,950 31,446 30,969 31,118 31,460 31,113 31,249 31,128 33,130 33,078 32,954 33,120 32,847 32,349 33,195 12 33,601 33,586 33,565 33,682 33,332 33,230 33,804 15 36,182 35,955 36,493 36,848 36,334 36,639 35,656 18 41,933 42,268 42,375 43,444 42,783 41,485 42,283 21 42,770 42,465 42,140 42,927 43,796 43,445 42,948 24 42,360 42,060 42,943 41,450 41,687 41,196 41,928 27 48,244 47,914 49,598 48,375 47,459 49,256 48,735 30 49,703 49,477 49,184 49,062 48,734 49,187 49,422 33 49,584 49,721 49,640 48,641 49,691 49,559 49,683 36 49,498 49,518 50,457 48,194 49,824 48,900 49,882 39 50,220 49,313 48,415 48,188 49,844 46,312 48,573 42 50,539 50,009 50,918 50,789 49,817 49,750 49,116 45 50,154 46,480 48,701 47,448 48,438 47,978 48,360 48 46,297 48,547 49,449 49,027 48,048 46,973 46,941 51 45,379 42,675 47,433 44,400 44,789 45,224 47,290 54 45,045 47,072 47,716 47,202 46,375 46,891 46,121 57 49,406 45,437 47,701 47,716 46,169 47,358 48,343 60 49,890 48,372 47,113 49,767 49,009 49,715 47,664 63 48,712 46,483 45,603 47,400 45,473 46,270 46,871 66 46,579 48,494 49,563 48,418 47,974 46,778 46,960 69 47,551 46,886 46,180 47,484 46,920 46,350 46,873 72 45,758 48,072 45,776 44,861 47,464 45,482 47,156 PL.13 Bảng 2.1c: Kết dự báo VMAX thử nghiệm MI với thành phần tổ hợp từ 015 đến 021 HẠN DỰ BÁO 015 016 017 018 019 020 021 28,393 28,393 28,393 28,393 28,393 28,393 28,393 28,430 28,030 28,238 28,237 28,241 28,286 28,157 31,252 31,093 31,161 31,380 31,140 31,109 31,060 33,128 33,292 32,650 33,408 32,967 33,527 33,043 12 33,574 33,809 33,533 33,614 33,639 33,860 33,285 15 36,104 36,419 36,013 36,574 36,465 36,976 36,374 18 42,561 42,184 42,420 42,900 41,645 42,172 41,654 21 43,432 43,421 43,278 42,763 43,030 42,784 42,444 24 42,255 42,605 41,300 42,556 40,368 43,240 42,444 27 49,406 48,774 48,817 47,838 47,031 48,228 48,562 30 49,172 48,247 48,854 49,165 48,983 48,983 48,630 33 49,237 49,273 49,809 48,926 49,140 49,642 50,920 36 49,800 50,438 49,496 48,782 50,241 48,866 48,668 39 49,148 49,359 49,002 48,956 47,608 48,699 49,154 42 48,940 51,168 49,498 47,828 48,795 48,498 49,738 45 47,439 49,345 49,049 48,132 48,357 48,580 46,482 48 45,900 47,722 47,126 48,949 46,323 46,870 49,536 51 46,993 45,459 44,696 46,261 47,217 47,689 44,443 54 44,144 43,431 47,965 46,852 46,500 47,805 47,249 57 48,602 45,110 46,687 46,927 47,252 49,077 46,981 60 48,124 46,881 48,737 48,095 48,854 47,664 48,589 63 49,189 48,280 47,929 46,868 47,428 47,082 46,347 66 48,485 47,031 49,510 47,553 46,487 47,880 46,193 69 47,787 48,419 49,487 47,601 46,423 48,741 47,120 72 46,697 47,411 47,870 46,070 46,503 46,965 46,677 PL.14 Bảng 2.2a: Kết dự báo VMAX thử nghiệm MP với thành phần tổ hợp từ 001 đến 007 HẠN DỰ BÁO 001 002 003 004 005 006 007 28,393 28,393 28,393 28,393 28,393 28,393 28,393 27,866 28,360 28,333 30,308 30,318 29,713 29,517 31,356 30,049 30,822 33,790 33,959 32,606 33,461 33,016 30,388 30,210 39,389 42,023 37,806 38,288 12 33,632 34,140 35,116 40,971 42,328 38,662 37,443 15 36,192 33,143 31,600 43,533 43,354 36,774 35,524 18 42,545 35,262 32,942 51,230 51,671 41,792 46,357 21 40,468 34,164 37,823 48,178 48,806 45,446 40,612 24 44,171 40,662 42,445 56,569 57,549 43,097 41,562 27 48,031 44,171 42,783 54,495 55,367 47,222 47,145 30 48,290 47,764 42,614 53,112 53,330 49,070 45,648 33 49,641 42,754 40,365 56,156 55,523 46,696 44,284 36 50,849 42,488 43,505 52,306 56,840 45,110 44,202 39 50,424 42,523 41,660 51,027 49,797 43,659 44,467 42 50,256 44,547 44,512 50,741 53,150 46,858 44,518 45 48,077 43,176 43,962 51,967 51,462 45,543 39,123 48 47,509 39,087 40,407 49,215 49,972 39,095 42,375 51 45,180 42,049 45,262 45,740 47,138 45,709 47,376 54 46,339 45,235 44,227 46,410 46,249 43,934 41,482 57 48,599 47,379 46,077 41,458 43,673 41,419 36,946 60 48,220 41,684 44,503 44,861 40,938 40,974 33,301 63 48,537 35,365 42,098 38,272 38,387 23,798 24,974 66 49,146 40,163 39,487 27,203 26,352 23,556 24,799 69 47,068 40,295 41,169 24,030 24,004 72 47,266 36,869 42,663 PL.15 Bảng 2.2b: Kết dự báo VMAX thử nghiệm MP với thành phần tổ hợp từ 008 đến 014 HẠN DỰ BÁO 008 009 010 011 012 013 014 28,393 28,393 28,393 28,393 28,393 28,393 28,393 29,995 29,808 29,039 29,092 30,189 30,247 29,529 33,258 33,598 32,319 32,231 33,365 33,319 32,103 40,192 40,111 38,335 38,588 38,906 40,719 34,757 12 44,510 46,754 38,422 34,968 42,065 41,760 37,894 15 39,878 40,578 34,164 37,450 42,667 40,832 36,663 18 41,516 39,883 40,231 39,991 42,450 40,158 38,306 21 46,321 43,902 36,327 37,665 46,670 46,149 40,226 24 40,952 42,822 39,023 31,540 44,512 45,624 43,213 27 42,149 44,991 36,388 33,789 45,436 48,867 38,851 30 50,705 54,452 37,857 33,193 50,697 53,615 38,083 33 51,215 49,243 33,772 34,224 52,652 53,122 37,760 36 46,620 45,923 31,171 35,206 49,096 50,236 36,113 39 41,561 41,473 33,247 32,255 46,483 50,782 33,937 42 44,218 45,377 31,875 30,799 44,595 42,700 32,395 45 45,270 43,046 28,214 32,293 39,980 41,706 34,850 48 44,796 44,620 32,196 28,902 41,004 40,702 34,045 51 45,560 43,350 31,145 25,373 42,917 42,752 30,755 54 45,766 40,455 29,914 25,686 42,234 42,271 33,279 57 40,625 39,685 26,854 23,852 39,967 40,578 29,294 60 38,570 40,455 25,480 23,973 37,373 36,326 31,607 63 40,076 39,538 25,699 24,860 36,991 35,091 27,850 66 38,663 35,530 27,563 24,944 27,873 27,316 24,572 69 33,379 33,754 25,119 29,603 27,537 26,199 72 PL.16 Bảng 2.2c: Kết dự báo VMAX thử nghiệm MP với thành phần tổ hợp từ 015 đến 021 HẠN DỰ BÁO 015 016 017 018 019 020 021 28,393 28,393 28,393 28,393 28,393 28,393 28,393 29,332 29,535 29,431 29,329 29,139 30,249 30,375 32,141 32,961 33,036 32,546 33,859 33,254 33,151 35,608 41,497 41,433 33,159 33,262 41,264 42,446 12 36,700 44,187 44,685 39,486 40,122 43,241 43,643 15 39,121 43,537 44,995 35,909 36,219 43,641 42,986 18 38,366 46,259 46,297 43,141 46,836 47,265 50,737 21 38,876 49,742 50,436 42,427 42,111 50,136 48,409 24 38,385 49,967 51,949 41,707 38,694 52,217 53,796 27 41,270 54,663 55,005 41,062 37,890 55,152 55,026 30 40,111 54,222 53,783 39,830 40,407 55,630 54,612 33 38,229 52,821 53,760 38,768 38,995 54,293 57,111 36 39,468 53,061 53,400 35,095 30,352 52,631 54,523 39 35,841 53,916 53,000 37,249 29,806 49,332 52,395 42 32,072 52,100 52,843 37,799 33,653 48,598 48,076 45 32,530 49,112 50,342 42,456 37,738 48,143 47,335 48 32,270 48,265 48,489 40,734 37,427 48,890 45,595 51 29,983 45,254 48,369 30,326 30,396 46,959 44,152 54 32,311 44,362 47,014 35,164 28,702 48,334 43,007 57 33,487 43,229 45,380 31,813 32,886 38,683 43,463 60 30,083 42,753 44,051 34,622 35,779 40,656 40,952 63 25,548 40,058 42,064 30,009 36,352 41,573 38,792 66 23,891 37,401 26,053 29,774 30,225 38,146 27,756 31,535 26,848 69 72 28,989 PL.17 Bảng 2.3a: Kết dự báo VMAX thử nghiệm PF với thành phần tổ hợp từ 001 đến 007 HẠN DỰ BÁO 001 002 003 004 005 006 007 28,393 28,393 28,393 28,393 28,393 28,393 28,393 28,178 28,270 28,317 28,306 28,316 28,313 28,326 31,158 31,169 31,334 31,491 31,132 31,435 31,474 32,815 32,885 32,935 33,000 33,027 32,977 32,865 12 33,231 33,390 33,276 33,450 33,218 33,771 33,400 15 36,669 36,255 36,394 36,459 36,562 36,783 36,193 18 42,076 42,178 41,325 40,936 41,876 41,410 40,927 21 41,157 41,264 40,643 40,611 41,342 41,206 41,532 24 42,123 41,792 42,234 42,642 42,270 42,494 42,380 27 45,904 46,875 47,599 46,891 48,762 47,781 48,341 30 49,300 49,672 48,778 48,634 48,360 48,211 49,629 33 49,799 48,328 50,481 50,080 50,613 50,284 50,302 36 50,651 48,871 51,257 50,008 49,823 49,861 49,814 39 49,428 49,347 50,249 49,219 50,661 50,408 49,908 42 50,477 50,784 50,702 50,451 51,295 51,111 49,879 45 48,427 50,647 46,935 51,365 45,861 48,813 47,663 48 48,631 48,919 45,698 49,717 47,309 48,362 49,018 51 46,843 46,941 45,786 48,475 47,158 45,570 46,135 54 47,379 47,754 46,779 47,439 44,171 49,203 45,486 57 45,887 47,641 47,959 49,371 48,055 47,488 47,549 60 45,674 47,047 48,698 46,614 47,158 47,173 47,774 63 47,517 46,364 46,050 47,513 46,710 47,685 47,965 66 47,079 49,002 47,150 47,480 46,202 48,320 48,123 69 46,612 46,751 47,504 47,509 47,616 47,967 46,149 72 45,432 45,199 47,063 46,208 47,119 47,052 45,437 PL.18 Bảng 2.3b: Kết dự báo VMAX thử nghiệm MP với thành phần tổ hợp từ 008 đến 014 HẠN DỰ BÁO 008 009 010 011 012 013 014 28,393 28,393 28,393 28,393 28,393 28,393 28,393 28,319 28,240 28,222 28,317 28,221 28,322 28,229 31,339 31,091 31,257 31,268 31,147 31,176 31,081 32,872 32,734 33,305 32,910 32,794 32,707 32,571 12 33,410 33,231 33,592 33,086 33,291 33,362 33,220 15 36,815 36,403 36,287 36,290 35,602 36,848 36,500 18 41,106 41,494 42,750 41,413 41,682 41,380 41,298 21 41,544 41,285 41,787 40,758 41,170 41,641 41,292 24 42,687 41,844 42,355 42,919 42,241 43,100 42,499 27 48,924 47,103 47,043 47,172 47,820 47,565 48,075 30 49,326 48,470 49,704 48,398 48,979 48,483 48,158 33 51,012 49,588 50,100 50,517 50,292 50,441 50,063 36 49,600 49,469 49,585 50,379 49,958 50,482 50,878 39 49,401 49,071 48,804 49,522 49,590 49,909 50,273 42 52,254 50,765 50,617 50,271 50,451 50,610 50,395 45 48,655 47,609 47,447 47,995 48,410 49,477 45,418 48 46,723 47,425 47,134 48,651 47,250 48,753 49,346 51 45,283 46,902 46,253 46,206 44,724 46,110 44,300 54 46,350 46,163 44,786 46,902 46,053 46,279 47,123 57 46,771 46,876 46,948 45,961 47,044 47,838 46,453 60 47,570 48,269 46,986 49,581 48,482 46,555 46,875 63 48,315 47,806 47,208 46,145 46,823 45,638 47,492 66 46,645 47,159 46,643 48,063 49,516 49,744 47,844 69 46,952 46,786 46,981 47,122 47,001 47,623 47,503 72 45,871 45,835 46,275 46,257 46,562 46,749 45,673 PL.19 Bảng 2.3c: Kết dự báo VMAX thử nghiệm PF với thành phần tổ hợp từ 015 đến 021 HẠN DỰ BÁO 015 016 017 018 019 020 021 28,393 28,393 28,393 28,393 28,393 28,393 28,393 28,403 28,412 28,359 28,238 28,252 28,242 28,328 31,545 31,333 31,144 31,359 31,387 31,369 31,235 33,115 33,489 32,531 32,985 33,036 32,873 32,802 12 33,351 33,470 33,282 33,537 33,369 34,016 33,532 15 36,520 36,478 36,486 36,758 36,694 36,026 36,641 18 41,762 41,437 41,703 42,193 41,945 41,407 41,497 21 40,800 40,676 41,179 40,598 41,875 40,819 41,877 24 42,144 41,580 41,979 41,728 42,092 42,005 42,591 27 48,578 47,117 48,569 47,419 48,160 48,371 47,996 30 49,012 49,167 48,784 48,218 48,768 49,064 49,811 33 50,744 49,845 50,280 50,530 50,870 50,278 50,184 36 50,741 49,703 50,387 50,355 50,213 49,653 50,486 39 48,906 49,615 49,387 49,764 50,782 50,176 49,910 42 50,536 50,456 50,978 50,306 50,303 49,658 51,254 45 47,232 45,070 49,398 47,860 49,018 46,540 49,278 48 50,019 46,943 49,123 45,852 47,656 47,077 47,286 51 43,887 48,148 47,574 48,123 48,158 46,627 45,479 54 46,685 45,744 44,465 45,034 45,585 47,972 47,459 57 46,118 49,465 44,185 45,087 48,401 49,381 48,528 60 44,648 48,322 48,623 47,082 48,457 48,175 48,327 63 46,538 46,770 47,827 46,239 47,874 47,270 48,429 66 46,997 48,000 46,564 49,865 50,223 48,506 48,487 69 45,627 47,479 46,361 48,168 46,931 46,611 47,640 72 45,426 46,254 46,474 45,626 46,832 48,035 47,085 PL.20 ... trường hợp dự báo, hoàn toàn khơng có xuất khơng khí lạnh Và chưa có nghiên cứu dự báo cường độ bão lọc Kalman tổ hợp điều kiện có khơng khí lạnh 1.2 PHƯƠNG PHÁP LỌC KALMAN TỔ HỢP ĐA VẬT LÝ Hiện... phương pháp lọc Kalman tổ hợp đa vật lý mơ hình WRF, hồn lưu qui mơ lớn trường hợp có bão hoạt động điều kiện khơng khí lạnh, chọn bão để thử nghiệm, áp dụng lọc Kalman tổ hợp đa vật lý để dự báo cường. .. dụng phương pháp lọc Kalman tổ hợp vào dự báo Trong báo “Ứng dụng phương pháp lọc Kalman tổ hợp vào dự báo quỹ đạo cường độ bão ngày” Trần Tân Tiến, phương pháp lọc Kalman tổ hợp (LETKF) kết hợp