Tổng quan nghiên cứu

Khu vực bờ biển Việt Nam trải dài hơn 3.200 km, nằm trọn trong vùng biển Tây Bắc Thái Bình Dương – một trong những ổ bão nhiệt đới hoạt động mạnh mẽ và phức tạp nhất trên thế giới. Thống kê giai đoạn 1961 - 2014 ghi nhận có tới 364 cơn bão và áp thấp nhiệt đới đổ bộ hoặc ảnh hưởng trực tiếp đến nước ta, trung bình mỗi năm bờ biển Việt Nam phải hứng chịu trên 7 cơn bão và áp thấp nhiệt đới. Việc xác định chính xác vị trí tâm và cường độ bão trong thời gian thực giữ vai trò quyết định đối với hiệu quả của công tác phòng chống thiên tai. Thực tế các đợt thiên tai cho thấy sự sai lệch trong phân tích ban đầu giữa các trung tâm khí tượng quốc tế thường rất lớn. Điển hình như cơn bão số 1 (Kujira) năm 2015 khi tiến gần đảo Bạch Long Vĩ, mức ước lượng cường độ giữa Hoa Kỳ (cấp 7) và Hàn Quốc (xấp xỉ cấp 10) chênh lệch tới gần 3 cấp bão, gây khó khăn lớn cho công tác chỉ đạo ứng phó.

Vấn đề nghiên cứu trọng tâm của luận văn là giải quyết sự phụ thuộc vào tính chủ quan của dự báo viên trong phương pháp Dvorak cổ điển đang áp dụng tại Việt Nam, đồng thời rút ngắn khoảng cách công nghệ kéo dài hơn 20 năm so với các trung tâm dự báo tiên tiến trên thế giới. Mục tiêu cụ thể là kiểm chứng, đánh giá toàn diện sai số vị trí và cường độ bão của phương pháp Dvorak cải tiến (Advanced Dvorak Technique - ADT) so với số liệu chuẩn Best track của Việt Nam. Nghiên cứu được triển khai trên toàn vùng Biển Đông với tập dữ liệu gồm 41 cơn bão hoạt động trong giai đoạn 6 năm từ 2010 đến 2015, tương ứng 632 chu kỳ quan trắc định kỳ 6 giờ/lần. Kết quả đánh giá cung cấp luận cứ khoa học tin cậy để đưa thuật toán ADT tự động vào quy trình nghiệp vụ thời gian thực, nâng cao độ chính xác của các bản tin bão phát nhanh hàng giờ.

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

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng lý thuyết xác định cường độ xoáy thuận nhiệt đới từ ảnh mây vệ tinh do Vernon Dvorak khởi xướng năm 1975, trải qua bốn cột mốc phát triển công nghệ gồm: phương pháp Dvorak cổ điển (DT), phương pháp Dvorak khách quan (ODT), phương pháp Dvorak khách quan có cải tiến (AODT) và phương pháp Dvorak cải tiến (ADT) do Trung tâm CIMSS thuộc Đại học Wisconsin (Hoa Kỳ) hoàn thiện. Ba khái niệm cốt lõi được sử dụng xuyên suốt gồm: Chỉ số T-number (đại lượng định lượng phản ánh mức độ phát triển cấu trúc mây), Chỉ số Cường độ hiện tại CI (Current Intensity - giá trị quy đổi sang tốc độ gió cực đại $V_{max}$ và áp suất mực biển cực tiểu $P_{min}$), và Kỹ thuật nhận dạng mẫu mây tự động bằng phân tích phổ biến đổi Fourier (FFT).

Mô hình phân loại của ADT chia cấu trúc xoáy thuận nhiệt đới thành 6 dạng mẫu mây 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 dày đặc trung tâm có nhiệt độ đồng đều (UNIFRM), khối mây trung tâm biến động lớn (IRRCDO) và dạng mắt bão (EYE). Bên cạnh đó, ADT tích hợp hai thuật toán mở rộng quan trọng là sơ đồ hiệu chỉnh áp suất mực biển theo vĩ độ (MSLP) và thuật toán ước tính bán kính gió mạnh cực đại (RMW) dựa trên trường nhiệt độ đỉnh mây hồng ngoại.

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

