Tổng quan nghiên cứu

Việt Nam sở hữu đường bờ biển dài hơn 3.260 km và hàng năm phải hứng chịu trung bình từ 5 đến 7 cơn bão cùng áp thấp nhiệt đới trực tiếp đổ bộ, gây thiệt hại kinh tế hàng nghìn tỷ đồng và đe dọa sinh mạng của hàng triệu cư dân ven biển. Trong bối cảnh biến đổi khí hậu toàn cầu khiến quỹ đạo và cường độ thiên tai ngày càng trở nên phức tạp, việc nâng cao năng lực dự báo sớm quỹ đạo và thời điểm bão tiếp bờ trở thành nhiệm vụ sống còn. Phần lớn các mô hình dự báo nghiệp vụ trước đây chỉ cung cấp kết quả tương đối tin cậy trong phạm vi hạn ngắn dưới 72 giờ, trong khi công tác sơ tán dân cư, gia cố đê điều và điều phối tàu thuyền trên biển đòi hỏi khoảng thời gian chuẩn bị dài hơn từ 3 đến 5 ngày (tương đương 72 đến 120 giờ).

Luận văn thạc sĩ khoa học chuyên ngành Khí tượng và Khí hậu học của tác giả Lã Thị Tuyết, dưới sự hướng dẫn khoa học của Giáo sư Tiến sĩ Trần Tân Tiến tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, đã tập trung giải quyết bài toán thử nghiệm dự báo vị trí và thời gian bão đổ bộ vào bờ biển Việt Nam trước 3 đến 5 ngày bằng mô hình số trị WRF.

Phạm vi nghiên cứu bao trùm toàn bộ dải bờ biển Việt Nam từ Bắc Bộ đến Nam Bộ với tập dữ liệu thực nghiệm gồm 7 cơn bão mạnh có quỹ đạo phức tạp giai đoạn 2005 - 2010. Ý nghĩa thực tiễn của công trình thể hiện qua việc đánh giá định lượng các chỉ tiêu sai số vị trí (MPE tính bằng km) và sai số thời gian tiếp bờ (tính bằng giờ), tạo cơ sở khoa học vững chắc giúp mở rộng thời gian cảnh báo thiên tai sớm thêm 48 đến 72 giờ so với các quy trình dự báo truyền thống.

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 động lực học khí quyển phi thủy tĩnh (Non-hydrostatic Atmospheric Dynamics) - cơ sở cốt lõi của lõi động lực học ARW trong mô hình nghiên cứu và dự báo thời tiết WRF. Bên cạnh đó, tác giả ứng dụng lý thuyết dòng dẫn đường lớp sâu (Deep Layer Mean Steering Flow Theory) để mô tả cơ chế di chuyển của hoàn lưu xoáy thuận nhiệt đới trong mối tương tác với các hệ thống hoàn lưu quy mô lớn như cao áp cận nhiệt đới Tây Thái Bình Dương và gió mùa Đông Á.

Hệ thống khái niệm chuyên ngành chủ đạo được sử dụng xuyên suốt công trình bao gồm:

  1. Mô hình số trị quy mô vừa WRF (Weather Research and Forecasting): Hệ thống mô hình dự báo thời tiết thế hệ mới tích hợp đầy đủ các quá trình động lực và tham số hóa vật lý đa dạng.
  2. Kỹ thuật ban đầu hóa xoáy nhân tạo (Bogus Vortex Initialization): Phương pháp tái tạo cấu trúc xoáy đối xứng trục và cân bằng động lực học nhằm bù đắp sự thiếu hụt dữ liệu quan trắc ngoài khơi xa.
  3. Sai số khoảng cách tâm bão trực tiếp (Direct Positional Error - DPE/PE): Khoảng cách trắc địa trên mặt cầu giữa tâm bão dự báo và tâm bão quan trắc thực tế.
  4. Lớp biên hành tinh (Planetary Boundary Layer - PBL) và thông lượng ẩm bề mặt: Yếu tố quyết định sự suy giảm năng lượng và biến đổi cấu trúc bão khi tiếp cận vùng đệm ven bờ.

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

