Tổng quan nghiên cứu

Khu vực Biển Đông trung bình mỗi năm đón nhận từ 10 đến 12 cơn bão và áp thấp nhiệt đới, gây ra những tổn thất nặng nề về người và tài sản với thiệt hại ước tính hàng nghìn tỷ đồng cho các tỉnh ven biển Việt Nam. Trong bối cảnh biến đổi khí hậu toàn cầu, các cơn bão có xu hướng gia tăng về tần suất cực đoan, cường độ tàn phá và quỹ đạo di chuyển ngày càng dị thường, phức tạp. Trước đây, các hệ thống cảnh báo khí tượng trong nước chủ yếu tập trung vào hạn dự báo ngắn từ 24 đến 72 giờ dựa trên phương pháp synop truyền thống hoặc mô hình đơn lẻ, vốn bộc lộ nhiều hạn chế khi thời gian dự báo kéo dài do sai số tích lũy từ trường khí quyển ban đầu.

Vấn đề cốt lõi đặt ra là nâng cao độ tin cậy và kéo dài hạn dự báo bão lên 120 giờ (tương đương 5 ngày) nhằm tạo khoảng thời gian vàng cho công tác chỉ đạo phòng chống thiên tai. Mục tiêu cụ thể của công trình là nghiên cứu thử nghiệm và đánh giá hệ thống 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 mô hình khí tượng quy mô vừa WRF (Weather Research and Forecasting), tích hợp các sản phẩm dự báo tổ hợp toàn cầu từ Trung tâm Dự báo Môi trường Quốc gia Hoa Kỳ (NCEP/NOAA).

Phạm vi nghiên cứu bao quát toàn bộ vùng Biển Đông và Tây Bắc Thái Bình Dương trong giai đoạn 2009 đến 2011 với tập dữ liệu thực nghiệm gồm 11 cơn bão nhiệt đới tiêu biểu. Ý nghĩa thực tiễn của đề tài được lượng hóa rõ rệt: việc kéo dài thời gian dự báo thêm 48 giờ giúp giảm thiểu hơn 25% rủi ro cho các đội tàu thuyền xa bờ, hạ thấp sai số dự báo quỹ đạo hạn 72 giờ xuống 262 km và kiểm soát sai số áp suất trung tâm bão trong khoảng 8 đến 16 mb, đóng góp một giải pháp công nghệ tiên tiến cho hệ thống cảnh báo sớm quốc gia.

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 lý thuyết hỗn độn khí quyển do Edward Lorenz đề xuất năm 1963, giải thích rằng trạng thái khí quyển là một hệ phi tuyến tính nhạy cảm với điều kiện ban đầu, dẫn đến hiệu ứng cánh bướm làm phân kỳ các kết quả dự báo sau một khoảng thời gian tích phân từ 10 đến 14 ngày. Bên cạnh đó, luận văn ứng dụng lý thuyết mô hình hóa động lực học khí quyển quy mô vừa và lý thuyết dự báo tổ hợp (Ensemble Forecasting) nhằm loại bỏ sai số ngẫu nhiên bằng cách tích phân nhiều kịch bản điều kiện ban đầu khác nhau.

Mô hình nghiên cứu nòng cốt là WRF phiên bản ARW (Advanced Research WRF) với cấu hình miền tính gồm 151x151 điểm lưới theo phương ngang, độ phân giải không gian 27 km và 31 tầng áp suất theo phương thẳng đứng. Tâm miền tính được định vị linh hoạt theo tọa độ tâm bão thực tế. Ba khái niệm chuyên ngành then chốt được chuẩn hóa bao gồm:

  • Phương pháp siêu tổ hợp (Superensemble): kỹ thuật gán trọng số hồi quy đa biến dựa trên dữ liệu quá khứ để tối ưu hóa dự báo vị trí và cường độ bão.
  • Sai số quỹ đạo khoảng cách (PE), sai số dọc quỹ đạo (ATE) và sai số ngang quỹ đạo (CTE) theo hướng di chuyển của bão.
  • Các chỉ số cường độ bão chuẩn: áp suất mực biển cực tiểu (MSLP, tính bằng mb) và vận tốc gió cực đại bề mặt (Vmax, tính bằng m/s).

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

Nguồn dữ liệu đầu vào sử dụng sản phẩm dự báo tổ hợp toàn cầu NCEP/NOAA gồm 21 thành phần, độ phân giải 1x1 độ, chứa 29 biến khí tượng tại 27 tầng thẳng đứng và cập nhật định kỳ mỗi 6 giờ. Dữ liệu kiểm chứng quỹ đạo và cường độ thực tế được trích xuất từ cơ sở dữ liệu Best Track của Trung tâm Khí tượng Nhật Bản (JMA) thông qua hệ thống Digital Typhoon.

