Tổng quan nghiên cứu

Khu vực Biển Đông mỗi năm đón nhận trung bình từ 10 đến 12 cơn bão và áp thấp nhiệt đới, trong đó khoảng 5 đến 7 cơn bão ảnh hưởng trực tiếp đến vùng biển và đất liền Việt Nam, gây thiệt hại nghiêm trọng về người và tài sản cho cộng đồng cư dân ven biển cũng như các hoạt động hàng hải. Mặc dù công nghệ dự báo khí tượng toàn cầu đã có những bước tiến vượt bậc, việc mô phỏng và dự báo chính xác quỹ đạo di chuyển của các xoáy thuận nhiệt đới trên Biển Đông vẫn đối mặt với nhiều thách thức lớn do tương tác đa quy mô phức tạp giữa hoàn lưu khí quyển, địa hình bán đảo Đông Dương và bề mặt đại dương.

Mục tiêu cốt lõi của nghiên cứu là khảo sát, đánh giá định lượng tác động của việc thay đổi độ phân giải không gian trong mô hình động lực khu vực WRF (Weather Research and Forecasting) phiên bản 3.1 đối với độ chính xác dự báo quỹ đạo bão trên Biển Đông. Đề tài tập trung so sánh hiệu năng giữa cấu hình lưới đơn độ phân giải thô 27 km và cấu hình lưới lồng tương tác hai chiều với độ phân giải mịn 9 km.

Nghiên cứu giới hạn phạm vi không gian bao phủ toàn bộ vùng Biển Đông (từ 5°S đến 35°N và 90°E đến 130°E đối với miền thô; từ 10°N đến 20°N và 106.5°E đến 116°E đối với miền lồng), với thời hạn dự báo chi tiết từ 00 giờ đến 72 giờ cho các cơn bão hoạt động trong giai đoạn 2006 - 2013. Kết quả của luận văn mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở định lượng để tối ưu hóa tài nguyên tính toán hiệu năng cao, giảm sai số vị trí tâm bão ở các hạn dự báo quan trọng 24 giờ và 48 giờ xuống từ 15% đến 25%, hỗ trợ đắc lực cho công tác cảnh báo thiên tai và bảo đảm an toàn hàng hải.

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 hệ phương trình nguyên thủy phi thủy tĩnh nén được của mô hình khí quyển WRF-ARW, tích hợp hệ tọa độ thẳng đứng theo khối lượng thủy tĩnh (tọa độ $\eta$) với 27 mực thẳng đứng, cho phép mô phỏng chính xác các quá trình động lực học và nhiệt động học trên địa hình phức tạp. Mô hình áp dụng sơ đồ tích phân thời gian Runge-Kutta bậc 3 nhằm duy trì tính ổn định tính toán khi giải các sóng âm và sóng trọng trường trong khí quyển.

Về mặt động lực học xoáy, nghiên cứu vận dụng lý thuyết dòng dẫn đường (Steering Flow Theory) dựa trên trường độ dày trung bình lớp sâu (Deep Layer Mean - DLM) và tương tác phi tuyến giữa hoàn lưu xoáy bão với trường môi trường quy mô synop. Bão di chuyển chịu sự chi phối chủ yếu của dòng dẫn đường ở rìa áp cao cận nhiệt đới Tây Thái Bình Dương kết hợp với nội lực sinh ra do hiệu ứng beta và sự phân tán sóng Rossby.

Kỹ thuật lưới lồng (Nested Grid) trên mạng lưới so le Arakawa-C với tỷ lệ bước lưới 3:1 (27 km xuống 9 km) được áp dụng để liên kết tương tác hai chiều giữa miền tính lớn và miền tính chi tiết, đảm bảo truyền tải đầy đủ các nhiễu động vật lý quy mô vừa và nhỏ vào cấu trúc hoàn lưu bão.

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

Nguồn dữ liệu đầu vào sử dụng số liệu phân tích và dự báo tái phân tích toàn cầu GFS (Global Forecast System) của Trung tâm Quốc gia về Dự báo Môi trường Hoa Kỳ (NCEP) với độ phân giải không gian $1^\circ \times 1^\circ$ cập nhật định kỳ 6 giờ một lần. Dữ liệu đối chứng kiểm nghiệm là tập số liệu quan trắc quỹ đạo thực tế (Best Track) từ Trung tâm Cảnh báo Bão Liên hợp (JTWC) và Cơ quan Khí tượng Nhật Bản (JMA).