Nguồn dữ liệu nghiên cứu sử dụng chuỗi ảnh từ vệ tinh địa tĩnh MTSAT-2 của Cơ quan Khí tượng Nhật Bản (JMA) trong suốt giai đoạn 2010 - 2015. Tập dữ liệu ảnh vệ tinh bao gồm 5 kênh phổ: kênh thị phổ VIS (độ phân giải không gian 1 km), kênh hồng ngoại nhiệt 1 và 2 (IR1, IR2), kênh hơi nước (WV) và kênh cận hồng ngoại (IR4) cùng có độ phân giải 4 km, với tần suất quét 30 phút/sản phẩm trên phạm vi tọa độ từ -5°N đến 40°N và từ 95°E đến 140°E định dạng NetCDF. Cỡ mẫu nghiên cứu gồm toàn bộ 41 cơn bão trên Biển Đông trong 6 năm (2010 - 2015) với tổng số 632 chu kỳ dự báo chuẩn 6 giờ/lần.

Phương pháp chọn mẫu áp dụng nguyên tắc lấy mẫu toàn bộ các chu kỳ khi tâm bão hoạt động hoàn toàn trên biển và chủ động loại bỏ các chu kỳ khi bão đã đổ bộ vào đất liền, do thuật toán ADT chỉ tối ưu hóa các tham số nhiệt động lực học trên mặt nước biển. Lý do lựa chọn hệ thống chỉ số thống kê gồm Sai số tuyệt đối trung bình (MAE), Độ lệch (BIAS), Căn bậc hai sai số bình phương trung bình (RMSE) và Hệ số tương quan Pearson (Cor) là nhằm đánh giá đa chiều cả về độ lớn sai số, xu hướng thiên lệch và mức độ đồng biến giữa kết quả tính toán tự động của ADT với số liệu Best track Việt Nam và phương pháp Dvorak thủ công trên hai cơn bão điển hình là Megi (2010) và Rammasun (2014).

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

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

Thứ nhất, về độ chính xác xác định vị trí tâm bão, kết quả kiểm chứng 632 chu kỳ chỉ ra rằng sai số vị trí phụ thuộc mật thiết vào độ hoàn thiện của cấu trúc mây. Dạng mắt bão (EYE) đạt độ chính xác cao nhất với sai số vị trí trung bình chỉ 28 km. Ngược lại, dạng lệch tâm (SHEAR) có sai số vị trí lớn nhất lên tới 44 km. Độ chính xác định vị tâm bão của ADT giảm dần theo thứ tự hình thái: dạng mắt bão EYE (28 km), dạng tâm nhúng EMBC (31 km), dạng khối mây đồng đều UNIFRM (35 km), dạng khối mây biến động IRRCDO (38 km), dạng băng cuốn CRVBND (41 km) và dạng lệch tâm SHEAR (44 km).

Thứ hai, về tốc độ gió cực đại ($V_{max}$), dạng băng cuốn (CRVBND) cho kết quả sát thực tế nhất với mức chênh lệch trung bình giữa ADT và Best track chỉ khoảng 4 kts (tương ứng tỷ lệ sai số 9,6%, sai số tuyệt đối trung bình MAE = 6,4 kts). Ngược lại, dạng mắt bão (EYE) có mức sai số lớn nhất khi ADT ước lượng tốc độ gió cao hơn Best track trung bình tới 24 kts (tỷ lệ sai số 27,6%, MAE = 24,1 kts). Các dạng mây UNIFRM, IRRCDO và EMBC có xu hướng ước lượng tốc độ gió cao hơn Best track từ 5 đến 10 kts.

Thứ ba, về trị số khí áp thấp nhất ($P_{min}$), phương pháp ADT đạt độ tin cậy rất cao trên toàn bộ 6 dạng mẫu mây với tỷ lệ sai số trung bình cực thấp, chỉ dao động trong khoảng từ 0,26% đến 0,40%. Chênh lệch giá trị trung bình giữa ADT và Best track duy trì ổn định dưới ngưỡng 4,5 mb (chỉ số MAE dao động từ 4,2 mb đến 5,8 mb).

Thảo luận kết quả

Khi biểu diễn dữ liệu thực nghiệm qua các biểu đồ phân tán (scatter plots) và phân tích tương quan hồi quy, đường hồi quy tuyến tính của dạng băng cuốn (CRVBND) gần như trùng khít và chạy song song với đường tham chiếu chuẩn 1:1, chứng minh thuật toán phân tích vòng xoắn logarit 10 độ của ADT hoạt động cực kỳ chính xác ở các cơn bão có cường độ trung bình (cấp 8 đến cấp 9). Hiện tượng ADT ước lượng gió cực đại cao hơn Best track khoảng 24 kts ở các cơn bão mạnh có dạng mắt (EYE) bắt nguồn từ việc bộ tham số hồi quy gốc của CIMSS được xây dựng dựa trên dữ liệu thám sát máy bay tại các vùng biển khơi rộng lớn ở Đại Tây Dương và Tây Bắc Thái Bình Dương. Khi bão tiến vào Biển Đông – một vùng biển hẹp có rào chắn địa hình phức tạp, hoàn lưu bão chịu ma sát đáy lớn làm tốc độ gió thực tế suy giảm nhanh hơn so với cấu trúc nhiệt độ đỉnh mây ghi nhận trên ảnh hồng ngoại.

