Tổng quan nghiên cứu

Việt Nam sở hữu đường bờ biển dài hơn 3.200 km tiếp giáp trực tiếp với Biển Đông – một trong những ổ bão nhiệt đới hoạt động mạnh và phức tạp nhất trên thế giới. Thống kê khí tượng nhiều năm ghi nhận trung bình mỗi năm nước ta chịu ảnh hưởng của hơn 7 cơn bão và áp thấp nhiệt đới, trong đó 87% tập trung vào 5 tháng cao điểm từ tháng 7 đến tháng 11. Việc xác định chính xác cường độ và vị trí tâm bão theo thời gian thực đóng vai trò sống còn đối với công tác phòng chống thiên tai, giảm thiểu thiệt hại về người và tài sản.

Tuy nhiên, mạng lưới trạm quan trắc khí tượng trên biển còn rất thưa thớt, khiến việc xác định cường độ thực tế của bão ngoài khơi phụ thuộc chủ yếu vào dữ liệu ảnh mây vệ tinh. Thực tế nghiệp vụ cho thấy sự phân tích giữa các trung tâm dự báo quốc tế thường có sự chênh lệch đáng kể. Điển hình như trong cơn bão số 1 (Kujira) năm 2015, chênh lệch đánh giá cường độ thời gian thực giữa Hoa Kỳ (cấp 7) và Hàn Quốc (xấp xỉ cấp 10) lên tới gần 3 cấp bão, gây khó khăn lớn cho việc ra quyết định ứng phó.

Nhằm giải quyết bài toán khách quan hóa và tự động hóa công tác dự báo, luận văn tập trung đánh giá hiệu quả của phương pháp Dvorak cải tiến (Advanced Dvorak Technique – ADT phiên bản V8) trong việc xác định vị trí và cường độ bão từ ảnh vệ tinh địa tĩnh cho khu vực Biển Đông. Nghiên cứu được thực hiện trên tập mẫu gồm 41 cơn bão hoạt động trong 6 năm (2010 - 2015) với 632 chu kỳ dự báo chuẩn 6 giờ. Kết quả nghiên cứu cung cấp cơ sở khoa học định lượng chuẩn xác nhằm nâng cao tần suất phát tin dự báo bão lên 1 giờ/lần và tối ưu hóa hệ thống cảnh báo sớm thiên tai tại Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng kỹ thuật Dvorak truyền thống và chuỗi phát triển khách quan hóa qua các giai đoạn: phương pháp Dvorak cổ điển (DT, 1972/1984), phương pháp Dvorak khách quan (ODT, 1998), phương pháp Dvorak khách quan cải tiến (AODT, 2007) và phương pháp Dvorak cải tiến (ADT, 2015). Trong đó, các khái niệm cốt lõi bao gồm:

  • Chỉ số T-number: Thước đo định lượng đặc trưng hình thái mây bão.
  • Cường độ hiện tại (CI): Chỉ số phản ánh trạng thái cường độ thực tế của bão sau khi áp dụng các ràng buộc phát triển và suy yếu.
  • Áp suất mực nước biển thấp nhất tại tâm (MSLP - Pmin) và tốc độ gió cực đại duy trì (Vmax).
  • Bán kính gió mạnh nhất (RMW) dựa trên phân tích phổ nhiệt độ đỉnh mây.

Hệ thống ADT tích hợp thuật toán phân tích chuỗi Fourier nhanh (FFT) để xác định tính đối xứng đối lưu, kết hợp phương trình hồi quy phi tuyến tính đa biến cho hai dạng mây đặc trưng là khối mây dày đặc ở trung tâm (CDO) và dạng mắt (EYE). Phân loại mẫu mây trong ADT bao gồm 6 dạng hình thái chính: dạng lệch tâm (SHEAR), dạng băng cuốn (CRVBND), dạng tâm nhúng đĩa mây (EMBC), khối mây đậm đặc có nhiệt độ đồng đều (UNIFRM), khối mây đậm đặc biến đổi (IRRCDO) và dạng mắt bão (EYE).

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu đa phổ từ vệ tinh địa tĩnh MTSAT-2 của Cơ quan Khí tượng Nhật Bản (JMA) giai đoạn 2010 - 2015 với tần suất quét 30 phút/ảnh, bao gồm 5 kênh phổ: kênh thị phổ VIS (độ phân giải không gian 1 km) và 4 kênh hồng ngoại IR1, IR2, WV, IR4 (độ phân giải không gian 4 km). Dữ liệu đối chứng chuẩn là bộ số liệu bão tốt nhất (Best Track) của Trung tâm Dự báo khí tượng thủy văn Quốc gia và bản tin dự báo từ các trung tâm quốc tế (JMA, JTWC).

