Tổng quan nghiên cứu

Khu vực Tây Bắc Thái Bình Dương và Biển Đông là một trong những cái nôi bão nhiệt đới lớn nhất thế giới, chiếm khoảng 30% tổng số cơn bão toàn cầu. Hàng năm, khu vực này ghi nhận trung bình từ 25 đến 30 cơn xoáy thuận nhiệt đới, trong đó khoảng 70% các cơn bão có áp suất tại tâm đạt từ 980 đến 1000 mb, và khoảng 10% phát triển thành siêu bão với áp suất giảm sâu xuống dưới 930 mb (tiêu biểu như cơn bão Tip năm 1979 đạt mức kỷ lục 870 mb).

Trường gió và trường khí áp là hai thành phần dữ liệu đầu vào mang tính quyết định đến độ chính xác của các mô hình thủy động lực học tính sóng và nước dâng do bão như SLOSH, Delft3D hay MIKE. Tuy nhiên, các mô hình trường gió truyền thống thường gặp khó khăn trong việc giải phương trình vi phân góc lệch hướng gió, dễ dẫn đến hiện tượng dao động nghiệm không ổn định ở vùng sát tâm bão.

Mục tiêu trọng tâm của nghiên cứu là xây dựng một mô hình số đơn giản, giải trực tiếp trường gió và trường áp từ các thông số bão dựa trên nền tảng mô hình SLOSH của Jelesnianski (1973), đồng thời xác định mối quan hệ định lượng giữa độ giảm áp tâm bão và vận tốc gió cực đại. Phạm vi nghiên cứu bao quát toàn bộ vùng biển Tây Bắc Thái Bình Dương và Biển Đông trong dải vĩ độ từ 5°B đến 45°B, sử dụng chuỗi số liệu lưu trữ bão lịch sử kéo dài nhiều thập kỷ. Kết quả nghiên cứu mang ý nghĩa thực tiễn lớn, giúp giảm thiểu sai số tính toán trường gió bão xuống dưới 5%, nâng cao độ tin cậy trong công tác cảnh báo thiên tai và thiết kế công trình cảng - biển.

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 xoáy thuận nhiệt đới và lý thuyết cân bằng gió gradient. Chuyển động của khối không khí trong bão chịu sự chi phối của ba lực cơ bản: lực gradient áp suất, lực ly tâm và lực Coriolis. Khi bỏ qua lực ma sát và lực Coriolis ở vùng gần tâm, dòng khí đạt trạng thái cân bằng cyclostrophic, cho phép thiết lập quan hệ giải tích giữa độ giảm áp tâm bão và vận tốc gió cực đại dưới dạng hàm mũ $W_m = C(p_n - p_c)^n$.

Khung lý thuyết của luận văn tích hợp mô hình bão tham số SLOSH của Jelesnianski (1973) cùng các mô hình phân bố khí áp và vận tốc gió kinh điển của Schloemer (1954), Holland (1980), Fujita (1952) và Takahashi (1939). Các khái niệm cốt lõi được áp dụng xuyên suốt bao gồm: vận tốc gió cực đại ($W_m$), bán kính đến vận tốc gió cực đại ($R_w$), kích thước cơn bão ($R_s$), góc lệch hướng tâm của vector gió ($\phi$), và áp suất ngoại biên khí quyển chuẩn ($p_n = 1010$ mb).

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

Nguồn dữ liệu chính của đề tài là tập dữ liệu bão tốt nhất (Best Track) được cung cấp bởi Trung tâm Khí tượng Chuyên ngành Bão Tokyo (RSMC Tokyo). Cỡ mẫu nghiên cứu bao gồm hàng trăm cơn bão thực tế hoạt động trên vùng Tây Bắc Thái Bình Dương và Biển Đông. Phương pháp chọn mẫu phân tầng theo không gian địa lý được áp dụng triệt để, chia toàn bộ dữ liệu thành các dải vĩ độ hẹp 5 độ (từ 5°B đến 45°B) và dải rộng 10 độ vĩ nhằm khảo sát chi tiết ảnh hưởng của lực Coriolis.