Mẫu nghiên cứu gồm 15 trường hợp bão điển hình hoạt động trên Biển Đông trong giai đoạn 2006 - 2013, được phân loại theo hai tiêu chí: cường độ (bão mạnh có vận tốc gió cực đại trên 32.5 m/s và bão yếu dưới 32.5 m/s) và thời gian hoạt động (bão đầu mùa từ tháng 5 đến tháng 8, bão cuối mùa từ tháng 9 đến tháng 12). Phương pháp chọn mẫu phân tầng có chủ đích này phản ánh đầy đủ các dạng quỹ đạo từ đơn giản đến phức tạp (như bão Chebi và bão Prapiroon năm 2006).

Quy trình phân tích số trị kết hợp các gói tham số hóa vật lý: vi vật lý WSM3, bức xạ sóng dài RRTM, bức xạ sóng ngắn Dudhia, lớp biên hành tinh YSU và đối lưu mây tích Kain-Fritsch. Thuật toán Downhill Simplex tối ưu hóa đa chiều được lập trình bằng Fortran 90 trong công cụ hậu xử lý để tự động định vị tâm áp suất cực tiểu và độ xoáy cực đại. Sai số khoảng cách giữa tâm bão dự báo và thực tế được tính toán thông qua công thức lượng giác mặt cầu (Great-Circle Distance), phân tách thành sai số dọc quỹ đạo (ATE) và sai số ngang quỹ đạo (CTE).

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

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

Kết quả thử nghiệm số trị cho thấy việc tăng độ phân giải không gian từ 27 km xuống 9 km thông qua kỹ thuật lưới lồng tương tác hai chiều mang lại sự cải thiện rõ rệt trong dự báo quỹ đạo bão ở tất cả các hạn thời gian từ 12 giờ đến 72 giờ:

  • Sai số vị trí trung bình toàn bộ tập mẫu ở hạn dự báo 24 giờ giảm từ 142.6 km ở lưới thô 27 km xuống còn 118.4 km khi sử dụng lưới lồng 9 km (cải thiện tương đương 16.97%). Ở hạn dự báo 48 giờ, sai số trung bình giảm từ 268.5 km xuống còn 214.2 km (mức cải thiện đạt 20.22%).
  • Đối với nhóm bão mạnh (vận tốc gió cực đại trên 32.5 m/s, ví dụ bão Chebi năm 2006), lưới lồng 9 km mô phỏng chính xác mắt bão và cấu trúc thành mắt bão, giúp sai số quỹ đạo ở hạn 48 giờ giảm mạnh từ 245.0 km xuống 182.3 km (giảm 25.59%).
  • Đối với nhóm bão có quỹ đạo ngoặt và đổi hướng phức tạp (ví dụ bão Prapiroon năm 2006), mô hình lưới lồng 9 km bắt kịp thời điểm chuyển hướng từ tây-tây bắc sang tây bắc sớm hơn 12 giờ so với lưới đơn 27 km, hạn chế hiện tượng bão trôi dạt quá nhanh về phía tây.
  • Phân tích sai số hệ thống cho thấy lưới thô 27 km có xu hướng dự báo bão di chuyển lệch về phía bắc và nhanh hơn thực tế trung bình khoảng 8.5% về tốc độ, trong khi lưới lồng 9 km hiệu chỉnh độ lệch phương vị, đưa sai số hướng về dưới 12 độ ở hạn 24 giờ.

Thảo luận kết quả

Sự vượt trội của cấu hình lưới lồng 9 km bắt nguồn từ khả năng giải quyết các quá trình động lực học quy mô trung gian (mesoscale) mà lưới 27 km thường làm nhẵn hoặc bỏ qua. Ở độ phân giải 9 km, građien khí áp ngang tại vùng lõi xoáy được tái tạo dốc hơn, cho phép mô hình tính toán chuẩn xác sự giải phóng ẩn nhiệt ngưng kết trong các dải mây đối lưu sâu quanh tâm bão. Điều này củng cố cấu trúc hoàn lưu xoáy thẳng đứng, giúp tâm xoáy bão duy trì liên kết chặt chẽ với dòng dẫn đường ở mực 500 hPa và 700 hPa.

+-------------------------------------------------------------------------+
|      BẢNG TỔNG HỢP SAI SỐ VỊ TRÍ TRUNG BÌNH (ĐƠN VỊ: KM)               |
+-------------------+----------+----------+----------+----------+---------+
| Hạn dự báo        | 12 giờ   | 24 giờ   | 36 giờ   | 48 giờ   | 72 giờ  |
+-------------------+----------+----------+----------+----------+---------+
| Lưới thô (27 km)  | 78.5     | 142.6    | 205.3    | 268.5    | 412.0   |
| Lưới lồng (9 km)  | 62.1     | 118.4    | 164.8    | 214.2    | 345.6   |
| Mức cải thiện (%) | 20.89%   | 16.97%   | 19.73%   | 20.22%   | 16.12%  |
+-------------------+----------+----------+----------+----------+---------+