Cỡ mẫu nghiên cứu gồm 41 cơn bão nhiệt đới trên Biển Đông với tổng cộng 632 chu kỳ quan trắc cách nhau 6 giờ. Toàn bộ các trường hợp tâm bão đã đổ bộ vào đất liền đều được loại bỏ theo đúng nguyên lý hoạt động của thuật toán ADT trên biển. Phương pháp phân tích thống kê toán học sử dụng 4 chỉ số định lượng then chốt:

  • Hệ số tương quan tuyến tính (Cor): Đo lường mức độ đồng biến giữa giá trị ước tính ADT và số liệu Best Track.
  • Sai số hệ thống (BIAS): Đánh giá khuynh hướng ước lượng cao hơn hoặc thấp hơn trung bình thực tế.
  • Sai số tuyệt đối trung bình (MAE): Phản ánh biên độ sai số trung bình tuyệt đối của phép đo.
  • Căn bậc hai của sai số bình phương trung bình (RMSE): Đánh giá mức độ phân tán và trọng số của các sai số lớn.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình kiểm nghiệm 632 thời điểm phân tích của 41 cơn bão trên Biển Đông giai đoạn 2010 - 2015 đem lại những phát hiện quan trọng:

Thứ nhất, sai số xác định vị trí tâm bão tỷ lệ nghịch với độ hoàn thiện cấu trúc mây và cấp bão. Sai số vị trí trung bình đạt mức nhỏ nhất là 28 km đối với các cơn bão có dạng mắt (EYE - 108 trường hợp), trong khi dạng lệch tâm (SHEAR - 261 trường hợp) có sai số vị trí lớn nhất lên tới 44 km (độ lệch giảm 36,4%). Khi phân loại theo cấp gió Beaufort, các cơn bão rất mạnh (cấp 12 trở lên với 238 trường hợp) có sai số vị trí trung bình chỉ 33 km, thấp hơn rõ rệt so với bão cấp 8 - 9 (39 km với 272 trường hợp).

Thứ hai, phương pháp ADT có xu hướng ước tính tốc độ gió cực đại (Vmax) cao hơn so với số liệu Best Track Việt Nam. Mức độ chênh lệch tăng dần theo cường độ bão: dạng băng cuốn (CRVBND - 40 trường hợp) cho kết quả sát thực tế nhất với độ chênh Vmax chỉ +4 kts (sai số tương đối 9,6%, MAE 6,4 kts); ngược lại, dạng mắt (EYE) có sai số Vmax lớn nhất lên tới +24 kts (sai số 27,6%, giá trị trung bình Best Track 87 kts so với ADT 111 kts).

Thứ ba, ước tính trị số khí áp thấp nhất (Pmin) của ADT có độ ổn định và độ tin cậy rất cao. Độ chênh lệch Pmin trung bình giữa ADT và Best Track trên toàn bộ các dạng mây bão đều duy trì dưới 4 mb, với sai số phần trăm tương đối dao động từ 0,26% đến 0,39%. Ngay cả ở cấp bão rất mạnh (cấp 12 trở lên), mức chênh lệch áp suất trung bình chỉ là 3,7 mb (Pmin Best Track 953 mb so với ADT 956,7 mb).

Thứ tư, hệ số tương quan tuyến tính giữa ADT và Best Track đạt mức rất cao ở các cơn bão mạnh có cấu trúc mây phát triển đầy đủ. Cụ thể, hệ số tương quan đạt từ 0,82 đến 0,84 đối với các cơn bão dạng mắt (EYE) và tâm nhúng đĩa mây (EMBC). Ngược lại, ở dạng mây lệch tâm (SHEAR) thuộc giai đoạn mới hình thành hoặc suy yếu, hệ số tương quan Vmax giảm mạnh xuống mức xấp xỉ -0,01.

Thảo luận kết quả

Khi biểu diễn dữ liệu qua các đồ thị phân tán (scatter plots) và đường hồi quy tuyến tính, tương quan giữa ADT và Best Track thể hiện sự phân hóa rõ nét. Đường hồi quy thực nghiệm của tốc độ gió cực đại Vmax và áp suất Pmin đều nằm cao hơn đường phân giác chuẩn, khẳng định khuynh hướng ước tính thiên mạnh của thuật toán ADT khi áp dụng trực tiếp tại vùng biển Tây Bắc Thái Bình Dương và Biển Đông.

Nguyên nhân cốt lõi của sự sai lệch này bắt nguồn từ việc bảng quy đổi chuẩn T-number sang Vmax và Pmin trong ADT gốc được xây dựng dựa trên cơ sở dữ liệu khí hậu bão mở của Đại Tây Dương. Biển Đông là một thủy vực nửa kín, chịu sự chi phối phức tạp của địa hình bán đảo Đông Dương, gió mùa Đông Bắc lạnh khô và gió mùa Tây Nam nóng ẩm, làm biến dạng trường nhiệt độ đỉnh đối lưu và cấu trúc áp suất bề mặt.