Cỡ mẫu nghiên cứu bao gồm 11 cơn bão hoạt động từ 5 ngày trở lên trong giai đoạn 2009-2011 với 30 trường hợp dự báo thử nghiệm độc lập. Phương pháp chọn mẫu là chọn mẫu có chủ đích, ưu tiên các cơn bão có quỹ đạo phức tạp, cường độ biến động mạnh trên Biển Đông để kiểm tra sức chịu tải của mô hình. Nhằm xây dựng hệ phương trình hồi quy cho phương pháp siêu tổ hợp, bộ số liệu phụ thuộc được mở rộng lên 90 mẫu bằng cách tích hợp kết quả dự báo tại các thời điểm lệch nhau 3 giờ, dựa trên giả thiết hệ số hồi quy ổn định trong khoảng 6 giờ.

Lý do lựa chọn phương pháp phân tích kết hợp giữa mô hình động lực WRF và kỹ thuật siêu tổ hợp là nhằm khắc phục triệt để hạn chế tài nguyên máy tính. Thay vì chạy toàn bộ 21 thành phần NCEP gây quá tải hệ thống, mô hình chỉ chọn 5 thành phần đại diện kết hợp với 3 sơ đồ tham số hóa đối lưu mây tích (Kain-Fritsch, Betts-Miller-Janjic, Grell-Devenyi) để tạo thành tổ hợp 15 kịch bản dự báo tối ưu. Toàn bộ chu trình tính toán và kiểm nghiệm độc lập được tiến hành xuyên suốt trong 24 tháng.

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

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

Thứ nhất, việc khảo sát số lượng thành phần trên cơn bão Conson chỉ ra rằng cấu hình 5 thành phần kết hợp 3 sơ đồ đối lưu cho sai số quỹ đạo hạn 120 giờ là 681 km, chỉ chênh lệch khoảng 10 km (tương đương 1,4%) so với phương án sử dụng đầy đủ 21 thành phần (680 km), trong khi tiết kiệm được hơn 65% thời gian xử lý dữ liệu.

Thứ hai, phương pháp tổ hợp trung bình đơn giản với 15 thành phần trên toàn bộ 30 trường hợp thử nghiệm đã ghi nhận sai số quỹ đạo trung bình đạt 100,5 km ở hạn 24 giờ, giảm xuống còn 262 km ở hạn 72 giờ và đạt 508 km ở hạn 120 giờ. Kết quả này vượt trội so với các nghiên cứu trước đây tại Việt Nam, giảm được 60 km sai số ở hạn 72 giờ (tương đương mức cải thiện 18,6% so với sai số 322 km của các mô hình truyền thống).

Thứ ba, phương pháp siêu tổ hợp khẳng định ưu thế vượt trội khi hạ sai số quỹ đạo hạn 24 giờ xuống 76,18 km và hạn 120 giờ xuống khoảng 300 đến 320 km khi sử dụng 12 thành phần hồi quy. Mức sai số này giảm hơn 37% so với phương pháp trung bình đơn giản ở cùng hạn 5 ngày. Ngay cả với các trạm dự báo có năng lực máy tính vừa phải, việc chỉ sử dụng 4 thành phần siêu tổ hợp vẫn duy trì được sai số hạn 120 giờ ở mức khoảng 400 km.

Thứ tư, đối với dự báo cường độ, sai số áp suất mực biển cực tiểu trung bình dao động từ 8,23 mb đến 16 mb, trong khi sai số vận tốc gió cực đại dao động từ -8,5 m/s đến 0 m/s. Đáng chú ý, sai số cường độ có dạng phân bố hình parabol, đạt độ chính xác cao nhất ở hai mốc hạn dự báo 24 giờ và 120 giờ.

Thảo luận kết quả

Nguyên nhân chính dẫn đến việc sai số vận tốc gió cực đại thường mang giá trị âm (mô hình dự báo gió yếu hơn thực tế) và sai số áp suất mang giá trị dương (áp suất dự báo cao hơn thực tế) là do bước lưới ngang 27 km chưa đủ mịn để mô phỏng trọn vẹn cấu trúc đối lưu sâu trong mắt bão, nơi có độ dốc khí áp cực lớn. Ngoài ra, sự thiếu hụt trạm thám không và phao quan trắc trên Biển Đông khiến trường ẩm và nhiệt ban đầu chưa phản ánh hoàn hảo tương tác biển - khí quyển.

Khi so sánh với các nghiên cứu quốc tế, kết quả sai số quỹ đạo 24 giờ (76,18 km) và 48 giờ (khoảng 213 km) của phương pháp siêu tổ hợp trong luận văn hoàn toàn cạnh tranh với các trung tâm lớn như JMA (147 km ở hạn 24 giờ) và ECMWF (214 km ở hạn 24 giờ).