Trên các bản đồ phân tích trường khí áp mặt biển (SLP) và địa thế vị mực 500 hPa trích xuất từ GrADs, miền lưới lồng 9 km thể hiện ranh giới áp cao cận nhiệt đới rõ nét, tránh hiện tượng xâm lấn giả của áp cao vào tâm bão. Khi đối chiếu với các công trình quốc tế, mức sai số 118.4 km ở hạn 24 giờ của cấu hình WRF 9 km hoàn toàn cạnh tranh với các trung tâm dự báo hàng đầu như JMA (khoảng 110 - 130 km) và vượt trội hơn đáng kể so với các mô hình chính áp truyền thống như WBAR hay LBAR (sai số 24 giờ thường dao động từ 150 km đến 180 km).

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

  • Ứng dụng cấu hình lưới lồng hai chiều 27 km - 9 km vào hệ thống dự báo nghiệp vụ hàng ngày: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần chuẩn hóa kịch bản chạy WRF với miền lồng bao trùm kinh độ 105°E - 120°E và vĩ độ 10°N - 22°N, hướng tới mục tiêu giảm sai số dự báo vị trí bão 24 giờ xuống dưới 100 km trong giai đoạn 2024 - 2026.
  • Đầu tư và nâng cấp hạ tầng siêu máy tính tính toán song song (HPC): Các cơ quan nghiên cứu và quản lý dự báo khí tượng thủy văn cần trang bị hệ thống cụm máy chủ hỗ trợ thư viện MPI với tối thiểu 128 đến 256 lõi xử lý, bảo đảm thời gian tích phân mô hình WRF lưới lồng cho hạn 72 giờ hoàn thành trong vòng 45 phút kể từ khi tiếp nhận dữ liệu biên toàn cầu.
  • Tích hợp kỹ thuật đồng hóa số liệu địa phương (WRF-Data Assimilation): Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu chủ trì triển khai hệ thống 3DVAR/EnKF, tích hợp dữ liệu trạm thám không vô tuyến ven biển, radar thời tiết thế hệ mới và số liệu gió vệ tinh scatt, nhằm cải thiện trường ban đầu trước mỗi mùa bão từ năm 2025.
  • Thiết lập quy trình hiệu chỉnh sai số thống kê hậu dự báo (Model Output Statistics - MOS): Phòng Bảo đảm Hàng hải và Cục Hàng hải Việt Nam phối hợp xây dựng module tự động nắn chỉnh quỹ đạo dự báo dựa trên cơ sở dữ liệu sai số lịch sử, nâng cao độ tin cậy của bản tin thời tiết chuyên dụng phục vụ cứu hộ cứu nạn trên biển.

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

  • Dự báo viên khí tượng tại các đài khí tượng thủy văn khu vực: Sử dụng luận văn như cẩm nang kỹ thuật thực hành thiết lập file cấu hình namelist.input, tối ưu hóa các tham số miền tính và vận hành mô hình WRF phục vụ tác nghiệp dự báo mùa bão.
  • Sĩ quan tham mưu tác chiến Hải quân, Cảnh sát biển và lực lượng Kiểm ngư: Khai thác phương pháp phân tích sai số ATE và CTE để xây dựng các kịch bản phòng tránh bão an toàn, xác định hành lang cơ động và neo đậu tàu thuyền cách xa tâm bão tối thiểu 150 - 200 hải lý.
  • Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Khí tượng học: Sử dụng tài liệu làm cơ sở tham khảo chuyên sâu về phương pháp số trị giải phương trình vi phân khí quyển, kỹ thuật ghép lưới lồng và giải thuật Downhill định vị cực trị khí áp.
  • Kỹ sư phát triển phần mềm mô phỏng môi trường và dữ liệu không gian GIS: Ứng dụng quy trình chuyển đổi định dạng dữ liệu NetCDF sang GrADs để phát triển các cổng thông tin trực quan hóa trường gió, sóng và quỹ đạo thiên tai thời gian thực.

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

Mô hình WRF có ưu điểm gì vượt trội so với các mô hình thế hệ cũ như MM5 hay HRM?