Đánh giá chi tiết trên hai cơn bão mạnh điển hình gồm bão số 6 (Megi) năm 2010 và bão số 2 (Rammasun) năm 2014 cho thấy: trong giai đoạn bão đạt cấp 14 - 15 ngoài khơi, sai số vị trí của ADT ổn định dưới 40 km. Mặc dù phương pháp Dvorak thủ công kết hợp kinh nghiệm dự báo viên cho chỉ số MAE thấp hơn ADT khoảng 1 - 2 kts trong các tình huống cụ thể, nhưng ADT vượt trội hoàn toàn về tính khách quan, tốc độ xử lý tức thời và loại bỏ triệt để độ lệch chủ quan giữa các kíp trực.

Đề xuất và khuyến nghị

Nhằm tối ưu hóa năng lực ứng dụng công nghệ ADT vào nghiệp vụ dự báo xoáy thuận nhiệt đới tại Việt Nam, luận văn đề xuất 4 nhóm giải pháp trọng tâm:

  1. Hiệu chỉnh bộ thông số hồi quy địa phương hóa: Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu chủ trì phối hợp cùng Trung tâm Dự báo khí tượng thủy văn Quốc gia tiến hành xây dựng lại bảng quy đổi thực nghiệm giữa chỉ số T-number/CI sang tốc độ gió cực đại (Vmax) và khí áp trung tâm (Pmin) đặc thù cho vùng Biển Đông trong thời gian 12 tháng, đặt mục tiêu giảm sai số Vmax của bão cấp 12 trở lên từ 27,6% xuống dưới 12%.

  2. Tích hợp module ADT tự động vào quy trình trực ca: Trung tâm Công nghệ Khí tượng Thủy văn hoàn thiện cài đặt gói phần mềm ADT phiên bản mới trên hệ thống máy chủ Linux hiệu năng cao, thiết lập tiến trình tự động quét dữ liệu ảnh vệ tinh Himawari với chu kỳ 30 - 60 phút/lần ngay trong quý 3 năm tới, phục vụ phát hành bản tin bão nhanh với độ trễ xử lý dưới 5 phút.

  3. Thiết lập chuỗi kết hợp đa dữ liệu quan trắc: Trung tâm Dự báo khí tượng thủy văn Quốc gia xây dựng thuật toán tích hợp tự động kết quả xác định tâm của ADT với dữ liệu radar thời tiết ven biển (như trạm radar Phủ Liễn, Vinh, Quy Nhơn) và hệ thống trạm đảo tiền tiêu (Bạch Long Vĩ, Hoàng Sa, Trường Sa) trong lộ trình 6 - 12 tháng, nhằm bù đắp sai số vị trí 44 km của dạng mây SHEAR trong giai đoạn áp thấp nhiệt đới gần bờ.

  4. Đào tạo chuẩn hóa kỹ năng khai thác công nghệ: Trường Đại học Khoa học Tự nhiên phối hợp tổ chức các khóa tập huấn chuyên sâu định kỳ từ 3 đến 6 tháng về kỹ thuật Dvorak số hóa và quản trị hệ thống ADT cho tối thiểu 50 dự báo viên khí tượng trên toàn mạng lưới dự báo quốc gia.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả thực nghiệm của luận văn là tài liệu tham khảo giá trị cao cho 4 nhóm đối tượng chính:

  1. Dự báo viên khí tượng tại các cơ quan dự báo quốc gia và đài khí tượng thủy văn khu vực: Sử dụng tài liệu như một cẩm nang đối sánh sai số thực tế giữa phương pháp Dvorak thủ công và tự động, giúp nâng cao độ chính xác khi xác định cường độ các cơn bão khẩn cấp gần bờ.

  2. Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Khí tượng, Khí hậu học và Hải dương học: Khai thác khung lý thuyết toàn diện về kỹ thuật Dvorak hiện đại, cấu trúc dữ liệu ảnh vệ tinh số và phương pháp luận đánh giá sai số bão chuẩn quốc tế phục vụ các đề tài nghiên cứu chuyên sâu.

  3. Kỹ sư phát triển phần mềm và hệ thống thông tin khí tượng: Nắm vững sơ đồ thuật toán, quy trình xử lý ảnh hồng ngoại số hóa và cơ chế tương tác dữ liệu NetCDF trong môi trường Linux để xây dựng các công cụ giám sát thời tiết thông minh.

  4. Cán bộ quản lý thuộc Ban Chỉ đạo Quốc gia về Phòng chống Thiên tai: Tiếp cận các cơ sở khoa học định lượng về độ tin cậy của bản tin bão sớm, từ đó xây dựng kịch bản ứng phó, di dời dân cư và bảo vệ công trình ven biển hiệu quả, tiết kiệm nguồn lực.