Nguồn dữ liệu đầu vào phục vụ mô hình bao gồm số liệu tái phân tích và dự báo toàn cầu AVN/GFS với độ phân giải thời gian 6 giờ một lần từ Trung tâm Dự báo Môi trường Quốc gia Hoa Kỳ (NCEP). Dữ liệu địa hình sử dụng tệp GTOPO30 với 10.801 điểm theo phương ngang và 5.400 điểm theo phương dọc, độ phân giải xấp xỉ 0,3 độ kinh vĩ.

Về quy mô mẫu và phương pháp chọn mẫu, nghiên cứu áp dụng phương pháp chọn mẫu chủ đích (purposive sampling) trên 50 ca dự báo thử nghiệm từ 7 cơn bão tiêu biểu giai đoạn 2005 - 2010 (gồm Damrey, Vicente, Washi năm 2005; Durian, Xangsane năm 2006; Francisco năm 2007 và Conson năm 2010). Tiêu chí lựa chọn là các cơn bão có thời gian hoạt động trên 3 ngày và đổ bộ trực tiếp vào đất liền. Sau khi sàng lọc các trường hợp mô hình dự báo bão tan ngoài biển, tổng số 35 trường hợp dự báo đổ bộ hợp lệ được đưa vào đánh giá chi tiết (17 trường hợp không cài xoáy và 18 trường hợp có cài xoáy giả).

Miền tính mô hình được thiết lập từ vĩ độ 5 độ Nam đến 35 độ Bắc và kinh độ 90 độ Đông đến 150 độ Đông với tâm lưới đặt tại 15 độ Bắc, 120 độ Đông, độ phân giải ngang 30 km và 23 tầng thẳng đứng. Mô hình tích hợp các sơ đồ vật lý tối ưu: vi vật lý Thompson, bức xạ sóng dài RRTM, mô hình đất Noah, lớp biên MRF và tham số hóa đối lưu Betts-Miller-Janjic (BMJ). Thuật toán xác định đổ bộ tự động hóa thực hiện quét vị trí tâm bão nội suy theo bước lặp thời gian 15 phút (0,25 giờ), kiểm tra độ cao địa hình ven biển để xác định chính xác thời điểm và tọa độ tiếp bờ đầu tiên.

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

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

Quá trình thử nghiệm và phân tích định lượng trên 35 trường hợp bão đổ bộ đã mang lại các phát hiện cốt lõi:

Thứ nhất, ở hạn dự báo trước 3 ngày (72 giờ), sai số vị trí đổ bộ trung bình giữa hai phương án có cài xoáy và không cài xoáy tương đương nhau, đạt giá trị khoảng 185,4 km. Đối với phương án không cài xoáy, sai số vị trí dao động trong khoảng từ 34,8 km đến 388,0 km; phương án cài xoáy dao động từ 33,0 km đến 367,6 km. Về thời gian, phương án không cài xoáy có xu hướng dự báo bão đổ bộ muộn hơn thực tế ở 75% số trường hợp (6 trên 8 ca), với độ lệch từ -17 giờ đến +18 giờ. Ngược lại, phương án cài xoáy dự báo bão đổ bộ sớm ở 80% trường hợp (4 trên 5 ca), độ lệch từ -3 giờ đến -21 giờ.

Thứ hai, ở hạn dự báo trước 4 ngày (96 giờ), phương án cài xoáy giả thể hiện ưu thế vượt trội khi giảm sai số vị trí trung bình khoảng 120 km so với phương án không cài xoáy. Sai số vị trí lớn nhất của phương án không cài xoáy lên tới 882,3 km (ghi nhận tại bão Conson), trong khi phương án cài xoáy khống chế sai số tối đa ở mức 412,3 km, tương đương mức giảm hơn 53% sai số cực đại. Đáng chú ý, toàn bộ 100% trường hợp cài xoáy đều dự báo bão đổ bộ sớm hơn thực tế với độ lệch thời gian dao động từ 3 đến 30 giờ.