Trong thực tế phân tích nghiệp vụ, các tập số liệu này được biểu diễn tối ưu thông qua hệ thống bản đồ mạng nhện (Spaghetti plots) trực quan hóa độ phân tán quỹ đạo bão, kết hợp bảng ma trận sai số ATE và CTE giúp chuyên viên nhận diện xu thế lệch phải hoặc lệch bắc của tâm bão dự báo so với quỹ đạo Best Track. Biểu đồ cột kép thể hiện tương quan giữa sai số áp suất cực tiểu và tốc độ gió cực đại cũng chứng minh rõ ràng khả năng bám sát chu kỳ sống của các cơn bão mạnh như Utor hay Conson.

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

Thứ nhất, chuẩn hóa và tích hợp quy trình dự báo siêu tổ hợp WRF 12 thành phần vào hệ thống nghiệp vụ của Trung tâm Dự báo Khí tượng Thủy văn Quốc gia trước quý 4 năm 2027. Mục tiêu là nâng cao độ chính xác của các bản tin cảnh báo bão sớm 5 ngày, giảm sai số vị trí tâm bão xuống dưới 280 km ở hạn 72 giờ và dưới 450 km ở hạn 120 giờ.

Thứ hai, triển khai ứng dụng sơ đồ đồng hóa số liệu ba chiều (3DVAR) tích hợp dữ liệu vệ tinh địa tĩnh độ phân giải cao và mạng lưới radar thời tiết ven biển vào mô hình WRF. Giải pháp này do Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu chủ trì thực hiện trong giai đoạn 2027-2028, hướng tới mục tiêu giảm ít nhất 15% sai số ban đầu của trường khí áp và độ ẩm trung tâm xoáy thuận.

Thứ ba, đầu tư nâng cấp hệ thống siêu máy tính tính toán song song (HPC) với cấu hình tối thiểu 128 lõi tính toán và 512 GB RAM tại các đài khí tượng thủy văn khu vực trọng điểm. Dự án cần hoàn thành trong vòng 18 tháng nhằm cho phép thu nhỏ độ phân giải miền tính từ 27 km xuống 9 km hoặc 3 km, giúp mô phỏng chân thực các biến động cường độ gió cực đại ven bờ.

Thứ tư, xây dựng bộ cơ sở dữ liệu số hóa đồng bộ về bão lịch sử và tổ chức 04 khóa đào tạo chuyên sâu về kỹ thuật siêu tổ hợp cho 60 dự báo viên trên toàn quốc từ nay đến năm 2028. Cơ quan quản lý thiên tai cấp tỉnh cần căn cứ vào biên độ sai số dự báo 120 giờ của mô hình để thiết lập kịch bản sơ tán dân cư trước 72 giờ thay vì 24 giờ như quy trình cũ.

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

Nhóm thứ nhất bao gồm các chuyên viên dự báo thời tiết tại các cơ quan khí tượng thủy văn quốc gia và khu vực. Luận văn cung cấp tài liệu kỹ thuật hoàn chỉnh để thiết lập cấu hình mô hình WRF, lựa chọn sơ đồ tham số hóa đối lưu phù hợp và khai thác sản phẩm NCEP toàn cầu phục vụ việc ra bản tin cảnh báo bão nghiệp vụ hạn 5 ngày.

Nhóm thứ hai là các nhà khoa học, nghiên cứu sinh và học viên cao học thuộc chuyên ngành Khí tượng học, Hải dương học và Khoa học Môi trường. Công trình là tài liệu tham khảo giá trị về phương pháp luận siêu tổ hợp đa mô hình, cách thức xử lý sai số thống kê ATE/CTE và mô phỏng số trị tương tác khí quyển quy mô vừa.

Nhóm thứ ba là Ban Chỉ đạo Quốc gia về Phòng chống thiên tai cùng các cơ quan quản lý cứu nạn cứu hộ. Dữ liệu đánh giá sai số quỹ đạo hạn 120 giờ cung cấp căn cứ khoa học xác thực để hoạch định phương án ứng phó, điều động lực lượng và chủ động khoanh vùng khu vực nguy hiểm trên biển từ sớm.

Nhóm thứ tư gồm các doanh nghiệp vận tải biển, khai thác dầu khí và nuôi trồng thủy sản xa bờ. Doanh nghiệp có thể ứng dụng các chỉ dẫn về sai số gió cực đại và áp suất cực tiểu để xây dựng quy trình quản trị rủi ro hàng hải, điều chỉnh hải trình tàu thuyền trước 5 ngày nhằm đảm bảo an toàn tuyệt đối cho người và phương tiện.

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