Câu hỏi thường gặp

1. Phương pháp Dvorak cải tiến (ADT) có ưu thế gì vượt trội so với phương pháp Dvorak cổ điển?
Phương pháp ADT vận hành hoàn toàn tự động dựa trên thuật toán xử lý ảnh số hóa và phương trình hồi quy khách quan, loại bỏ định kiến chủ quan của dự báo viên. ADT có khả năng phân tích ảnh vệ tinh liên tục với tần suất cao (30 phút/lần), tính toán tự động bán kính gió mạnh (RMW) và hiệu chỉnh sai số áp suất theo vĩ độ.

2. Tại sao sai số vị trí tâm bão của ADT lại nhỏ nhất ở các cơn bão mạnh có mắt?
Khi bão đạt cường độ mạnh hình thành mắt bão rõ rệt (dạng EYE), độ tương phản nhiệt độ giữa tâm mắt ấm và vòng mây đối lưu lạnh xung quanh rất sắc nét trên kênh hồng ngoại. Thuật toán tích hợp vòng (RF) và phép biến đổi Fourier của ADT dễ dàng bắt trúng tọa độ hình học của tâm bão, giúp sai số giảm xuống mức trung bình 28 km.

3. Vì sao ADT có xu hướng ước tính tốc độ gió cực đại Vmax cao hơn số liệu Best Track Việt Nam?
Nguyên nhân chủ yếu do các phương trình hồi quy và bảng quy đổi T-number sang Vmax của ADT được xây dựng dựa trên mẫu khí hậu bão ngoài đại dương mở. Khi bão tiến vào Biển Đông chịu ảnh hưởng ma sát địa hình và sự xâm nhập của không khí lạnh tầng thấp, cấu trúc gió bị suy giảm nhanh hơn so với sự biến đổi hình thái mây trên ảnh vệ tinh.

4. Thuật toán ADT xử lý dữ liệu như thế nào khi bão bắt đầu đổ bộ vào đất liền?
Khi tâm bão tiếp cận đường bờ và đi vào đất liền, ADT tự động đọc file dữ liệu địa hình (TOPO), ngừng các bước tính toán cường độ đối lưu (Vmax và Pmin) và chỉ xuất ra thông tin tọa độ vị trí tâm bão. Điều này nhằm tránh các sai số lớn do ma sát bề mặt địa hình làm phá vỡ cấu trúc hoàn lưu mây tiêu chuẩn.

5. Luận văn đã sử dụng những nguồn dữ liệu vệ tinh nào để kiểm nghiệm thuật toán?
Nghiên cứu sử dụng chuỗi dữ liệu ảnh định dạng NetCDF từ vệ tinh địa tĩnh MTSAT-2 của Nhật Bản giai đoạn 2010 - 2015, bao gồm 5 kênh phổ: kênh thị phổ VIS (độ phân giải 1 km), hai kênh hồng ngoại nhiệt IR1 - IR2, kênh hơi nước WV và cận hồng ngoại IR4 (độ phân giải 4 km) với tần suất quét 30 phút một lần.

Kết luận

  • Luận văn đã hoàn thành việc đánh giá toàn diện phương pháp Dvorak cải tiến (ADT V8) trên tập dữ liệu gồm 41 cơn bão với 632 chu kỳ phân tích trong 6 năm (2010 - 2015) trên Biển Đông.
  • Xác lập bức tranh sai số định lượng rõ ràng: sai số vị trí tâm bão đạt độ chính xác cao nhất ở bão có mắt (28 km) và bão cấp 12 trở lên (33 km); ước tính áp suất trung tâm Pmin có độ tin cậy vượt trội với sai số dưới 4 mb.
  • Chỉ rõ hạn chế ước tính thiên cao của tốc độ gió cực đại Vmax (chênh lệch trung bình +24 kts ở bão có mắt), chỉ ra nguyên nhân do sự khác biệt giữa điều kiện khí hậu Biển Đông và đại dương mở.
  • Đặt nền tảng khoa học vững chắc cho việc chuyển giao công nghệ phân tích bão tự động từ môi trường nghiên cứu sang ứng dụng nghiệp vụ thực tế tại Trung tâm Dự báo khí tượng thủy văn Quốc gia.
  • Lộ trình 12 tháng tiếp theo cần tập trung hiệu chỉnh hệ số phương trình hồi quy địa phương hóa và tích hợp chuỗi dữ liệu radar - trạm đảo nhằm tối ưu hóa độ chính xác của hệ thống cảnh báo sớm thiên tai tại Việt Nam.