Thứ ba, ở hạn dự báo trước 5 ngày (120 giờ), xu hướng đảo chiều xuất hiện khi phương án không cài xoáy đạt độ chính xác vị trí tốt hơn phương án cài xoáy khoảng 100 km. Phương án không cài xoáy ghi nhận 60% trường hợp có sai số vị trí dưới 190 km (điển hình như bão Damrey đạt sai số chỉ 56,9 km), trong khi phương án cài xoáy ghi nhận sai số cực đại lên đến 628,5 km. Về thời gian, 80% ca cài xoáy tiếp tục dự báo bão đến sớm từ 4 đến 16 giờ.

Thứ tư, phân tích ca bão mạnh điển hình Xangsane (2006) với sức gió vùng gần tâm đạt cấp 15: Ở hạn 4 ngày, cả hai phương án cho kết quả rất sát thực tế với sai số vị trí chỉ 88,7 km (cài xoáy) và 94,6 km (không cài xoáy), thời gian đổ bộ sớm từ 3 đến 4 giờ. Ở hạn 5 ngày, sai số vị trí của hai phương án tăng lên lần lượt là 497,4 km và 483,3 km, thời gian đổ bộ sớm từ 11 đến 16 giờ.

Thảo luận kết quả

Nguyên nhân của sự khác biệt giữa hai phương án bắt nguồn từ bản chất động lực học của quá trình ban đầu hóa. Việc cài xoáy giả ở hạn 4 ngày giúp mô hình tái tạo trường xoáy ban đầu hoàn chỉnh và đối xứng hơn khi bão còn cách xa đất liền hàng nghìn kilomet, từ đó tương tác hiệu quả với dòng dẫn đường quy mô lớn. Tuy nhiên, khi thời gian tích phân kéo dài đến 120 giờ (5 ngày), các sai số tích lũy từ trường khí áp quy mô lớn và sự mất cân bằng giữa xoáy nhân tạo với trường môi trường xung quanh khiến xoáy bão di chuyển nhanh hơn thực tế, làm gia tăng sai số vị trí và dẫn đến hiện tượng bão đổ bộ sớm.

Khi đối chiếu với các trung tâm khí tượng quốc tế lớn, mức sai số vị trí hạn 5 ngày của mô hình WRF trong nghiên cứu (dao động từ 370 km đến 490 km) hoàn toàn tương đồng với sai số dự báo hạn 5 ngày của Trung tâm Dự báo Bão Hoa Kỳ NHC (khoảng 265 hải lý, tương đương 490 km) và Cơ quan Khí tượng Nhật Bản JMA (trung bình 451 km giai đoạn 2005 - 2007). Toàn bộ dữ liệu này có thể được biểu diễn trực quan thông qua các bảng thống kê sai số trực tiếp (DPE), biểu đồ sai số trung bình (MPE) phân cấp theo hạn 72, 96, 120 giờ và đồ thị phân bố độ lệch thời gian tiếp bờ.

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