Đối với hai cơn bão điển hình là Megi (năm 2010) và Rammasun (năm 2014), biểu đồ biến thiên theo trục thời gian cho thấy đường diễn biến áp suất và tốc độ gió của ADT nắm bắt rất nhạy các giai đoạn bão tăng cường độ đột ngột (Rapid Intensification). So sánh với phương pháp Dvorak thủ công, ADT loại bỏ hoàn toàn hiện tượng trễ thông tin từ 6 đến 12 giờ do tâm lý e ngại thay đổi cấp bão quá nhanh của dự báo viên, khẳng định ưu thế vượt trội của thuật toán tự động khách quan.

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

Thứ nhất, tích hợp trực tiếp module tính toán ADT tự động vào hệ thống dự báo số trị thời gian thực tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia. Cần thiết lập đường truyền số liệu tự động từ các vệ tinh địa tĩnh thế hệ mới với chu kỳ quét liên tục nhằm xuất bản kết quả vị trí và cường độ bão với tần suất 60 phút/lần, hoàn thành triển khai thử nghiệm diện rộng trong giai đoạn 2024 - 2025.

Thứ hai, hiệu chỉnh cục bộ bộ hệ số phương trình hồi quy chuyển đổi T-number sang tốc độ gió cực đại ($V_{max}$) cho riêng khu vực Biển Đông. Viện Khoa học Khí tượng Thủy văn và Hải văn cần phối hợp với các trường đại học áp dụng hệ số suy giảm từ 10% đến 15% đối với các cơn bão đạt cấp 12 trở lên có dạng mắt (EYE), đặt mục tiêu giảm sai số MAE của tốc độ gió xuống dưới 12 kts trước mùa bão năm 2025.

Thứ ba, phát triển thuật toán lọc nhiễu và định vị tâm bổ trợ cho các trường hợp mây bão dạng lệch tâm (SHEAR). Cần kết hợp dữ liệu ảnh hồng ngoại với dữ liệu mạng lưới radar thời tiết ven biển và các trạm quan trắc hải đảo (như đảo Bạch Long Vĩ, Song Tử Tây) để giảm sai số vị trí tâm bão ở giai đoạn áp thấp nhiệt đới từ 44 km xuống dưới 30 km.

Thứ tư, chuẩn hóa chương trình tập huấn và chuyển giao công nghệ phân tích ảnh vệ tinh tự động cho 100% dự báo viên bão tại các Đài Khí tượng Thủy văn khu vực. Dự báo viên cần nắm vững quy tắc phối hợp giữa kết quả tự động ADT và kinh nghiệm phân tích địa phương nhằm nâng cao chất lượng bản tin bão khẩn cấp phục vụ cộng đồng từ năm 2026.

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

Nhóm 1 - Dự báo viên bão tại các cơ quan khí tượng thủy văn: Sử dụng luận văn như tài liệu hướng dẫn kỹ thuật chi tiết để nắm rõ bản chất thuật toán ADT, nhận diện ưu điểm và hạn chế của 6 dạng mẫu mây, từ đó đưa ra quyết định hiệu chỉnh cường độ gió và khí áp bão một cách chuẩn xác trong các ca trực nghiệp vụ.

Nhóm 2 - Học viên cao học và nhà nghiên cứu chuyên ngành Khí tượng, Hải dương học: Luận văn cung cấp phương pháp luận chuẩn mực về đánh giá sai số mô hình với bộ dữ liệu thực nghiệm 41 cơn bão và 632 chu kỳ quan trắc, là tài liệu tham khảo giá trị cho các công trình nghiên cứu viễn thám và động lực học bão nhiệt đới.

Nhóm 3 - Ban Chỉ đạo Quốc gia về Phòng, chống thiên tai và các cơ quan ứng phó khẩn cấp: Tiếp cận cơ sở khoa học về mức độ sai lệch của các nguồn số liệu bão quốc tế và trong nước, phục vụ việc phân tích rủi ro, lập kịch bản sơ tán dân cư ven biển và điều hành phương án bảo vệ đê điều, hồ đập kịp thời.

Nhóm 4 - Kỹ sư phát triển phần mềm và công nghệ môi trường: Khai thác cấu trúc xử lý dữ liệu viễn thám NetCDF, thuật toán phân tích Fourier (FFT) và logic tự động hóa trên hệ điều hành Linux để xây dựng các nền tảng giám sát thời tiết thông minh tại Việt Nam.

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