Phương pháp siêu tổ hợp khác biệt như thế nào so với phương pháp trung bình đơn giản? Phương pháp trung bình đơn giản chia đều trọng số 1/N cho mọi thành phần tham gia, dẫn đến rủi ro bị lệch nếu có thành phần dự báo kém. Ngược lại, phương pháp siêu tổ hợp sử dụng hồi quy đa biến có lọc trên dữ liệu quá khứ để gán trọng số cao cho các thành phần dự báo chính xác, giúp hạ sai số hạn 120 giờ từ 508 km xuống còn khoảng 300 đến 320 km.

Tại sao nghiên cứu chỉ lựa chọn 5 thành phần NCEP thay vì toàn bộ 21 thành phần toàn cầu? Thực nghiệm trên bão Conson cho thấy sai số quỹ đạo hạn 5 ngày của phương án 5 thành phần (681 km) và phương án 21 thành phần (680 km) chỉ chênh nhau vỏn vẹn 10 km. Việc chọn 5 thành phần kết hợp 3 sơ đồ đối lưu giúp tiết kiệm hơn 65% thời gian tính toán mà vẫn đảm bảo độ tán và chất lượng dự báo tối ưu.

Nguyên nhân nào khiến mô hình WRF thường dự báo cường độ bão yếu hơn so với thực tế? Hiện tượng này xuất phát từ độ phân giải không gian 27 km chưa đủ mịn để nắm bắt vùng lõi xoáy mạnh nhất quanh mắt bão. Kết quả thống kê cho thấy sai số áp suất cực tiểu Pmin thường cao hơn thực tế từ 8 đến 16 mb, trong khi vận tốc gió cực đại Vmax dự báo thấp hơn thực tế từ 0 đến 8,5 m/s.

Cấu hình miền tính 151x151 điểm lưới có đáp ứng tốt cho việc theo dõi bão trên Biển Đông không? Miền tính 151x151 điểm lưới với bước lưới 27 km bao phủ diện tích đủ rộng cho toàn bộ Biển Đông và vùng biển lân cận phía đông Philippines. Cấu hình 31 tầng áp suất đảm bảo mô phỏng chi tiết dòng dẫn đường tầng đối lưu giữa, giúp mô hình tích phân ổn định trong suốt 120 giờ dự báo.

Các trung tâm dự báo có hệ thống tính toán vừa phải có thể ứng dụng kết quả này được không? Hoàn toàn khả thi. Nghiên cứu đã chứng minh rằng phương án siêu tổ hợp chỉ cần 4 thành phần đã đạt sai số quỹ đạo hạn 120 giờ ở mức 400 km. Mức sai số này hoàn toàn đáp ứng tiêu chuẩn nghiệp vụ và tương thích tốt với năng lực của các máy chủ tính toán tầm trung.

Kết luận

  • Luận văn đã tiên phong thử nghiệm thành công mô hình WRF tích hợp sản phẩm tổ hợp toàn cầu NCEP/NOAA để dự báo quỹ đạo và cường độ bão hạn 5 ngày trên Biển Đông.
  • Xác lập cấu hình tối ưu gồm 5 thành phần tổ hợp kết hợp 3 sơ đồ đối lưu (Kain-Fritsch, Betts-Miller-Janjic, Grell-Devenyi), giúp tiết kiệm hơn 65% chi phí tính toán so với việc dùng toàn bộ 21 thành phần.
  • Phương pháp siêu tổ hợp 12 thành phần đạt độ chính xác vượt bậc với sai số quỹ đạo hạn 24 giờ là 76,18 km và hạn 120 giờ khoảng 300 đến 320 km, cải thiện hơn 37% so với phương pháp trung bình đơn giản.
  • Nhận diện và lượng hóa rõ ràng đặc tính sai số cường độ bão theo dạng hình parabol, với sai số áp suất dao động từ 8 đến 16 mb và sai số gió cực đại từ -8,5 đến 0 m/s.
  • Cung cấp cơ sở khoa học và bộ phương trình hồi quy vững chắc để chuyển giao công nghệ vào thực tiễn dự báo nghiệp vụ khí tượng thủy văn tại Việt Nam.

Đóng góp quan trọng nhất của công trình là mở rộng thành công hạn dự báo bão chi tiết lên 120 giờ với sai số được kiểm soát chặt chẽ. Kế hoạch tiếp theo trong 12 đến 24 tháng tới là mở rộng kiểm nghiệm trên các mùa bão mới và tích hợp công nghệ đồng hóa xoáy giả 3DVAR. Các cơ quan quản lý và đơn vị nghiên cứu hãy chủ động ứng dụng khung phương pháp siêu tổ hợp này để nâng cao hiệu quả phòng ngừa thiên tai và bảo vệ an toàn lãnh hải quốc gia.