Dựa trên các kết luận thực nghiệm, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược nhằm tối ưu hóa công tác dự báo bão hạn vừa tại Việt Nam:

  1. Ứng dụng quy trình cài xoáy giả có chọn lọc theo hạn dự báo: Tổng cục Khí tượng Thủy văn cần ban hành văn bản hướng dẫn nghiệp vụ chuẩn hóa, quy định áp dụng mô-đun cài xoáy giả cho các bản tin dự báo bão hạn 72 đến 96 giờ, đồng thời duy trì trường phân tích nguyên bản hoặc tích hợp đồng hóa số liệu cho hạn 120 giờ nhằm giảm sai số định vị xuống dưới 300 km, thực hiện hoàn thành trước năm 2027.

  2. Triển khai hệ thống dự báo tổ hợp quy mô vừa (Ensemble Prediction System): Trung tâm Dự báo Khí tượng Thủy văn Quốc gia phối hợp cùng các viện nghiên cứu phát triển tổ hợp dự báo từ 20 đến 30 thành phần vật lý trong giai đoạn 2026 - 2028. Giải pháp này giúp xác định xác suất vùng đổ bộ và hạn chế rủi ro dự báo lệch hướng quá 200 km ở hạn 5 ngày.

  3. Tăng cường độ phân giải lưới tính và đồng hóa dữ liệu vi ba: Nâng cấp mô hình WRF lên hệ thống lưới lồng phân giải cao 9 km và 3 km cho khu vực Biển Đông, tích hợp đồng hóa số liệu 3DVAR/4DVAR từ radar ven biển và vệ tinh địa tĩnh, đặt mục tiêu giảm sai số thời gian đổ bộ xuống dưới mức 6 giờ vào năm 2028.

  4. Mở rộng mạng lưới trạm quan trắc hải văn và thám không biển khơi: Bộ Tài nguyên và Môi trường cần đầu tư bổ sung từ 5 đến 10 trạm thám không vô tuyến tự động và hệ thống phao biển thông minh tại khu vực quần đảo Hoàng Sa và Trường Sa trước năm 2030, giúp giảm thiểu từ 15% đến 20% sai số trường gió ban đầu của các mô hình số trị.

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

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

  1. Chuyên viên dự báo khí tượng thủy văn tại Trung ương và các Đài khu vực: Tham khảo cấu hình vật lý mô hình WRF (sơ đồ Thompson, BMJ, MRF) và quy tắc hiệu chỉnh sai số sớm/muộn để vận hành nghiệp vụ dự báo bão hạn 3 đến 5 ngày, nâng cao độ tin cậy của các bản tin cảnh báo thiên tai khẩn cấp.

  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 và Hải dương học: Sử dụng toàn bộ khung lý thuyết động lực học xoáy bão, thuật toán nội suy địa hình GTOPO30 và hệ thống chỉ tiêu đánh giá MPE/ME làm tài liệu giảng dạy, nghiên cứu chuyên sâu về mô hình hóa khí quyển.

  3. Ban Chỉ đạo Quốc gia về Phòng chống thiên tai và Ủy ban Tìm kiếm cứu nạn các cấp: Ứng dụng các ngưỡng sai số vị trí (khoảng 180 đến 400 km) và xu thế bão đến sớm từ 4 đến 16 giờ để lập kịch bản ứng phó, di dời dân cư và neo đậu tàu thuyền sớm trước 96 giờ nhằm giảm thiểu tối đa thiệt hại nhân mạng.

  4. Các kỹ sư công nghệ và đơn vị phát triển giải pháp phần mềm thời tiết: Kế thừa thuật toán xác định điểm đổ bộ tự động dựa trên ma trận độ cao địa hình để xây dựng các nền tảng số hóa theo dõi đường đi bão thời gian thực phục vụ cộng đồng.

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

Mô hình WRF có xu hướng dự báo bão đổ bộ sớm hay muộn hơn thực tế?

Trong phần lớn các thử nghiệm, mô hình WRF có xu thế dự báo bão đổ bộ sớm hơn so với quan trắc thực tế. Đặc biệt ở phương án tích hợp cài xoáy giả, có tới hơn 80% trường hợp bão tiếp bờ sớm từ 3 đến 22 giờ. Điều này bắt nguồn từ việc mô phỏng gia tốc di chuyển của xoáy bão trong mô hình thường cao hơn vận tốc thực tế ngoài biển.

Khi nào nên áp dụng mô-đun cài xoáy giả trong mô hình WRF?

Kỹ thuật cài xoáy giả phát huy hiệu quả tối ưu nhất ở hạn dự báo 4 ngày (96 giờ), giúp giảm sai số vị trí trung bình khoảng 120 km và khống chế sai số cực đại dưới 412,3 km. Ngược lại, ở hạn 5 ngày (120 giờ), không nên lạm dụng cài xoáy đơn thuần vì có thể làm tăng sai số vị trí lên trên 600 km do mất cân bằng động lực kéo dài.