Phương pháp phân tích trọng tâm là phương pháp sai phân hữu hạn hiện kết hợp thuật toán lặp trực tiếp đại số. Bằng cách bỏ qua số hạng đạo hàm bậc cao theo bán kính của góc lệch hướng gió, mô hình cho phép tìm lời giải trực tiếp tại từng điểm tính toán mà không cần tích phân toàn miền bằng phương pháp Runge-Kutta bậc 4 như trước đây. Lý do lựa chọn giải pháp này là nhằm triệt tiêu hoàn toàn hiện tượng dao động toán học gần tâm bão và tối ưu hóa thời gian tính toán. Toàn bộ quá trình nghiên cứu và xử lý số liệu được thực hiện trong thời gian 12 tháng, từ tháng 09/2011 đến tháng 09/2012.

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

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

Thứ nhất, mô hình trường gió hiệu chỉnh đã giải quyết triệt để sự mất ổn định số trị của mô hình SLOSH gốc. Lời giải lặp trực tiếp giúp đường cong phân bố góc lệch hướng gió và trường khí áp biến thiên trơn tru từ biên ngoài $R_s = 410$ km vào tận tâm bão $r = 0$, triệt tiêu hoàn toàn các bước nhảy dao động sai lệch ở vùng bán kính dưới 47 km.

Thứ hai, kết quả khảo sát bán kính gió cực đại ($R_w$) khẳng định tính hiệu quả của hai phương án tính toán: sử dụng giá trị trung bình $R_w = 47$ km (theo kinh nghiệm của Hsu và Yan) và công thức hàm phụ thuộc vĩ độ của H. Rahn. Cả hai phương án đều cho trường áp suất mô phỏng bám rất sát số liệu thực tế, với sai số tương đối về áp suất tâm nhỏ hơn 5%.

Thứ ba, nghiên cứu đã lượng hóa thành công mối quan hệ giữa độ giảm áp tâm bão $(p_n - p_c)$ và vận tốc gió cực đại $W_m$ cho từng đới 5 độ vĩ độ. Hệ số kinh nghiệm $C$ dao động trong khoảng từ 3.5 đến 7.5 và số mũ $n$ biến thiên từ 0.45 đến 0.65 tùy thuộc vào vị trí địa lý.

Thứ tư, tại khu vực Biển Đông, vận tốc gió cực đại tương ứng với cùng một mức áp suất tâm có giá trị thấp hơn khoảng 5% đến 8% so với toàn vùng Tây Bắc Thái Bình Dương. Hiện tượng này bắt nguồn từ sự tiêu hao năng lượng do ma sát địa hình phức tạp của các quần đảo và đường bờ biển bao bọc.

Thảo luận kết quả

Dữ liệu nghiên cứu được biểu diễn trực quan qua các hệ thống đồ thị tương quan phân tán và đường hồi quy phi tuyến theo từng dải 5 độ vĩ độ, kết hợp cùng bảng thống kê chi tiết các tham số thực nghiệm. Khi so sánh với các công thức kinh điển trên thế giới, kết quả của luận văn thể hiện sự tương đồng cao với nghiên cứu của Atkinson và Holliday (1977, $C = 3.644$) ở các vĩ độ thấp, đồng thời tiệm cận với mô hình của Takahashi (1939, $C = 6.5$) tại các vĩ độ cao.

Nguyên nhân chính dẫn đến sự thay đổi của hệ số $C$ và số mũ $n$ theo vĩ độ là do gia tăng thông số Coriolis $f = 2\Omega \sin\varphi$. Khi bão di chuyển từ vùng nhiệt đới (khoảng 10°B - 15°B) lên các vĩ độ cao hơn (30°B - 40°B), lực Coriolis lớn hơn làm gia tăng độ lệch dòng khí và biến đổi cấu trúc động lực học của bão. Việc phân chia chi tiết theo đới 5 độ vĩ giúp mô hình số thể hiện quy luật vật lý rõ ràng và chính xác hơn nhiều so với việc áp dụng một công thức đồng nhất cho toàn đại dương.

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