Phương pháp ADT có ưu điểm vượt trội nào so với kỹ thuật Dvorak truyền thống? Phương pháp ADT tự động hóa 100% quy trình tính toán từ định vị tâm bão đến nhận dạng mẫu mây bằng thuật toán biến đổi Fourier. Kỹ thuật này loại bỏ hoàn toàn tính chủ quan và sự chênh lệch trình độ giữa các dự báo viên, đồng thời cho phép trích xuất kết quả cường độ bão tự động sau mỗi 30 phút thay vì phải đo đạc thủ công mất nhiều thời gian.

Tại sao dạng mây lệch tâm (SHEAR) lại có sai số vị trí tâm bão lớn nhất (khoảng 44 km)? Dạng SHEAR thường xuất hiện khi xoáy thuận mới hình thành hoặc chịu gió đứt thẳng đứng lớn làm khối mây đối lưu bị đẩy lệch xa khỏi tâm xoáy tầng thấp. Do tâm hoàn lưu bị che khuất và thiếu cấu trúc xoắn đồng tâm rõ nét, thuật toán định vị quang học gặp khó khăn trong việc nhận diện điểm tâm, dẫn đến sai số vị trí trung bình lên tới 44 km.

Nguyên nhân nào khiến ADT ước lượng tốc độ gió cực đại ở dạng mắt (EYE) cao hơn Best track khoảng 24 kts? Công thức hồi quy của ADT được xây dựng dựa trên mẫu bão ngoài đại dương mở. Khi bão đi vào Biển Đông, tác động ma sát đáy biển nông và hoàn lưu lục địa làm giảm tốc độ gió bề mặt nhanh hơn tốc độ suy giảm nhiệt độ đỉnh mây, khiến thuật toán ước lượng tốc độ gió cao hơn số liệu thực tế khoảng 27,6%.

Nguồn dữ liệu vệ tinh nào đã được sử dụng để kiểm nghiệm thuật toán trong nghiên cứu? Luận văn sử dụng nguồn dữ liệu ảnh vệ tinh địa tĩnh MTSAT-2 của Cơ quan Khí tượng Nhật Bản với 5 kênh phổ (VIS, IR1, IR2, WV, IR4), độ phân giải từ 1 km đến 4 km. Toàn bộ chuỗi dữ liệu gồm 41 cơn bão giai đoạn 2010 - 2015 với 632 chu kỳ tính toán đã được chạy nghiệm thu đầy đủ trên môi trường Linux.

Ý nghĩa thực tiễn lớn nhất của việc áp dụng thuật toán ADT tại Việt Nam là gì? Ý nghĩa lớn nhất là hiện đại hóa công nghệ dự báo bão, cung cấp công cụ khách quan để xác định vị trí và cường độ bão liên tục phục vụ bản tin bão nhanh phát 1 giờ/lần. Điều này giúp giảm độ trễ thông tin, hỗ trợ đắc lực cho công tác chỉ đạo ứng phó thiên tai giảm thiểu thiệt hại kinh tế.

Kết luận

  • Luận văn đã thực hiện đánh giá toàn diện và hệ thống độ chính xác của phương pháp Dvorak cải tiến (ADT) trên chuỗi dữ liệu thực nghiệm gồm 41 cơn bão với 632 chu kỳ quan trắc trong 6 năm (2010 - 2015) trên khu vực Biển Đông.
  • Kết quả khẳng định ADT có độ tin cậy vượt trội trong việc xác định vị trí bão dạng mắt (EYE) với sai số chỉ 28 km và ước tính khí áp cực tiểu ($P_{min}$) với sai số trung bình luôn duy trì dưới 4,5 mb trên tất cả 6 dạng mẫu mây.
  • Luận văn chỉ ra rõ các giới hạn cần lưu ý: ADT ước lượng tốc độ gió cực đại ($V_{max}$) cao hơn thực tế khoảng 24 kts ở các cơn bão mạnh có mắt và sai số vị trí còn lớn ở dạng lệch tâm (khoảng 44 km).
  • Công trình đã đóng góp cơ sở khoa học vững chắc giúp rút ngắn khoảng cách công nghệ hơn 20 năm của ngành khí tượng Việt Nam, mở đường cho việc tự động hóa hoàn toàn quy trình phân tích bão thời gian thực.
  • Kế hoạch tiếp theo cần tập trung hiệu chỉnh bộ tham số hồi quy cho riêng Biển Đông trên chuỗi dữ liệu vệ tinh phân giải cao thế hệ mới. Độc giả và các đơn vị chuyên môn có thể khai thác trực tiếp các bộ số liệu thống kê trong luận văn để phục vụ công tác nghiên cứu và triển khai nghiệp vụ dự báo bão thực tế.