Thuật toán xác định vị trí và thời gian bão đổ bộ trong luận văn hoạt động ra sao?

Thuật toán kết hợp tọa độ tâm bão mô phỏng với tập dữ liệu độ cao địa hình GTOPO30 độ phân giải 0,3 độ. Bằng phương pháp nội suy chuyển động thẳng đều theo bước thời gian 15 phút giữa hai mốc dự báo 6 giờ, mô hình kiểm tra độ cao địa hình tại tâm bão nội suy; thời điểm độ cao chuyển từ giá trị âm sang dương chính là thời khắc bão tiếp bờ.

Sai số dự báo bão của mô hình WRF tại Việt Nam có tiệm cận chuẩn quốc tế không?

Hoàn toàn tiệm cận. Sai số vị trí trung bình hạn 3 ngày của mô hình WRF đạt khoảng 185,4 km, và hạn 5 ngày dao động từ 370 đến 490 km. Mức sai số này tương đương với kết quả dự báo của Trung tâm Dự báo Bão Hoa Kỳ NHC (khoảng 490 km) và Cơ quan Khí tượng Nhật Bản JMA (khoảng 451 km).

Nghiên cứu đã thử nghiệm trên những cơn bão cụ thể nào?

Nghiên cứu đã thử nghiệm trên 7 cơn bão thực tế đổ bộ vào các vùng bờ biển khác nhau của Việt Nam giai đoạn 2005 - 2010, bao gồm bão Damrey, Vicente, Washi (năm 2005), Xangsane, Durian (năm 2006), Francisco (năm 2007) và Conson (năm 2010), bảo đảm tính đại diện cao về mặt địa lý và hình thế khí tượng.

Kết luận

Công trình nghiên cứu đã khẳng định tính khả thi và độ tin cậy của việc ứng dụng mô hình số trị WRF trong dự báo bão hạn vừa trước 3 đến 5 ngày phục vụ phòng chống thiên tai tại bờ biển Việt Nam.

Năm đóng góp và kết luận khoa học quan trọng nhất của luận văn bao gồm:

  • Khẳng định mô hình WRF đạt độ chính xác cao ở hạn dự báo 3 ngày với sai số vị trí trung bình khoảng 185,4 km cho cả hai phương án thử nghiệm.
  • Chứng minh hiệu quả vượt trội của mô-đun cài xoáy giả ở hạn 4 ngày khi giúp giảm 120 km sai số vị trí trung bình và giảm 53% sai số cực đại so với phương án thông thường.
  • Làm rõ quy luật dự báo bão đổ bộ sớm ở hơn 80% trường hợp có cài xoáy với độ lệch thời gian phổ biến từ 3 đến 16 giờ.
  • Xây dựng và kiểm chứng thành công thuật toán tự động hóa xác định vị trí và thời khắc bão tiếp bờ dựa trên tệp địa hình GTOPO30 với bước lặp 15 phút.
  • Thiết lập cơ sở dữ liệu thực nghiệm chi tiết từ 35 ca bão đổ bộ làm tiền đề chuẩn hóa cấu hình mô hình WRF nghiệp vụ cho toàn ngành khí tượng thủy văn.

Trong giai đoạn 2026 - 2030, việc tiếp tục phát triển hệ thống dự báo tổ hợp và tích hợp đồng hóa dữ liệu vệ tinh đa phổ sẽ là bước đi then chốt nhằm nâng cao hơn nữa độ chính xác của các bản tin cảnh báo bão sớm. Các cơ quan quản lý, nhà nghiên cứu và đơn vị tác nghiệp cần khai thác ngay các kết quả này để hoàn thiện quy trình ứng phó rủi ro thiên tai, bảo vệ an toàn tính mạng và tài sản cho cộng đồng ven biển.