Bảng 3.14. Sai số khoảng cách dự báo tổ hợp 5 thành phần CMA (số liệu độc lập) (số liệu độc lập)
Hạn dự báo
Sai số khoảng cách (km)
PP trung bình đg PP sử dụng phương sai của sai số PP siêu tổ hơ ̣p
24h 169.6 172.19 129.56
48h 349.5 385.15 324.91
72h 442.8 446.25 448.48
96h 533.6 540.21 502.73
120h 599.1 587.91 535.51
Hình 3.6. Biểu đồ sai số khoảng cách (km) dự báo tổ hợp 5 tp CMA ( số liệu độc lập) ( số liệu độc lập)
Nhận xét: Từ bảng so sánh (Bảng 3.14 ) và biểu đồ (Hình 3.6) kết quả sai
số khoảng cách trên mẫu độc lập tương tự kết quả trên mẫu phụ thuộc giữa ba phương pháp tìm hệ số của tổ hợp ta thấy đối với mẫu số liệu nào thì phương pháp siêu tổ hơ ̣p cũng cho sai số khoảng cách (ở tất cả các hạn dự báo) thấp nhất. Kết quả này chứng tỏ rằng khi tổ hợp bằng phương pháp siêu tở hơ ̣p ta có thể tìm được các hệ số phù hợp cho các thành phần tham gia tổ hợp dựa vào tương quan của các nhân tố dự báo với nhau và tương quan của các nhân tố dự báo với yếu tố dự báo.
c. Kết quả thử nghiệm với 1 thành phần của CMA và thành phần của ECMWF
Sai số khoảng cách đánh giá trên chuỗi số liệu phụ thuộc:
Bảng 3.15. Sai số khoảng cách dự báo tổ hợp 1 thành phần của CMA và thành phần ECMWF (số liệu phụ thuộc)
Hạn dự báo
Sai số khoảng cách (km)
PP trung bình đg PP sử dụng phương sai của sai số PP siêu tổ hợp
24h 106.62 91.49 77.53
48h 171.02 156.2 135.11
72h 226.41 194.98 188.9
96h 305.2 289.69 269.32
120h 409.95 396.27 384.66
Hình 3.7. Biểu đờ sai số khoảng cách (km) dự báo tổ hợp 1 tp CMA và thành phần ECMWF ( số liệu phụ thuộc)
Nhận xét: Phân tích sai số khoảng cách khi tổ hợp một thành phần dự báo từ
hai trung tâm trên mẫu số liệu phụ thuộc (Bảng 3.15, Hình 3.7), kết quả sai số khơng cịn cao như các trường hợp tổ hợp trước: ở hạn 24h sai số cho 3 phương pháp tổ hợp là từ 77km đến hơn 100km; đến hạn 48h từ 135km đến 171km; đến hạn 72h từ 188km đến 226km; đến hạn 96h từ 269km đến 305km; 120h từ 384km đến hơn 400km. Đồng thời sai số của cả 3 phương pháp tổ hợp cũng không chênh lệch nhau nhiều và sắp xếp theo thứ tự tăng dần về sai số thì phương pháp siêu tở hợp ln đứng đầu cịn phương pháp trung bình đơn giản xếp cuối cùng, ở giữa là phương pháp dùng phương sai của sai số.
Sai số khoảng cách đánh giá trên chuỗi số liệu độc lập:
Bảng 3.16. Sai số khoảng cách dự báo tổ hợp 1 thành phần của CMA và thành phần ECMWF (số liệu độc lập)
Hạn dự báo
Sai số khoảng cách (km)
PP trung bình đg PP sử dụng phương
sai của sai số PP siêu tổ hợp
24h 135.14 123.42 121.74
48h 219.52 209.46 207.48
72h 322.09 345.02 319.68
96h 398.14 412.67 453.63
120h 511.67 517.43 485.9
Hình 3.8. Biểu đờ sai số khoảng cách (km) dự báo tổ hợp 1 tp CMA và thành phần ECMWF ( số liệu độc lập)
Nhận xét: Phân tích sai số khoảng cách khi tổ hợp một thành phần dự báo
từ hai trung tâm trên mẫu số liệu độc lập (Bảng 3.16, Hình 3.8) kết quả cho thấy sai số cao hơn so với mẫu phụ thuộc. Cũng như ở Mẫu phụ thuộc, Phương pháp siêu tổ hơ ̣p cho sai số khoảng cách thấp nhất ở các hạn dự báo 24h, 48h, 72h và 120h tuy nhiên ở hạn dự báo 96h sai số khoảng các của tổ hợp dự báo theo phương pháp này lại lớn nhất và sai số của phương pháp trung bình đơn giản lại cho kết quả sai số nhỏ nhất.
d. Kết quả thử nghiệm với 15 thành phần của CMA và thành phần của ECMWF
Sai số khoảng cách đánh giá trên chuỗi số liệu phụ thuộc:
Bảng 3.17. Sai số khoảng cách dự báo tổ hợp 15 thành phần của CMA và thành phần ECMWF (số liệu phụ thuộc)
Hạn dự báo
Sai số khoảng cách (km)
PP trung bình đg PP sử dụng phương sai
của sai số PP siêu tổ hợp
24h 112.17 125.32 60.75
48h 163.27 163.69 101.73
72h 242.29 238.35 165.22
96h 345.55 321.23 219.68
Hình 3.9. Biểu đồ sai số khoảng cách (km) dự báo tổ hợp 15 tp CMA và thành phần ECMWF ( số liệu phụ thuộc)
Sai số khoảng cách đánh giá trên chuỗi số liệu độc lập:
Bảng 3.18. Sai số khoảng cách dự báo tổ hợp 15 thành phần của CMA và thành phần ECMWF (số liệu độc lập)
Hạn dự báo
Sai số khoảng cách (km)
PP trung bình đg
PP sử dụng phương sai của sai
số
PP siêu tổ hợp
24h 150.21 143.38 137.16
48h 285.13 225.52 285.56
72h 389.88 371.84 369.21
96h 565.02 521.77 498.63
Hình 3.10. Biểu đồ sai số khoảng cách (km) dự báo tổ hợp 15 tp CMA và thành phần ECMWF ( số liệu độc lập)
Nhận xét: Phân tích các kết quả sai số khoảng cách khi tổ hợp 15 thành
phần dự báo của CMA và thành phần của ECMWF trên mẫu số liệu phụ thuộc (Hình 3.9) cho thấy sai số khoảng cách của phương pháp siêu tổ hơ ̣p nhỏ đi rõ rệt, ở hạn 24h chỉ có hơn 60km; 48h: hơn 100km; 72h: hơn 160km; 96h: hơn 210km; 120h: chỉ có hơn 230km. các sai số này rất nhỏ so với hai phương pháp còn lại. tuy nhiên khi đánh giá trên tập mẫu độc lập (Hình 3.10) thì sai số khoảng cách của phương pháp siêu tổ hơ ̣p lại tăng lên rất cao gần bằng với sai số của phương pháp dùng phương sai của sai số ở các hạn dự báo, và ở hạn 120h lên đến 650km hơn hẳn 420 km so với sai số của hạn dự báo này ở mẫu phụ thuộc. Điều này chứng tỏ phương trình dự báo tổ hợp 15 thành phần của CMA và 1 thành phần của ECMWF chưa mang lại kết quả tốt ở mẫu độc lập, đã có các kết quả cho sai số tương đối cao.
e. Kết quả thử nghiệm với 5 thành phần của CMA và thành phần của ECMWF.
Bảng 3.19. Sai số khoảng cách dự báo tổ hợp 6 thành phần của 2 trung tâm ( Số liệu phụ thuộc)
Hạn dự báo
Sai số khoảng cách (km)
PP trung bình đg PP sử dụng phương
sai của sai số PP siêu tổ hợp
24h 138.29 107.31 71.81
48h 207.22 160.76 130.49
72h 286.91 221.95 187
96h 421.04 306.69 263.52
120h 451.8 385.52 312.58
Hình 3.11. Biểu đờ sai số khoảng cách (km) dự báo tổ hợp 6 thành phần của 2 trung tâm (số liệu phụ thuộc)
Nhận xét: Đồ thị (Hình 3.11) trên đây thể hiện rất rõ sự chênh lệch về kết
quả sai số dự báo khi tổ hợp các kết quả dự báo bằng ba phương pháp khác nhau. Cụ thể, trong khi sai số khoảng cách của phương pháp siêu tổ hợp ở hạn 72h chỉ có khoảng 187km thì phương pháp dùng phương sai của sai số lên đến 221km, cịn phương pháp lấy trung bình đơn giản lên đến hơn 286km, mỗi phương pháp chênh lệch nhau khoảng 50km; 120h cũng chênh lệch nhau khoảng 60km. Như vậy
phương pháp siêu tổ hơ ̣p là phương pháp tối ưu luôn cho sai số nhỏ nhất như vậy đây là phương pháp nên lựa chọn.
Các giá trị sai số: ở hạn 24h là 71.81km; 48h- 130.49km; 72h- 187km; 6h- 283.52km; 120h- 312.58km đã chỉ ra rằng khi tổ hợp các kết quả dụ báo của cả hai trung tâm bằng phương pháp hồi quy chúng ta sẽ có được một dự báo khá tốt, tương đương với kết quả của dự báo bão trước đây là sai số khoảng 350km cho hạn dự báo 120h.
Sai số khoảng cách đánh giá trên chuỗi số liệu độc lập:
Bảng 3.20. Sai số khoảng cách dự báo tổ hợp 6 thành phần của 2 trung tâm (Số liệu độc lập)
Hạn dự báo
Sai số khoảng cách (km)
PP trung bình đg PP sử dụng phương
sai của sai số PP siêu tổ hơ ̣p
24h 154.94 126.09 110.6
48h 299.73 275.16 159.37
72h 393.25 352.33 238.5
96h 462.07 452.91 343.23
120h 581.05 590.56 362.01
Nhận xét: Đồ thị: Hình 3.12, kết quả này cũng cho thấy tương tự với kết
quả của tập mẫu phụ thuộc.
So sánh kết quả giữa hai tâ ̣p mẫu phu ̣ thuô ̣c và đô ̣c lâ ̣p ta thấy , với phương pháp siêu tổ hơ ̣p sai số trên hai tập mẫu chênh lệch nhau không nhiều lắm, cao nhất là ở hạn 96h khoảng 80km, thấp nhất ở hạn 48h khoảng 30km; ở hạn 120h cũng chênh lệch ít vào khoảng 50km. Như vậy chất lượng dự báo của số liệu độc lập và số liệu phụ thuộc là tương đối gần nhau , và các kết quả dự báo đều khá tốt.
3.4. So sánh các kết quả dự báo có sử dụng phƣơng pháp hồi quy tuyến tính:
Hình 3.13. Đồ thị sai số khoảng cách của các trường hợp dự báo trên tập mẫu phụ thuộc
Hình 3.14. Đồ thị sai số khoảng cách của các trường hợp dự báo trên tập mẫu độc lập
Tổng hợp các kết quả dự báo ở trên ta có được một số nhận xét sau đây:
Trên tập mẫu phụ thuộc (Hình 3.13)
- Khi ta sử dụng kết quả dự báo của từng trung tâm để xây dựng phương
trình dự báo thì kết quả dự báo ln cho sai số lớn hơn so với khi ta tổ hợp các kết quả dự báo lại để xây dựng phương trình.
- Sai số dự báo của CMA ở những hạn dự báo ngắn hạn 24h, 48h thấp hơn
so với sai số dự báo của ECMWF, còn đối với những dự báo hạn 96h,120h thì ngược lại. Điều đó chứng tỏ, dự báo của CMA chỉ có chất lượng tốt ở những hạn dự báo ngắn như 24h, 48h; cịn dự báo của ECMWF lại có chất lượng tốt hơn ở hạn dự báo 96h và 120h.
- Tổ hợp các kết quả của trung tâm CMA kết quả sai số khoảng cách có giảm đi so với khi khơng tổ hợp, ví dụ ở hạn 120h: CMA cho sai số khoảng hơn 500km, CMA(5) cho sai số khoảng gần 450km, ở CMA(15) sai số chỉ còn gần 400km.
- Tổ hợp các kết quả dự báo của cả hai trung tâm CMA- ECMWF đã cho kết quả khả quan hơn, cụ thể từng trường hợp thử nghiệm tổ hợp được thể hiện rất rõ trên hình vẽ: Tổ hợp CMA(15)- ECMWF cho kết quả sai số thấp nhất; CMA(1)- ECMWF và CMA(5)- ECMWF cho kết quả sai số gần bằng nhau ở các hạn dự báo 24h, 48h, 72h, 96h; cịn với hạn 120h thì sai số khoảng cách của CMA(5)-ECMWF nhỏ hơn khoảng gần 100km. Sai số khi tổ hợp của CMA(5)- ECMWF và CMA(15)- ECMWF sai số hạn 120h chỉ khoảng từ 200km đến 300km, trong khi đó sai số khi dự báo thử nghiệm từ các mơ hình như từ kết quả của nhóm nghiên cứu trường ĐH Khoa học tự nhiên hay từ một số nghiên cứu của nước ngoài ở hạn này là khoảng 350km, như vậy kết quả này là tương đối khả quan.
Trên tập mẫu độc lập (Hình 3.14 )
- Cũng cho kết quả tương tự như số liệu phụ thuộc tuy nhiên sai số của dự
báo CMA tăng lên khá nhiều còn sai số của dự báo ECMWF tăng lên không nhiều lắm, điều này chứng tỏ kết quả dự báo của ECMWF có tính ổn định hơn các dự báo của CMA.
- Ở tập mẫu độc lập có sự thay đổi rất lớn về kết quả sai số khoảng cách so với mẫu phụ thuộc, sai số cao nhất lại thuộc về dự báo tổ hợp CMA(15), sau đó đến CMA, đến CMA(5), đến CMA(15)-ECMWF, đến ECMWF, đến CMA(1)-ECMWF và cuối cùng sai số nhỏ nhất lại là của CMA(5)-ECMWF. dự báo tổ hợp CMA(15)- ECMWF khơng cịn cho kết quả sai số thấp nhất nhất nữa điều này có thể cho ta thấy khi tổ hợp nhiều kết quả dự báo của CMA không phải lúc nào kết quả dự báo cũng cho sai số thấp, như vậy dự báo này là không ổn định.
- Ở tập mẫu độc lập cũng giống với mẫu phụ thuộc, dự báo tổ hợp CMA(5)-ECMWF tất cả các dự báo của cả hai trung tâm vẫn cho kết quả tốt nhất và sai số ở 120 h là khoảng 360km, tương đương với kết quả dự báo đánh gía sai số trung bình của khu vực Tây Bắc Thái Bình Dương là khoảng 300-350km (Trung Tâm cảnh báo bão Mỹ) và trong dự báo thử nghiệm của TH.S Hoàng Thị Thủy (2013) trong luận văn thạc sĩ: “Thử nghiệm dự báo tổ hợp quỹ đạo và cường độ bão
hạn 5 ngày trên khu vực Biển Đông bằng WRF sử dụng sản phẩm tổ hợp toàn cầu”
kết quả sai số dự báo ở 120h cũng khoảng 360km
Từ các nhận xét trên cũng có thể kết luận rằng phương trình dự báo tối ưu nhất là phương trình tổ hợp CMA(5)- ECMWF với phương pháp tổ hợp là thực hiện hồi quy tuyến tính (HQTT) đa biến có chọn lọc hay cịn gọi là phương pháp siêu tổ hợp.
3.5. Đánh giá kết quả dự báo CMA(5)- ECMWF trên một số cơn bão điển hình
Trong tập mẫu số liệu độc lập dùng để đánh giá kết quả dự báo gồm có các cơn bão mạnh, cơn bão yếu, cơn bão hoạt động đầu mùa, cơn bão hoạt động cuối mùa. Để đánh giá sâu hơn về chất lượng các phương trình dự báo tơi lựa chọn hai cơn bão điển hình cơn bão Nari (10- 15/10/2013) cơn bão cuối mùa , cơn bão Rammasun (13- 20/7/2014) cơn bão đầu mùa để tiến hành phân tích và đánh giá.
Khảo sát bão Nari:
Diễn biến quỹ đạo thực của cơn bão Nari (10- 15/10/2013):
Nari được đánh giá là cơn bão có diễn biến phức tạp cường độ rất mạnh, hình thành ngày (8/10/2013) từ một xốy thuận nhiệt đới trên đảo Lu- dơng (Philipin), ngày 10/10 cơn bão này di chuyển nhanh về phía Philipin; sau đó vượt qua quần đảo Lu-Dông đi vào biển Đông với hướng di chuyển khá ổn định.
Hồi 01 giờ ngày 11/10, vị trí tâm bão ở vào khoảng 15,3 độ Vĩ Bắc; 125,4 độ Kinh Đông, cách đảo Lu-Dông (Philippin) khoảng 400km về phía Đơng
Hồi 01 giờ ngày 12/10, vị trí tâm bão ở vào khoảng 15,6 độ Vĩ Bắc; 120,5 độ Kinh Đơng, trên khu vực phía Tây Nam đảo Lu-Dông (Philippin).
Hồi 04 giờ ngày 14/10, vị trí tâm bão ở vào khoảng 15,6 độ Vĩ Bắc; 111,8 độ Kinh Đơng, trên khu vực phía Nam quần đảo Hồng Sa. Cách bờ biển các tỉnh Quảng Trị - Quảng Ngãi khoảng 330km về phía Đơng.
Từ 2 giờ đến 7 giờ sáng ngày 15/10, bão Nari đổ bộ vào TP.Đà Nẵng và tỉnh Quảng Nam, trở thành trận cuồng phong lớn, kéo dài nhất từ trước tới nay với sức gió giật cấp 11, cấp 12, giật cấp 13.
Dự báo quỹ đạo bão Nari:
Thời điểm bắt đầu thực hiện sự báo là 00Z ngày 10/10/2013 bão Nari có vị trí : 15.0 N – 127.80 E. Quá trình dự báo (Hình 3.16), ở giai đoạn đầu khi chưa vượt qua đảo Ludong- Philipin cơn bão di chuyển chậm hơn so với dự báo thực. Ngày 11/10/2013 tâm bão dự báo cách tâm bão thực 63,2 km (về phía bắc ); đến ngày 12/10/2013 tâm bão dự báo cách tâm bão thực khoảng 180km (về phía bắc ); Giai đoạn sau, khi đi vào khu vực biển Đông cơn bão di chuyển nhanh hơn so với dự báo thực. Ngày 13/10/2013 tâm bão dự báo cách tâm bão thực 162km cũng về phía bắc; tuy nhiên đến ngày 14/10/2013 tâm bão dự báo cách tâm bão thực khá xa - sai số khoảng cách lên đến 420km lệch hẳn vào đảo Hải Nam Trung Quốc. Ngày 15/10/2013 bão dự báo cũng đổ bộ vào đất liền nhưng sớm hơn so với bão thực khoảng cách giữa tâm dự báo và tâm thực là khoảng 138km . Tóm lại , khi khảo sát phương trình dự báo trên cơn bão Nari các hạn dự báo: 24h, 28h, 72h, 120h sai số khoảng cách tương đối nhỏ tuy nhiên ở hạn dự báo 96h sai số lại tăng lên cao điều này có thể do một số dự báo trong thành phần tổ hợp của cơn bão này cho kết quả không tốt. Như vâ ̣y, nếu xét cả quá trình dự báo đồng thời với mu ̣c đích dự báo ở đây là dự báo bão ha ̣n 5 ngày thì kết quả này có thể cho là tương đối khả quan.
Khảo sát bão Rammasun:
Diễn biến quỹ đạo thực của cơn bão Rammasun (13- 20/7/2014):
Rammasun- Thần sấm là một cơn bão rất mạnh, di chuyển nhanh và có diễn biến phức tạp, bắt đầu từ Đông Nam Philippines từ ngày 13/7 và đã vượt qua quần đảo Hoàng Sa, đảo Hải Nam, vịnh Bắc Bộ và tiến vào Quảng Ninh- Việt Nam. Cụ