Nghiên cứu ứng dụng một số phản ứng hạt nhân gây bởi chùm hạt tích điện trên máy gia tốc tĩnh điện trong phân tích Nghiên cứu ứng dụng một số phản ứng hạt nhân gây bởi chùm hạt tích điện trên máy gia tốc tĩnh điện trong phân tích luận văn tốt nghiệp thạc sĩ
ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN -*** - VÕ VĂN HÒA NGHIÊN CỨU ỨNG DỤNG PHƯƠNG PHÁP THỐNG KÊ SAU MƠ HÌNH TỔ HỢP (EMOS) VÀO DỰ BÁO NHIỆT ĐỘ VÀ ĐIỂM SƯƠNG BỀ MẶT Ở VIỆT NAM LUẬN ÁN TIẾN SĨ NGÀNH KHÍ TƯỢNG VÀ KHÍ HẬU HỌC Hà Nội - 2015 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN -*** - VÕ VĂN HÒA NGHIÊN CỨU ỨNG DỤNG PHƯƠNG PHÁP THỐNG KÊ SAU MƠ HÌNH TỔ HỢP (EMOS) VÀO DỰ BÁO NHIỆT ĐỘ VÀ ĐIỂM SƯƠNG BỀ MẶT Ở VIỆT NAM Chuyên ngành: Khí tượng khí hậu học Mã số: 62.44.02.22 LUẬN ÁN TIẾN SĨ NGÀNH KHÍ TƯỢNG VÀ KHÍ HẬU HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC TS Bùi Minh Tăng GS.TS Phan Văn Tân Hà Nội - 2015 Lời cam đoan Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết trình bày luận án trung thực chưa công bố cơng trình khác Tác giả luận án Võ Văn Hịa Lời cảm ơn Luận án hồn thành Trung tâm Dự báo khí tượng thủy văn Trung ương Khoa Khí tượng Thủy văn Hải dương học, Trường Đại học Khoa học Tự Nhiên, Hà Nội Tác giả xin bày tỏ lòng biết ơn chân thành tới hai nhà khoa học hết lòng động viên, tận tình giúp đỡ quan tâm tới bước nghiên cứu luận án Để thực luận án, tác giả giúp đỡ thời gian điều kiện nghiên cứu thuận lợi từ Trung tâm Dự báo khí tượng thủy văn Trung ương, nơi tác giả hỗ trợ tài liệu, số liệu nghiên cứu, trang thiết bị tính tốn lưu trữ số liệu Lời tri ân tác giả muốn gửi tới đồng nghiệp Trung tâm Dự báo khí tượng thủy văn Trung ương chia sẻ, giúp đỡ góp ý số nội dung nghiên cứu luận án Tác giả không quên quan tâm, chăm sóc, sẻ chia buồn vui giúp đỡ qua bao khó khăn người bạn đời hai nhỏ Lòng biết ơn sâu nặng tác giả xin gửi cha mẹ, người ban cho tác giả sống, dưỡng nuôi suốt thời thơ ấu định hướng khoa học đường theo đuổi suốt đời tác giả Tác giả Mục lục Lời cam đoan Lời cảm ơn Mục lục Danh sách ký hiệu chữ viết tắt Danh mục bảng Danh mục hình vẽ, đồ thị 12 MỞ ĐẦU 14 CHƯƠNG TỔNG QUAN VỀ CÁC PHƯƠNG PHÁP NGHIÊN CỨU HIỆU CHỈNH THỐNG KÊ CHO DỰ BÁO TỔ HỢP 21 1.1 Khái niệm vai trò dự báo tổ hợp 21 1.2 Sự cần thiết phải hiệu chỉnh dự báo tổ hợp 24 1.3 Tổng quan nghiên cứu hiệu chỉnh thống kê cho dự báo tổ hợp 27 1.3.1 Ngoài nước 27 1.3.2 Trong nước 34 CHƯƠNG MÔ TẢ SỐ LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 37 2.1 Đặt toán 37 2.2 Nghiên cứu lựa chọn ứng dụng phương pháp thống kê để nâng cao chất lượng dự báo trung bình tổ hợp xác suất từ hệ thống SREPS 40 2.2.1 Các phương pháp thống kê tác động đến dự báo trung bình tổ hợp 40 2.2.2 Các phương pháp thống kê tác động đến dự báo trung bình tổ hợp dự báo xác suất 43 2.2.3 Ứng dụng phương pháp EMOS cho hệ thống SREPS 49 2.3 Mô tả tập số liệu nghiên cứu 53 2.4 Phương pháp xử lý số liệu 58 2.5 Phương pháp đánh giá 60 2.5.1 Các số đánh giá chất lượng dự báo trung bình tổ hợp 61 2.5.2 Các số đánh giá chất lượng dự báo xác suất 62 CHƯƠNG KẾT QUẢ THỬ NGHIỆM VÀ ĐÁNH GIÁ 67 3.1 Đánh giá khả khử sai số hệ thống 67 3.2 Đánh giá khả cải thiện chất lượng dự báo trung bình tổ hợp 74 3.2.1 Yếu tố nhiệt độ khơng khí 74 3.2.2 Yếu tố nhiệt độ điểm sương 78 3.2.3 Yếu tố nhiệt độ tối cao ngày 81 3.2.4 Yếu tố nhiệt độ tối thấp ngày 83 3.2.5 Nhận định chung 86 3.3 Đánh giá khả cải thiện chất lượng dự báo xác suất 88 3.3.1 Yếu tố nhiệt độ khơng khí 89 3.3.2 Yếu tố nhiệt độ điểm sương 94 3.3.3 Yếu tố nhiệt độ tối cao ngày 98 3.3.4 Yếu tố nhiệt độ tối thấp ngày 102 3.3.5 Nhận định chung 106 3.4 Xây dựng qui trình dự báo nghiệp vụ đánh giá kết thử nghiệm nghiệp vụ từ năm 2011-2013 109 3.4.1 Xây dựng qui trình dự báo nghiệp vụ cho hệ thống EMOS dựa phương pháp NGR_EM 109 3.4.2 Kết đánh giá thử nghiệm nghiệp vụ từ năm 2011 đến 2013 112 KẾT LUẬN 122 KIẾN NGHỊ VỀ CÁC HƯỚNG NGHIÊN CỨU TIẾP THEO 124 DANH MỤC CÁC CƠNG TRÌNH NGHIÊN CỨU CỦA TÁC GIẢ CĨ LIÊN QUAN ĐẾN LUẬN ÁN 125 TÀI LIỆU THAM KHẢO 127 PHỤ LỤC Danh sách 174 trạm sử dụng nghiên cứu 139 PHỤ LỤC Kết đánh giá cho T2m hạn dự báo 48h 72h 142 PHỤ LỤC Kết đánh giá cho Td2m hạn dự báo 48h 72h 148 PHỤ LỤC Kết đánh giá cho Tmax hạn dự báo 48h 72h 154 PHỤ LỤC Kết đánh giá cho Tmin hạn dự báo 48h 72h 160 PHỤ LỤC Kết đánh giá thử nghiệm nghiệp vụ từ 2011-2013 166 Danh sách ký hiệu chữ viết tắt Hệ số hồi qui thứ i trọng số thứ i ANN Mạng thần kinh nhân tạo BCES Hiệu chỉnh sai số phương pháp trọng số hàm mũ BCLR Hiệu chỉnh sai số phương pháp hồi quy tuyến tính BCMA Hiệu chỉnh sai số phương pháp trung bình trượt BMA Phương pháp trung bình Bayes BIAS Sai số hệ thống BS Chỉ số đánh giá xác suất BSS Chỉ số đánh giá kỹ dự báo xác suất CMC Cơ quan khí tượng Canađa CRPS Chỉ số xác suất hạng liên tục DBTH Dự báo tổ hợp DWD Tổng cục thời tiết Cộng hoà Liên bang Đức ECMWF Trung tâm dự báo thời tiết hạn vừa Châu Âu EMES Trung bình tổ hợp có trọng số tính theo hàm mũ EMMV Trung bình tổ hợp có trọng số tính theo phương sai sai số EMOS Hiệu chỉnh thống kê cho dự báo tổ hợp F Yếu tố dự báo Fi Thành phần dự báo tổ hợp thứ i GEM Mơ hình tồn cầu CMC GME Mơ hình tồn cầu DWD GFS Mơ hình tồn cầu NCEP GSM Mơ hình tồn cầu JMA HTTH Hệ thống dự báo tổ hợp JMA Cơ quan khí tượng Nhật Bản KF Phương pháp lọc Kalman KTTV Khí Tượng Thủy Văn LR Hồi quy Logistic ME Sai số trung bình MOS Thống kê sau mơ hình N Dung lượng mẫu tổng số dự báo thành phần NCEP Trung tâm dự báo môi trường quốc gia Mỹ NOGAPS Mơ hình tồn cầu Hải qn Mỹ (US Navy) nnk Những người khác NGR_EM Hồi quy Gauss không kỳ vọng quan hệ tuyến tính với trung bình tổ hợp NGR_ER Hồi quy Gauss khơng kỳ vọng trung bình có trọng số dự báo thành phần NGR_EP Tương tự NGR_ER có loại bỏ trọng số âm NWP Dự báo thời tiết số trị O Giá trị quan trắc Raw Dự báo trung bình tổ hợp xác suất trực tiếp từ hệ thống dự báo tổ hợp thời tiết hạn ngắn (SREPS) RMSE Sai số quân phương SREPS Hệ thống dự báo tổ hợp thời tiết hạn ngắn TBTH Trung bình tổ hợp TTDBTƯ Trung Tâm Dự Báo khí tượng thủy văn Trung Ương Td2m Nhiệt độ điểm sương độ cao mét Tmax Nhiệt độ tối cao ngày Tmin Nhiệt độ tối thấp ngày T2m Nhiệt độ khơng khí độ cao mét Danh mục bảng Bảng 2.1: Các kỹ thuật thống kê thử nghiệm để nâng cao chất lượng dự báo TBTH xác suất hệ thống SREPS (dấu tích ngụ ý có tác động) 50 Bảng 2.2: Cấu hình chi tiết mơ hình NWP khu vực lựa chọn để xây dựng hệ thống SREPS cho khu vực Việt Nam 55 Bảng 2.3: Danh sách nguồn số liệu thu thập phục vụ nghiên cứu 57 Bảng 3.1: Kết tính toán số đánh giá ME (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố T2m với hạn dự báo 24h (các giá trị bôi đậm ngụ ý không đem lại cải thiện sai số) 75 Bảng 3.2: Kết tính tốn số đánh giá MAE (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố T2m với hạn dự báo 24h 76 Bảng 3.3: Kết tính tốn số đánh giá RMSE (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố T2m với hạn dự báo 24h 77 Bảng 3.4: Kết tính tốn số đánh giá ME (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố Td2m với hạn dự báo 24h 79 Bảng 3.5: Kết tính tốn số đánh giá MAE (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố Td2m với hạn dự báo 24h 79 Bảng 3.6: Kết tính tốn số đánh giá RMSE (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố Td2m với hạn dự báo 24h 80 Bảng 3.7: Kết tính tốn số đánh giá ME (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố Tmax với hạn dự báo 24h 81 Bảng 3.8: Kết tính tốn số đánh giá MAE (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố Tmax với hạn dự báo 24h 82 Bảng 3.9: Kết tính toán số đánh giá RMSE (độ C) dự báo TBTH từ Raw 10 phương pháp EMOS cho yếu tố Tmax với hạn dự báo 24h 83 Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCLR 2.75 2.81 2.64 2.65 2.83 2.31 1.75 1.99 1.41 BCES 2.74 2.77 2.63 2.69 2.81 2.14 1.7 1.91 1.46 BCKF 2.94 2.84 2.82 2.85 3.19 2.43 1.97 1.66 NGR_ER 4.19 4.12 3.86 4.27 4.34 3.35 2.52 2.81 2.11 NGR_EP 2.71 2.7 2.65 2.55 2.68 2.25 1.66 2.02 1.51 NGR_EM 2.72 2.79 2.64 2.59 2.74 2.21 1.72 1.95 1.45 Raw 4.79 4.31 3.23 2.67 3.16 3.01 2.95 3.19 2.37 Bảng P4.7: Kết tính tốn số đánh giá CRPS dự báo xác suất từ Raw phương pháp EMOS cho yếu tố Tmax với hạn dự báo 48h (các giá trị bôi đậm ngụ ý không đem lại cải thiện chất lượng dự báo) Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 1.41 1.41 1.35 1.37 1.48 1.2 0.94 1.08 0.81 BCLR 1.65 1.68 1.52 1.59 1.73 1.36 1.03 1.23 0.88 BCES 1.51 1.51 1.44 1.48 1.57 1.22 0.95 1.07 0.8 BCKF 1.65 1.67 1.58 1.6 1.82 1.33 1.06 1.13 0.9 NGR_ER 2.44 2.52 2.35 2.47 2.49 1.98 1.5 1.7 1.29 NGR_EP 1.39 1.41 1.36 1.35 1.45 1.23 0.92 1.11 0.87 NGR_EM 1.37 1.41 1.33 1.35 1.43 1.18 0.9 1.07 0.77 Raw 3.36 2.79 2.08 1.41 1.9 2.03 2.09 1.48 Bảng P4.8: Tương tự bảng P4.7 cho hạn dự báo 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 1.5 1.53 1.46 1.47 1.55 1.21 0.96 1.08 0.83 BCLR 1.75 1.82 1.65 1.7 1.83 1.43 1.08 1.26 0.9 BCES 1.61 1.62 1.56 1.56 1.66 1.24 0.98 1.1 0.83 157 BCKF 1.74 1.73 1.7 1.69 1.95 1.4 1.12 1.21 0.94 NGR_ER 2.64 2.62 2.38 2.64 2.72 2.12 1.58 1.78 1.35 NGR_EP 1.52 1.52 1.48 1.44 1.52 1.27 0.93 1.14 0.85 NGR_EM 1.49 1.56 1.44 1.44 1.52 1.22 0.93 1.08 0.79 Raw 3.3 2.78 2.1 1.49 1.91 1.98 2.02 2.13 1.51 Bảng P4.9: Tương tự bảng P4.7 cho số Độ phủ 90.48% Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 0.68 0.71 0.65 0.72 0.66 0.66 0.68 0.73 0.77 BCLR 0.38 0.38 0.4 0.38 0.4 0.49 0.45 0.49 0.45 BCES 0.63 0.67 0.61 0.69 0.63 0.64 0.67 0.72 0.76 BCKF 0.55 0.52 0.54 0.54 0.47 0.65 0.77 0.74 0.84 NGR_ER 0.46 0.45 0.44 0.44 0.45 0.44 0.44 0.45 0.45 NGR_EP 0.82 0.81 0.8 0.81 0.81 0.77 0.77 0.75 0.75 NGR_EM 0.87 0.85 0.84 0.83 0.86 0.82 0.82 0.83 0.84 Raw 0.32 0.47 0.48 0.77 0.62 0.44 0.45 0.45 0.65 Bảng P4.10: Tương tự bảng P4.7 cho số Độ phủ 90.48% hạn 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 0.69 0.71 0.67 0.75 0.68 0.69 0.7 0.76 0.77 BCLR 0.38 0.37 0.39 0.37 0.39 0.49 0.43 0.48 0.44 BCES 0.64 0.68 0.63 0.71 0.67 0.68 0.69 0.74 0.78 BCKF 0.58 0.57 0.57 0.58 0.52 0.69 0.78 0.76 0.86 NGR_ER 0.44 0.43 0.44 0.42 0.44 0.43 0.43 0.44 0.43 NGR_EP 0.8 0.79 0.78 0.8 0.81 0.76 0.77 0.74 0.75 NGR_EM 0.85 0.82 0.82 0.82 0.84 0.82 0.82 0.83 0.83 Raw 0.35 0.5 0.51 0.8 0.65 0.49 0.47 0.45 0.65 158 Bảng P4.11: Tương tự bảng P4.7 cho số Độ rộng 90.48% Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 4.61 4.96 4.34 5.27 4.89 3.92 3.39 4.04 3.84 BCLR 2.46 2.39 2.32 2.34 2.63 2.72 1.78 2.52 1.54 BCES 4.55 4.9 4.3 5.22 4.9 3.93 3.36 4.05 3.79 BCKF 4.62 4.59 4.52 4.62 4.34 4.67 5.59 5.3 6.08 NGR_ER 4.02 3.95 3.55 3.83 4.08 3.25 2.35 2.91 2.18 NGR_EP 6.48 6.26 5.83 6.12 6.62 5.01 3.71 4.33 3.19 NGR_EM 7.34 7.08 6.42 6.69 7.42 5.6 4.18 5.07 3.71 Raw 5.58 6.02 4.99 6.23 5.83 4.47 4.27 4.85 5.19 Bảng P4.12: Tương tự bảng P4.7 cho số Độ rộng 90.48% hạn 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 5.15 5.56 4.88 5.98 5.57 4.31 3.68 4.32 4.07 BCLR 2.51 2.42 2.37 2.36 2.78 2.84 1.78 2.5 1.53 BCES 5.11 5.54 4.9 5.99 5.64 4.39 3.7 4.4 4.09 BCKF 5.73 5.61 5.71 6.0 5.55 5.68 6.27 6.19 6.62 NGR_ER 4.16 3.91 3.58 3.91 4.33 3.32 2.45 2.86 2.16 NGR_EP 6.62 6.31 5.97 6.29 7.04 5.08 3.69 4.23 3.2 NGR_EM 7.4 7.18 6.64 6.9 7.8 5.74 4.24 5.04 3.76 Raw 6.12 6.58 5.5 6.8 6.44 4.86 4.49 5.12 5.31 159 PHỤ LỤC Kết đánh giá cho Tmin hạn dự báo 48h 72h Bảng P5.1: Kết tính tốn số đánh giá ME dự báo xác suất từ Raw 10 phương pháp EMOS cho yếu tố Tmax với hạn dự báo 48h (các giá trị bôi đậm ngụ ý không đem lại cải thiện chất lượng dự báo) Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB EMLR 0.36 0.41 0.82 0.64 0.65 0.29 0.19 0.26 0.31 EMES 0.16 0.07 0.14 -0.08 -0.16 -0.04 0.03 0.07 -0.02 EMMV 0.14 0.05 0.14 -0.1 -0.17 -0.03 0.05 0.08 -0.02 BCMA 0.04 0.04 -0.15 -0.15 -0.05 0.02 0.06 -0.04 BCLR 0.36 0.38 0.35 0.41 0.32 0.28 0.21 0.21 0.1 BCES 0.01 0.0 0.02 -0.03 -0.04 -0.01 0.02 0.01 0.02 BCKF 0.0 0.01 -0.05 -0.05 -0.08 -0.12 -0.14 -0.2 -0.11 NGR_ER 0.25 0.36 0.36 0.44 0.37 0.19 0.19 0.14 0.11 NGR_EP 0.2 0.27 0.23 0.28 0.23 0.14 0.17 0.07 0.07 NGR_EM 0.29 0.38 0.37 0.44 0.32 0.23 0.18 0.15 0.08 Raw -0.17 0.2 1.27 1.0 0.72 0.9 0.79 0.51 1.34 Bảng P5.2: Tương tự bảng P5.1 cho hạn dự báo 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB EMLR 0.43 0.61 0.18 0.92 0.37 0.28 0.28 0.25 0.13 EMES 0.12 0.05 0.11 -0.07 -0.18 -0.08 0.02 0.08 -0.01 EMMV 0.08 0.02 0.11 -0.12 -0.2 -0.07 0.02 0.08 -0.01 BCMA 0.0 -0.02 0.03 -0.16 -0.18 -0.07 0.01 0.04 -0.04 BCLR 0.44 0.47 0.41 0.5 0.38 0.32 0.22 0.24 0.12 BCES 0.0 0.01 0.01 -0.04 -0.06 -0.02 0.02 0.01 0.03 BCKF -0.06 -0.06 -0.15 -0.14 -0.17 -0.17 -0.19 -0.2 -0.13 NGR_ER 0.32 0.42 0.38 0.51 0.27 0.21 0.17 0.13 160 0.18 Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EP 0.23 0.29 0.23 0.31 0.2 0.15 0.17 0.07 0.04 NGR_EM 0.31 0.41 0.38 0.48 0.33 0.25 0.19 0.15 0.12 Raw -0.12 0.29 1.33 1.11 0.81 0.91 0.78 0.49 1.31 Bảng P5.3: Tương tự bảng P5.1 cho số MAE Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB EMLR 1.82 2.05 2.32 2.98 1.92 1.57 1.61 1.37 1.28 EMES 1.33 1.44 1.54 1.45 1.32 0.99 0.8 0.83 0.82 EMMV 1.33 1.43 1.53 1.45 1.32 0.99 0.79 0.82 0.81 BCMA 1.18 1.3 1.36 1.4 1.25 0.96 0.79 0.77 0.81 BCLR 1.26 1.35 1.42 1.43 1.23 0.94 0.8 0.8 0.81 BCES 1.22 1.34 1.39 1.45 1.28 0.98 0.79 0.76 0.8 BCKF 1.39 1.55 1.69 1.67 1.48 1.24 1.12 0.97 1.08 NGR_ER 1.67 1.91 2.01 2.08 1.67 1.34 1.15 1.15 1.22 NGR_EP 1.19 1.33 1.4 1.38 1.19 0.95 0.82 0.83 0.87 NGR_EM 1.23 1.34 1.41 1.4 1.2 0.95 0.8 0.79 0.82 Raw 1.91 1.85 1.85 1.61 1.55 1.32 1.13 1.06 1.44 Bảng P5.4: Tương tự bảng P5.1 cho số MAE hạn 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB EMLR 1.34 2.09 2.35 2.62 1.98 1.73 1.46 1.51 1.25 EMES 1.37 1.48 1.59 1.51 1.35 1.05 0.83 0.85 0.84 EMMV 1.36 1.46 1.58 1.5 1.35 1.04 0.83 0.84 0.83 BCMA 1.19 1.32 1.39 1.43 1.25 0.81 0.77 0.82 BCLR 1.31 1.4 1.48 1.5 1.27 0.99 0.82 0.82 0.83 BCES 1.24 1.37 1.43 1.48 1.28 1.02 0.82 0.77 0.83 161 Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCKF 1.45 1.56 1.73 1.73 1.55 1.35 1.19 1.14 NGR_ER 1.74 1.99 2.07 2.15 1.74 1.42 1.22 1.21 1.25 NGR_EP 1.23 1.37 1.45 1.42 1.23 0.98 0.85 0.85 0.91 NGR_EM 1.26 1.37 1.46 1.47 1.24 0.98 0.83 0.81 0.84 Raw 1.91 1.87 1.89 1.69 1.59 1.37 1.13 1.05 1.43 Bảng P5.5: Tương tự bảng P5.1 cho số RMSE Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB EMLR 2.02 2.67 2.91 2.83 2.64 2.03 1.95 1.68 1.72 EMES 1.8 1.9 2.02 1.8 1.69 1.3 1.01 1.09 1.03 EMMV 1.79 1.89 2.02 1.81 1.68 1.29 1.01 1.08 1.03 BCMA 1.61 1.72 1.8 1.74 1.6 1.25 1.01 1.02 BCLR 1.72 1.85 1.92 1.82 1.64 1.27 1.03 1.08 1.02 BCES 1.66 1.79 1.85 1.82 1.66 1.27 1.01 0.99 1.03 BCKF 1.88 2.03 2.2 2.12 1.96 1.62 1.46 1.3 1.39 NGR_ER 2.29 2.61 2.73 2.77 2.27 1.8 1.51 1.54 1.6 NGR_EP 1.67 1.83 1.9 1.77 1.61 1.27 1.06 1.11 1.11 NGR_EM 1.7 1.85 1.91 1.8 1.62 1.27 1.03 1.07 1.04 Raw 2.39 2.4 2.42 2.08 1.99 1.69 1.45 1.34 1.77 Bảng P5.6: Tương tự bảng P5.1 cho số RMSE hạn 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB EMLR 2.09 2.82 2.91 2.93 2.67 2.05 1.89 1.83 1.66 EMES 1.83 1.93 2.09 1.89 1.71 1.37 1.04 1.1 1.05 EMMV 1.8 1.92 2.08 1.89 1.71 1.36 1.04 1.09 1.05 BCMA 1.61 1.74 1.85 1.79 1.61 1.3 1.03 1.01 1.04 162 Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCLR 1.78 1.92 2.02 1.92 1.71 1.34 1.06 1.12 1.05 BCES 1.68 1.81 1.9 1.88 1.67 1.33 1.04 1.01 1.05 BCKF 1.95 2.03 2.24 2.23 2.07 1.74 1.56 1.29 1.45 NGR_ER 2.39 2.69 2.83 2.83 2.31 1.93 1.59 1.62 1.66 NGR_EP 1.7 1.86 2.05 1.83 1.64 1.32 1.09 1.13 1.15 NGR_EM 1.73 1.89 1.98 1.89 1.68 1.33 1.07 1.1 1.07 Raw 2.38 2.42 2.46 2.17 2.03 1.74 1.45 1.34 1.76 Bảng P5.7: Kết tính tốn số đánh giá CRPS dự báo xác suất từ Raw phương pháp EMOS cho yếu tố Tmin với hạn dự báo 48h (các giá trị bôi đậm ngụ ý không đem lại cải thiện chất lượng dự báo) Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 0.89 1.07 1.08 0.95 0.74 0.61 0.57 0.62 BCLR 1.05 1.15 1.2 1.23 1.04 0.8 0.69 0.68 0.7 BCES 0.93 1.04 1.1 1.13 0.99 0.75 0.61 0.56 0.62 BCKF 1.07 1.21 1.31 1.29 1.13 0.92 0.8 0.72 0.78 NGR_ER 1.42 1.63 1.72 1.78 1.42 1.13 0.97 0.97 1.02 NGR_EP 0.91 1.01 1.05 1.02 0.89 0.71 0.61 0.62 0.65 NGR_EM 0.91 1.05 1.03 0.89 0.7 0.58 0.58 0.59 Raw 1.49 1.43 1.48 1.21 1.18 1.02 0.87 0.79 1.12 Bảng P5.8: Tương tự bảng P5.7 cho hạn dự báo 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 0.89 1.08 1.09 0.94 0.76 0.62 0.57 0.63 BCLR 1.08 1.18 1.24 1.28 1.07 0.85 0.71 0.7 0.72 BCES 0.93 1.04 1.12 1.14 0.98 0.77 0.63 0.57 0.64 163 BCKF 1.08 1.19 1.31 1.29 1.15 0.97 0.85 0.74 0.82 NGR_ER 1.5 1.71 1.77 1.85 1.49 1.21 1.03 1.03 1.06 NGR_EP 0.94 1.04 1.1 1.05 0.93 0.74 0.63 0.64 0.68 NGR_EM 0.93 1.03 1.09 1.08 0.92 0.73 0.61 0.6 0.61 Raw 1.46 1.43 1.49 1.26 1.2 1.05 0.87 0.78 1.11 Bảng P5.9: Tương tự bảng P5.7 cho số Độ phủ 90.48% Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 0.71 0.66 0.6 0.59 0.64 0.68 0.67 0.77 0.67 BCLR 0.4 0.38 0.38 0.32 0.36 0.39 0.36 0.39 0.34 BCES 0.7 0.65 0.6 0.58 0.64 0.68 0.67 0.77 0.67 BCKF 0.72 0.65 0.63 0.59 0.66 0.77 0.87 0.91 0.9 NGR_ER 0.4 0.41 0.4 0.39 0.41 0.42 0.44 0.42 0.44 NGR_EP 0.74 0.75 0.76 0.78 0.77 0.78 0.78 0.75 0.75 NGR_EM 0.79 0.79 0.8 0.81 0.81 0.82 0.84 0.83 0.86 Raw 0.54 0.62 0.56 0.69 0.65 0.64 0.68 0.74 0.6 Bảng P5.10: Tương tự bảng P5.7 cho số Độ phủ 90.48% hạn 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 0.75 0.71 0.64 0.65 0.7 0.7 0.67 0.79 0.68 BCLR 0.43 0.41 0.4 0.34 0.39 0.38 0.36 0.41 0.35 BCES 0.75 0.7 0.63 0.64 0.7 0.69 0.68 0.8 0.69 BCKF 0.76 0.7 0.68 0.67 0.71 0.8 0.86 0.92 0.9 NGR_ER 0.39 0.39 0.4 0.38 0.4 0.41 0.42 0.41 0.43 NGR_EP 0.72 0.74 0.75 0.78 0.75 0.77 0.75 0.74 0.75 NGR_EM 0.79 0.79 0.79 0.8 0.81 0.81 0.82 0.82 0.85 Raw 0.59 0.65 0.6 0.73 0.69 0.66 0.68 0.77 0.62 164 Bảng P5.11: Tương tự bảng P5.7 cho số Độ rộng 90.48% Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 3.24 3.32 3.26 3.07 2.51 2.08 2.58 2.17 BCLR 1.89 1.84 1.95 1.67 1.58 1.28 0.98 1.16 0.96 BCES 3.27 3.3 3.21 3.05 2.51 2.08 2.55 2.17 BCKF 4.67 4.53 4.5 4.38 4.34 4.92 5.74 6.1 6.18 NGR_ER 2.1 2.42 2.46 2.57 2.14 1.79 1.63 1.55 1.75 NGR_EP 3.27 3.83 4.08 4.28 3.51 2.92 2.53 2.36 2.58 NGR_EM 3.86 4.28 4.47 4.69 3.96 3.23 2.86 2.74 3.08 Raw 3.84 4.11 3.67 4.34 3.76 3.04 2.69 3.15 3.31 Bảng P5.12: Tương tự bảng P5.7 cho số Độ rộng 90.48% hạn 72h Khu vực nghiên cứu Phương pháp EMOS TB VB ĐB ĐBBB BTB TTB NTB TN NB BCMA 3.63 3.75 3.4 3.76 3.53 2.78 2.2 2.73 2.28 BCLR 2.1 2.07 2.22 1.93 1.81 1.34 1.2 0.99 BCES 3.68 3.76 3.43 3.76 3.55 2.82 2.23 2.76 2.31 BCKF 5.4 5.26 5.35 5.43 5.24 5.76 6.23 6.75 6.71 NGR_ER 2.12 2.42 2.55 2.56 2.15 1.8 1.63 1.55 1.73 NGR_EP 3.27 3.8 4.09 4.25 3.49 2.93 2.52 2.34 2.65 NGR_EM 3.89 4.32 4.58 4.74 3.99 3.26 2.86 2.75 3.1 Raw 4.23 4.57 4.1 4.85 4.19 3.33 2.8 3.33 3.43 165 PHỤ LỤC Kết đánh giá thử nghiệm nghiệp vụ từ 2011-2013 Bảng P6.1: Kết tính toán số đánh giá ME dự báo TBTH cho yếu tố T2m hạn 48h từ phương pháp NGR_EM Raw năm 2011-2013 Năm 2011 2012 2013 Phương pháp Khu vực nghiên cứu TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM -0.12 -0.06 0.32 -0.05 -0.17 -0.25 0.01 -0.03 0.04 Raw -0.98 -0.48 0.13 0.48 0.49 0.58 -0.51 -0.4 0.07 NGR_EM 0.02 0.17 0.13 -0.01 -0.02 0.08 0.1 0.27 Raw -1.5 -0.96 -0.08 0.17 0.01 0.15 -0.45 -0.33 0.06 NGR_EM -0.12 -0.05 0.15 0.07 -0.05 -0.09 0.03 0.03 0.06 Raw -1.35 -0.96 0.12 0.12 0.06 0.26 -0.29 -0.27 0.08 Bảng P6.2: Tương tự bảng P6.1 cho số MAE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 0.88 1.03 1.22 1.03 0.96 0.98 0.8 0.71 0.89 Raw 1.99 1.71 1.37 1.18 1.28 1.27 1.01 1.14 0.88 NGR_EM 0.9 0.92 1.1 1.12 0.91 0.86 0.73 0.73 0.79 Raw 2.35 1.95 1.41 1.28 1.37 1.17 0.95 1.12 0.88 NGR_EM 1.2 1.14 1.28 1.15 0.99 0.95 0.82 0.85 0.79 Raw 2.23 1.92 1.51 1.28 1.41 1.2 0.98 1.08 0.87 Bảng P6.3: Tương tự bảng P6.1 cho số RMSE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 1.16 1.3 1.53 1.29 1.27 1.3 1.04 0.9 1.14 Raw 2.45 2.13 1.74 1.47 1.61 1.62 1.25 1.47 1.11 NGR_EM 1.25 1.21 1.42 1.39 1.2 1.16 0.98 0.98 1.03 Raw 2.9 2.44 1.83 1.61 1.75 1.49 1.21 1.45 1.12 NGR_EM 1.6 1.51 1.67 1.48 1.32 1.25 1.08 1.13 1.03 Raw 2.81 2.44 1.95 1.63 1.79 1.51 1.25 1.4 1.11 166 Bảng P6.4: Kết tính tốn số đánh giá ME dự báo TBTH cho yếu tố Td2m hạn 48h từ phương pháp NGR_EM Raw năm 2011-2013 Năm 2011 2012 2013 Phương pháp Khu vực nghiên cứu TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM -0.17 -0.25 -0.14 -0.32 -0.37 -0.14 -0.04 -0.03 0.22 Raw -1.06 -1.19 -0.57 -0.98 -1.33 -0.91 -0.52 0.16 -0.76 NGR_EM 0.07 0.1 0.14 0.18 0.13 0.1 0.21 Raw -1.49 -1.61 -0.83 -1.32 -1.68 -1.13 -0.7 0.06 -0.97 0.13 0.12 0.15 0.05 0.08 0.13 -1.78 -2.04 -1.17 -1.85 -1.96 -1.1 -0.7 0.15 -1.06 NGR_EM Raw 0.1 -0.04 Bảng P6.5: Tương tự bảng P6.4 cho số MAE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 0.83 1.01 0.93 1.08 0.96 0.93 0.79 0.69 0.69 Raw 1.97 1.41 1.57 1.81 1.62 1.07 1.14 1.23 NGR_EM 0.89 0.99 1.11 1.23 1.01 0.85 0.78 0.67 0.69 Raw 2.31 2.33 1.66 1.93 2.06 1.57 1.12 1.1 1.36 NGR_EM 1.22 1.24 1.3 1.39 1.2 0.95 0.79 0.81 0.74 Raw 2.52 2.65 1.99 2.21 2.23 1.47 1.1 1.08 1.39 Bảng P6.6: Tương tự bảng P6.4 cho số RMSE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 1.11 1.29 1.2 1.32 1.23 1.24 1.04 0.92 0.91 Raw 2.43 2.49 1.89 2.0 2.18 1.97 1.34 1.51 1.63 NGR_EM 1.29 1.33 1.52 1.69 1.46 1.14 1.07 0.89 0.91 Raw 2.83 2.89 2.26 2.46 2.5 1.93 1.43 1.42 1.79 NGR_EM 1.66 1.68 1.77 1.86 1.6 1.25 1.05 1.13 0.98 Raw 3.15 3.37 2.73 2.99 2.84 1.83 1.41 1.44 1.84 167 Bảng P6.7: Kết tính tốn số đánh giá ME dự báo TBTH cho yếu tố Tmax hạn 48h từ phương pháp NGR_EM Raw năm 2011-2013 Năm 2011 2012 2013 Phương pháp Khu vực nghiên cứu TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM -0.38 -0.55 -0.21 -0.51 -0.51 -0.26 -0.03 -0.07 -0.07 Raw -2.69 -1.41 -1.35 0.72 -0.8 -2.19 -2.94 -2.41 -1.34 NGR_EM 0.39 0.17 0.2 0.08 0.35 0.51 0.32 0.28 Raw -3.12 -2.3 -1.67 0.0 -1.16 -2.35 -2.83 -2.04 -1.25 NGR_EM 0.08 0.09 0.11 0.04 0.09 0.16 0.14 0.14 0.14 Raw -2.82 -2.08 -1.27 0.4 -0.8 -2.08 -2.63 -2.1 -1.08 0.25 Bảng P6.8: Tương tự bảng P6.7 cho số MAE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 2.01 2.02 1.88 2.1 2.28 1.72 1.13 1.22 1.03 Raw 3.96 3.38 2.7 2.23 2.68 3.11 3.03 2.71 1.87 1.81 1.73 1.67 1.91 1.58 1.12 1.24 0.98 Raw 4.15 3.7 2.87 2.18 2.63 2.97 2.96 2.42 1.87 NGR_EM 2.24 2.2 2.21 2.18 2.24 1.7 1.11 1.39 1.09 Raw 3.98 3.59 2.91 2.16 2.55 2.69 2.79 2.5 1.76 NGR_EM Bảng P6.9: Tương tự bảng P6.7 cho số RMSE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 2.58 2.56 2.43 2.64 2.92 2.26 1.46 1.53 1.31 Raw 4.69 4.14 3.32 2.88 3.31 3.6 3.43 3.18 2.28 NGR_EM 2.55 2.4 2.28 2.18 2.47 2.19 1.49 1.62 1.34 Raw 4.87 4.43 3.48 2.8 3.28 3.44 3.31 2.91 2.29 NGR_EM 2.88 2.81 2.82 2.82 2.88 2.28 1.51 1.83 1.5 Raw 4.67 4.33 3.52 2.8 3.22 3.13 3.18 2.99 2.18 168 Bảng P6.10: Kết tính tốn số đánh giá ME dự báo TBTH cho yếu tố Tmin hạn 48h từ phương pháp NGR_EM Raw năm 2011-2013 Năm 2011 2012 2013 Phương pháp Khu vực nghiên cứu TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM -0.12 -0.08 0.13 -0.07 -0.14 -0.24 -0.06 -0.04 0.01 Raw -0.38 0.03 0.65 0.89 0.67 0.78 0.18 0.16 1.03 NGR_EM 0.02 0.17 0.13 -0.01 -0.02 0.08 0.1 0.27 Raw -1.5 -0.96 -0.08 0.17 0.01 0.15 -0.45 -0.33 0.06 NGR_EM -0.06 0.09 0.23 0.22 0.05 -0.05 -0.02 -0.01 0.06 Raw -0.33 -0.03 1.01 0.98 0.61 0.64 0.38 0.27 1.09 Bảng P6.11: Tương tự bảng P6.10 cho số MAE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 1.04 1.23 1.35 1.35 1.18 1.07 0.9 0.71 0.9 Raw 1.96 1.86 1.6 1.56 1.58 1.49 1.15 1.27 NGR_EM 1.11 1.12 1.18 1.26 1.18 0.94 0.84 0.72 0.86 Raw 2.21 1.99 1.71 1.61 1.68 1.32 0.96 1.18 1.29 NGR_EM 1.39 1.35 1.5 1.51 1.29 1.03 0.88 0.9 0.86 Raw 2.12 1.94 1.92 1.72 1.67 1.38 1.05 1.2 1.3 Bảng P6.12: Tương tự bảng P6.10 cho số RMSE Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 1.38 1.57 1.76 1.74 1.55 1.42 1.15 0.92 1.16 Raw 2.43 2.33 2.03 1.99 1.9 1.28 1.43 1.59 NGR_EM 1.55 1.57 1.65 1.61 1.6 1.37 1.07 0.95 1.09 Raw 2.73 2.5 2.28 2.16 2.16 1.7 1.24 1.49 1.62 NGR_EM 1.82 1.78 1.95 1.73 1.4 1.15 1.2 1.11 Raw 2.66 2.47 2.47 2.24 2.15 1.78 1.36 1.51 1.63 169 Bảng P6.13: Kết tính tốn số CRPS dự báo xác suất cho yếu tố T2m hạn 48h từ phương pháp NGR_EM Raw năm 2011-2013 Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 0.65 0.75 0.88 0.75 0.7 0.71 0.59 0.51 0.66 Raw 1.62 1.35 1.07 0.86 0.97 0.98 0.77 0.9 0.66 NGR_EM 0.66 0.68 0.81 0.81 0.67 0.64 0.54 0.53 0.58 Raw 1.92 1.55 1.09 0.97 1.04 0.9 0.73 0.89 0.67 NGR_EM 0.89 0.85 0.94 0.84 0.74 0.7 0.61 0.62 0.58 Raw 1.83 1.54 1.17 0.96 1.07 0.92 0.78 0.84 0.66 Bảng P6.14: Tương tự bảng P6.13 cho yếu tố Td2m Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 0.6 0.72 0.67 0.77 0.68 0.67 0.57 0.5 0.51 Raw 1.55 1.55 1.05 1.16 1.37 1.21 0.78 0.87 0.93 NGR_EM 0.67 0.73 0.82 0.92 0.77 0.63 0.58 0.49 0.51 Raw 1.84 1.83 1.27 1.45 1.59 1.17 0.83 0.85 1.01 NGR_EM 0.91 0.93 0.97 1.04 0.9 0.7 0.58 0.6 0.55 2.11 1.52 1.71 1.69 1.1 0.83 0.81 1.04 Raw Bảng P6.15: Tương tự bảng P6.13 cho yếu tố Tmax Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 1.45 1.47 1.38 1.53 1.66 1.27 0.83 0.88 0.76 Raw 3.47 2.83 2.19 1.74 2.14 2.57 2.62 2.16 1.47 NGR_EM 1.49 1.36 1.28 1.25 1.41 1.2 0.83 0.91 0.73 Raw 3.67 3.09 2.33 1.71 2.11 2.43 2.54 1.9 1.47 NGR_EM 1.65 1.63 1.62 1.61 1.65 1.26 0.84 1.04 0.84 Raw 3.51 3.05 2.42 1.75 2.09 2.16 2.46 1.98 1.39 170 Bảng P6.16: Tương tự bảng P6.13 cho yếu tố Tmin Năm 2011 2012 2013 Khu vực nghiên cứu Phương pháp TB VB ĐB ĐBBB BTB TTB NTB TN NB NGR_EM 0.75 0.88 0.97 0.98 0.85 0.77 0.66 0.51 0.66 Raw 1.6 1.49 1.28 1.21 1.24 1.19 0.8 0.9 1.03 NGR_EM 0.83 0.83 0.87 0.91 0.88 0.71 0.61 0.52 0.62 Raw 1.8 1.6 1.37 1.27 1.32 1.04 0.76 0.93 1.05 NGR_EM 1.02 1.11 1.11 0.95 0.77 0.65 0.66 0.63 Raw 1.75 1.57 1.56 1.33 1.31 1.08 0.84 0.94 1.04 171 ... thực tế 33 1.3.2 Trong nước Tại Việt Nam, NWP cịn giai đoạn bước đầu tiếp thu cơng nghệ nghiên cứu ứng dụng Do đó, DBTH giai đoạn bước đầu tìm hiểu thử nghiệm Các nghiên cứu ứng dụng DBTH tập trung... phục ngun nhân mơ hình NWP, rõ ràng cần phải đầu tư nghiên cứu cải tiến mơ hình từ động lực, vật lý phương pháp số Để cải tiến trường ban đầu, cần phải nghiên cứu ứng dụng tốn đồng hóa số liệu,... giản hàm phân bố xác suất rút từ phân bố (2.16) Giả định k ngày số liệu khứ độc lập, hàm cần tìm tích k hàm phân bố xác suất Cực đại hàm cho ta tập hệ số cần tìm ứng với tập hệ số nói trên, xác