Thứ nhất, tích hợp ngay thuật toán trường gió hiệu chỉnh vào hệ thống phần mềm dự báo sóng và nước dâng quốc gia. Viện Khoa học Khí tượng Thủy văn và Hải văn cùng các trung tâm dự báo cần cập nhật module này trong vòng 6 đến 12 tháng tới, nhằm nâng cao độ chính xác dự báo độ cao sóng và mực nước dâng do bão lên thêm 15% đến 20%.

Thứ hai, chuẩn hóa và áp dụng bộ công thức quan hệ áp suất - vận tốc gió theo từng đới vĩ độ vào công tác tính toán tải trọng công trình biển. Bộ Xây dựng và Cục Hàng hải Việt Nam nên rà soát, đưa các hệ số thực nghiệm này vào quy chuẩn thiết kế đê kè, bến cảng và công trình ngoài khơi trong lộ trình 2024-2026, giúp tối ưu hóa chi phí đầu tư xây dựng khoảng 10% mà vẫn đảm bảo hệ số an toàn.

Thứ ba, xây dựng cơ sở dữ liệu động số hóa 100% các tham số bão trên Biển Đông. Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần chủ trì số hóa toàn bộ dữ liệu lịch sử bão, phân loại chi tiết các hệ số $C$ và $n$ đặc thù cho từng phân vùng biển Việt Nam, hoàn thành trong thời gian 18 tháng.

Thứ tư, đẩy mạnh ứng dụng công nghệ viễn thám vệ tinh kết hợp mô hình số. Các trường đại học và viện chuyên ngành cần ứng dụng dữ liệu ảnh vệ tinh đa phổ để xác định nhanh bán kính gió cực đại $R_w$ theo thời gian thực, rút ngắn chu kỳ cập nhật dữ liệu đầu vào cho mô hình gió xuống mức 3 giờ một lần.

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

Nhóm kỹ sư thiết kế công trình thủy công và dầu khí ngoài khơi: Luận văn cung cấp phương pháp xác định trường áp và trường gió chính xác, giúp tính toán chính xác tải trọng sóng tác động lên thân giàn khoan, đê chắn sóng bến cảng với chu kỳ lặp thiết kế 50 đến 100 năm.

Nhóm chuyên viên dự báo khí tượng thủy văn hàng hải: Dễ dàng ứng dụng các công thức tương quan thực nghiệm để ước tính nhanh vận tốc gió bão cực đại từ số liệu áp suất tâm quan trắc qua thám không hoặc vệ tinh, phục vụ cảnh báo sớm cho tàu thuyền trên biển.

Nhóm nhà nghiên cứu và giảng viên chuyên ngành Kỹ thuật Biển: Tài liệu là tài khảo cứu giá trị về phương pháp số, kỹ thuật xử lý bất ổn định toán học trong mô hình SLOSH và phương pháp hồi quy tham số khí tượng động lực học.

Nhóm cơ quan quản lý nhà nước về phòng chống thiên tai ven biển: Có cơ sở khoa học vững chắc để xây dựng các kịch bản ngập lụt, quy hoạch không gian bờ biển và bảo vệ an toàn cho hàng triệu cư dân thuộc 28 tỉnh, thành phố ven biển Việt Nam.

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

Mô hình gió bão trong luận văn có điểm gì cải tiến nổi bật so với mô hình SLOSH gốc của Jelesnianski? Mô hình gốc sử dụng phương pháp Runge-Kutta bậc 4 giải phương trình vi phân góc lệch hướng gió nên thường bị dao động số ở vùng sát tâm bão. Nghiên cứu đã loại bỏ số hạng đạo hàm bậc cao, chuyển thành bài toán lặp đại số trực tiếp từng điểm. Giải pháp này giúp trường gió trơn tru tuyệt đối và rút ngắn hơn 30% thời gian tính toán.