WRF được thiết kế theo cấu trúc mã nguồn mở module hóa hiện đại, sử dụng hệ phương trình phi thủy tĩnh nén hoàn toàn với tọa độ khối lượng $\eta$. Khác với MM5 hay HRM bị giới hạn ở một số sơ đồ tham số hóa cố định, WRF cung cấp hàng chục tùy chọn vật lý vi mô, bức xạ và lớp biên tiên tiến, đồng thời hỗ trợ tối ưu hóa tính toán song song quy mô lớn qua MPI, giúp xử lý các hiện tượng quy mô vừa với độ chính xác cao hơn.

Tại sao việc tăng độ phân giải từ 27 km xuống 9 km lại làm giảm sai số quỹ đạo bão?

Việc thu hẹp khoảng cách lưới từ 27 km xuống 9 km giúp mô hình mô tả chân thực cấu trúc dòng thăng, dòng giáng thẳng đứng và gradient khí áp ngang quanh mắt bão. Độ phân giải cao giải quyết tốt hơn sự bất đối xứng của trường gió và quá trình giải phóng ẩn nhiệt ngưng kết trong mây đối lưu, từ đó mô phỏng chính xác lực nội tại và tương tác của bão với dòng dẫn đường quy mô lớn.

Kỹ thuật lưới lồng tương tác hai chiều (2-way nesting) hoạt động như thế nào?

Trong cơ chế 2-way nesting, miền lưới thô (27 km) cung cấp các giá trị điều kiện biên bên ngoài cho miền lưới mịn (9 km) tại mỗi bước thời gian. Ngược lại, sau khi miền lưới mịn hoàn thành tích phân các quá trình vật lý chi tiết, các giá trị trạng thái khí tượng tại các điểm lưới trùng khớp sẽ được phản hồi và ghi đè cập nhật lại cho miền lưới thô, đảm bảo tính cân bằng động lực học toàn hệ thống.

Thuật toán Downhill Simplex thực hiện nhiệm vụ gì trong công tác hậu xử lý?

Thuật toán Downhill là phương pháp tìm cực trị hàm đa chiều không cần sử dụng đạo hàm. Trong nghiên cứu này, thuật toán so sánh giá trị khí áp và độ xoáy tại ba đỉnh của tam giác biến dạng trên trường lưới mô phỏng để tự động hội tụ chính xác về vị trí tâm bão có áp suất thấp nhất, giúp loại bỏ hoàn toàn các sai số chủ quan so với việc quan sát bằng mắt thường trên bản đồ synop.

Cần lưu ý những hạn chế nào khi vận hành mô hình WRF độ phân giải cao cho Biển Đông?

Hạn chế lớn nhất là chi phí tài nguyên tính toán tăng theo lũy thừa bậc ba khi giảm bước lưới, đòi hỏi hệ thống máy chủ cấu hình mạnh để bảo đảm thời gian phát bản tin dự báo. Bên cạnh đó, chất lượng dự báo vẫn phụ thuộc chặt chẽ vào độ chính xác của trường phân tích ban đầu từ các mô hình toàn cầu như GFS; nếu trường ban đầu thiếu dữ liệu quan trắc biển, mô hình có thể phát sinh sai lệch pha ban đầu.

Kết luận

  • Luận văn đã thiết lập thành công mô hình WRF 3.1 với cấu hình lưới lồng tương tác hai chiều 27 km - 9 km, chứng minh tính khả thi và độ tin cậy cao trong mô phỏng bão nhiệt đới trên khu vực Biển Đông.
  • Việc tăng độ phân giải ngang lên 9 km giúp giảm sai số vị trí tâm bão trung bình khoảng 24.2 km ở hạn 24 giờ và 54.3 km ở hạn 48 giờ so với lưới thô 27 km.
  • Mô hình thể hiện năng lực vượt trội đối với các cơn bão có cường độ mạnh và cấu trúc quỹ đạo phức tạp, bám sát các diễn biến ngoặt dòng do tác động của áp cao cận nhiệt đới.
  • Luận văn đóng góp quy trình công nghệ hoàn chỉnh từ tiền xử lý WPS, tích phân ARW đến hậu xử lý tự động định vị tâm bão bằng giải thuật Downhill viết trên Fortran 90.
  • Định hướng giai đoạn 2025 - 2027 cần mở rộng thử nghiệm lưới lồng di động 3 km kết hợp đồng hóa dữ liệu radar và thám không biển để hoàn thiện hệ thống cảnh báo sớm bão phục vụ quốc phòng và dân sinh. Hãy tham khảo toàn văn công trình nghiên cứu để nắm bắt chi tiết các bộ tham số cài đặt và bảng số liệu kiểm nghiệm thực tế.