Tại sao phải phân chia dữ liệu bão thành từng đới 5 độ vĩ độ để phân tích? Lực Coriolis phụ thuộc trực tiếp vào sin của vĩ độ địa lý. Khi bão di chuyển từ vĩ độ 5°B lên 45°B ở Tây Bắc Thái Bình Dương, lực Coriolis tăng đáng kể làm thay đổi trạng thái cân bằng gió gradient. Phân chia theo đới 5 độ giúp xác định chính xác sự biến thiên của hệ số $C$ và số mũ $n$, tránh sai số khi áp dụng đại trà.

Việc giả định bán kính gió cực đại trung bình bằng 47 km có gây sai số lớn không? Theo các nghiên cứu thực nghiệm của Hsu và Yan, giá trị trung bình $R_w = 47$ km đại diện rất tốt cho các cơn bão ở Tây Bắc Thái Bình Dương. Kiểm chứng trong luận văn cho thấy sai số giữa việc cố định $R_w = 47$ km và tính theo công thức động của H. Rahn là không đáng kể, đều cho độ lệch áp suất tâm dưới 5%.

Mối quan hệ giữa áp suất tâm và vận tốc gió tại Biển Đông có gì khác biệt so với ngoài khơi Thái Bình Dương? Do Biển Đông là vùng biển nửa kín, được bao bọc bởi nhiều khối lục địa và đảo lớn, lực ma sát bề mặt làm suy giảm động năng của bão. Ở cùng một mức độ giảm áp tâm bão, vận tốc gió cực đại ghi nhận trên Biển Đông thường thấp hơn từ 5% đến 8% so với khu vực ngoài khơi Tây Bắc Thái Bình Dương.

Mô hình trường gió hiệu chỉnh này cần những thông số đầu vào tối thiểu nào để hoạt động? Mô hình chỉ cần các thông số bão cơ bản bao gồm: vận tốc gió cực đại (hoặc áp suất tâm), tọa độ vĩ độ tâm bão, bán kính gió cực đại (mặc định 47 km), kích thước cơn bão ($R_s = 410$ km) và áp suất khí quyển ngoại biên ($p_n = 1010$ mb).

Kết luận

  • Xây dựng thành công thuật toán lặp trực tiếp giải góc lệch hướng gió bão, khắc phục hoàn toàn hiện tượng dao động toán học gần tâm bão trong mô hình SLOSH.
  • Thiết lập hệ thống phương trình thực nghiệm liên hệ giữa độ giảm áp tâm bão và vận tốc gió cực đại cho toàn bộ dải vĩ độ từ 5°B đến 45°B khu vực Tây Bắc Thái Bình Dương.
  • Chứng minh rõ nét quy luật biến thiên có hệ thống của hệ số $C$ và số mũ $n$ theo vĩ độ địa lý dưới tác động của lực Coriolis.
  • Xác định sự khác biệt về đặc trưng động lực học trường gió trên Biển Đông so với vùng biển mở Thái Bình Dương với mức chênh lệch vận tốc gió từ 5% đến 8%.
  • Đóng góp nguồn dữ liệu và công cụ tính toán trường gió tin cậy cho các bài toán mô phỏng sóng, nước dâng và thiết kế công trình cảng - bờ biển tại Việt Nam.

Công trình nghiên cứu mở ra giải pháp đột phá trong tính toán khí tượng hải văn công trình. Các đơn vị tư vấn thiết kế, cơ quan khí tượng và doanh nghiệp hàng hải hãy ứng dụng ngay mô hình này vào thực tiễn để nâng cao chất lượng dự báo và đảm bảo an toàn tuyệt đối cho các công trình ven biển.