Tổng quan luận án
Luận án tiến sĩ "Nghiên cứu phát triển và ứng dụng sơ đồ ban đầu hóa xoáy ba chiều cho mục đích dự báo chuyển động bão ở Việt Nam" do nghiên cứu sinh Bùi Hoàng Hải thực hiện năm 2008 tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội, thuộc chuyên ngành Khí tượng học (mã số 62 44 87 01), dưới sự hướng dẫn của PGS. Phan Văn Tân và PGS. Nguyễn Đăng Quế.
Tính cấp thiết và khoảng trống nghiên cứu
Bão nhiệt đới là một trong những hiện tượng thời tiết nguy hiểm gây thiệt hại nặng nề thông qua gió mạnh, mưa lớn, lũ lụt và nước dâng. Bối cảnh biến đổi khí hậu toàn cầu kèm theo xu thế gia tăng cường độ và tính tàn phá của bão (Emanuel, 2005) đặt ra yêu cầu nâng cao năng lực dự báo nghiệp vụ. Tại thời điểm nghiên cứu, quy trình dự báo bão ở Việt Nam chủ yếu dựa vào phương pháp synop và mô hình thống kê CLIPER (Climatology and Persistence — mô hình không kỹ năng làm mốc so sánh), với hạn dự báo nghiệp vụ chỉ dừng ở 24 giờ. Các mô hình số trị khu vực tuy đã bắt đầu được thử nghiệm nhưng độ chính xác còn hạn chế.
Mô hình chính áp WBAR (Weber, 2001) tuy có kỹ năng nhất định đối với các cơn bão mạnh có quỹ đạo ổn định nhưng lại bộc lộ sai số lớn khi bão chịu sự chi phối của các hệ thống thời tiết tà áp phức tạp (Lê Công Thành, 2004). Trong khi đó, mô hình khu vực phân giải cao HRM (Kiều Thị Xin và cộng sự, 2002) có cấu trúc khí quyển ba chiều đầy đủ nhưng vẫn gặp sai số quỹ đạo đáng kể, nhất là với các cơn bão cường độ yếu. Nguyên nhân kỹ thuật cốt lõi là mạng lưới trạm quan trắc synop truyền thống trên các vùng biển nhiệt đới rất thưa thớt, khiến trường phân tích quy mô lớn toàn cầu không biểu diễn chính xác vị trí tâm, cấu trúc và cường độ của xoáy bão thực tế. Khoảng trống nghiên cứu được luận án xác định là việc thiếu một sơ đồ ban đầu hóa xoáy ba chiều có cơ sở động lực cân bằng vững chắc, phù hợp để tích hợp vào mô hình số trị tà áp nhằm nâng hạn dự báo quỹ đạo bão tại Việt Nam lên 48 giờ.
Mục tiêu nghiên cứu
Mục đích tổng quát của luận án là xây dựng một sơ đồ ban đầu hóa xoáy ba chiều và áp dụng cho một mô hình số trị nhằm góp phần nâng cao chất lượng dự báo quỹ đạo bão ở Việt Nam tới hạn 48 giờ.
Các nhiệm vụ cụ thể được văn bản xác định gồm:
- Nghiên cứu, phát triển sơ đồ ban đầu hóa xoáy ba chiều mới dựa trên lý thuyết xoáy cân bằng, có thể áp dụng cho mô hình dự báo với số liệu thực hoặc nghiên cứu lý tưởng.
- Áp dụng sơ đồ ban đầu hóa xoáy mới vào mô hình khu vực phân giải cao HRM và phát triển mô hình này thành phiên bản chuyên dụng HRM_TC (vừa có chức năng dự báo bão vừa dự báo thời tiết tổng quát).
- Khảo sát thực nghiệm độ nhạy để xác lập bộ tham số khả dụng tối ưu cho sơ đồ ban đầu hóa xoáy của HRM_TC đối với các cơn bão hoạt động trên khu vực Biển Đông.
- Thử nghiệm dự báo quỹ đạo trên các trường hợp bão thực tế bằng mô hình HRM_TC để kiểm chứng định lượng hiệu quả của sơ đồ.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Cấu trúc động lực ba chiều của xoáy thuận nhiệt đới và chuyển động quỹ đạo bão trong mô hình số trị.
- Phạm vi không gian: Khu vực Biển Đông và Tây Bắc Thái Bình Dương, tập trung vào các cơn bão có ảnh hưởng trực tiếp hoặc gián tiếp đến vùng biển và đất liền Việt Nam. Miền tính toán của mô hình HRM bao phủ toàn bộ Biển Đông và các vùng phụ cận.
- Phạm vi thời gian và tư liệu: Dữ liệu trường phân tích toàn cầu GME của Cơ quan Khí tượng Đức (DWD); số liệu quan trắc chỉ thị bão thực tế trong giai đoạn 2003–2006 (khảo sát chi tiết các cơn bão: Koni 2003, Imbudo 2003, Conson 2004, Chanthu 2004, Kaitak 2005, Chanchu 2006); các hạn dự báo đánh giá gồm 12h, 24h, 36h và 48h.
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu lý thuyết và ứng dụng về bão
Luận án điểm lại có hệ thống các giai đoạn phát triển của lý thuyết chuyển động xoáy thuận nhiệt đới:
- Khái niệm dòng dẫn và mô hình chính áp: Kasahara (1957) sử dụng mô hình chính áp không phân kỳ, tách hàm dòng thành dòng quy mô lớn và nhiễu động xoáy địa phương, chứng minh tâm xoáy di chuyển gần vận tốc dòng dẫn nhưng lệch trái 90 độ so với gradient xoáy tuyệt đối trên mặt beta. Sanders và cộng sự (1975), George và Gray (1976), Chan và Gray (1982) nghiên cứu thực nghiệm xác định mực dòng dẫn tối ưu (700 hPa cho tốc độ, 500 hPa cho hướng, và trung bình lớp sâu 500–700 hPa hoặc 850–200 hPa). Khái niệm dòng trung bình lớp sâu (Deep Layer Mean — DLM) là cơ sở cho các mô hình chính áp như VICBAR (DeMaria và cộng sự, 1992) và WBAR (Weber, 2001).
- Độ nhạy của chuyển động bão với kích thước và cường độ: DeMaria (1985) và Fiorino và Elsberry (1989) chỉ ra quỹ đạo bão nhạy với vị trí ban đầu và kích thước hoàn lưu ngoài (bán kính 300–800 km và cường độ gió ngoài 300–1000 km) hơn là cường độ gió cực đại tại tâm. Smith và Ulrich (1990), Smith và Weber (1993) phát triển lý thuyết giải tích về thành phần phi đối xứng xoáy (xoắn beta — beta-gyres) điều khiển hướng di chuyển của bão.
- Mô hình hóa ba chiều và quá trình ẩm: Ooyama (1969) xây dựng mô hình 3 lớp đối xứng trục; Frank và Ritchie (1999) dùng mô hình MM5 chỉ ra vai trò của độ đứt gió đứng và ma sát lớp biên đối với tổ chức đối lưu; Kimball và Evans (2002) dùng mô hình phi thủy tĩnh khảo sát biến đổi bán kính gió cực đại ($R_{max}$); Zhu và cộng sự (2001) cùng Nguyễn Chi Mai và cộng sự (2002) phát triển các mô hình ba chiều tối thiểu 3 mực khảo sát dòng siêu gradient lớp biên.
Các phương pháp ban đầu hóa xoáy bão
Văn bản phân tích ba nhóm kỹ thuật chính trong bài toán ban đầu hóa xoáy:
- Xây dựng xoáy nhân tạo:
- Phương pháp động lực: Kurihara và cộng sự (1993, mô hình MMM của GFDL) tích phân một phiên bản đối xứng trục từ trạng thái tĩnh để tạo xoáy cân bằng, điều chỉnh dần gió tiếp tuyến về trường mục tiêu theo thời gian.
- Phương pháp kinh nghiệm: Iwasaki và cộng sự (1987, mô hình bão của JMA) sử dụng công thức khí áp bề mặt Fujita (1952), cấu trúc lõi ấm và hệ phương trình chuyển động trong tọa độ trụ để giải lặp trường gió; Davidson và Puri (1992), Davidson và cộng sự (1993, LAPS), Davidson và Weber (2000, TC-LAPS); Low-Nam (2001, NCAR-AFWA cho MM5) sử dụng profile Rankine sửa đổi; Smith (2005) đề xuất phương pháp giải tích thiết lập xoáy ba chiều cân bằng từ phân bố gió tiếp tuyến định trước.
- Phương pháp phân tích tách xoáy:
- Kurihara và cộng sự (1993) áp dụng toán tử lọc làm trơn đa bước trên lưới kinh vĩ và hệ tọa độ cực để tách trường nền, trường nhiễu phi xoáy và xoáy phân tích.
- Weber và Smith (1995) sử dụng phương pháp lọc Barnes (1964) kết hợp biến đổi Fourier nhanh (FFT) tần thấp để tách trường môi trường quy mô lớn, quy mô nhỏ và phân tích Fourier phương vị tách thành phần đối xứng (sóng số 0) và phi đối xứng (sóng số 1).
- Low-Nam (2001) giải phương trình Poisson đối với hàm dòng và thế vận tốc để bóc tách thành phần xoáy không phân kỳ và phân kỳ.
- Phương pháp kết hợp xoáy và đồng hóa số liệu:
- Cài xoáy trực tiếp bằng hàm trọng số cosine (Mathur, 1991) hoặc hàm suy giảm tuyến tính theo bán kính (Davidson và cộng sự, 1993).
- Đồng hóa dữ liệu giả biến phân bốn chiều 4DVar (Zou và Xiao, 2000) hoặc ba chiều 3DVar (Xiao và cộng sự, 2006) trên mô hình MM5.
Nghiên cứu trong nước và vị trí của luận án
Ở Việt Nam, Trịnh Văn Thư và Krishnamurti (1992) đã thử nghiệm ban đầu hóa xoáy Rankine cho mô hình nước nông một mực dự báo bão Betty (1987) và Dan (1989). Sau năm 2000, mô hình chính áp WBAR được khảo sát và cải tiến qua các công trình của Phan Văn Tân và cộng sự (2002), Bùi Hoàng Hải và Phan Văn Tân (2002), Võ Văn Hòa (2005), và Nguyễn Thị Minh Phương (2003, 2005, 2007 với phiên bản BARO_PH).
Bên cạnh mô hình chính áp, mô hình khu vực phân giải cao HRM được chuyển giao và ứng dụng nghiệp vụ qua đề tài của Kiều Thị Xin và cộng sự (2002). Lê Công Thành (2004) khẳng định HRM có kỹ năng dự báo bão ổn định hơn các mô hình chính áp khi bão đổi hướng. Các thử nghiệm trên mô hình ETA (Trần Tân Tiến và cộng sự, 2004) và MM5 (Hoàng Đức Cường, 2004; Phan Văn Tân và Bùi Hoàng Hải, 2003; Đặng Thị Hồng Nga và cộng sự, 2006) cũng khẳng định sự cần thiết của việc cài xoáy.
Bảng 1: Sai số vị trí tâm bão mùa bão 2005 của một số mô hình nghiệp vụ (trích từ văn bản)
| Mô hình |
Hạn 12h (km) |
Hạn 24h (km) |
Hạn 36h (km) |
Hạn 48h (km) |
| JTWC |
82,6 |
134,8 |
196,4 |
260,2 |
| JMA |
86,2 |
145,0 |
215,8 |
291,0 |
| CLIPER |
95,8 |
214,1 |
358,0 |
485,3 |
| BARO_PH |
83,5 |
138,2 |
210,5 |
295,4 |
Vị trí của luận án: Luận án không đi theo hướng tinh chỉnh mô hình chính áp đơn giản (WBAR/BARO_PH) vốn bị giới hạn vật lý, cũng không áp dụng các sơ đồ kinh nghiệm thô sơ (như Rankine gián đoạn). Thay vào đó, tác giả lựa chọn tiếp thu lý thuyết xoáy cân bằng giải tích của Smith (2005) và kỹ thuật lọc Fourier của Weber và Smith (1995), xây dựng một sơ đồ ban đầu hóa xoáy ba chiều hoàn chỉnh và tích hợp trực tiếp vào mô hình khu vực tà áp HRM.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung phân tích tách xoáy ba chiều
Luận án áp dụng phương pháp phân tích xoáy của Smith và Weber (1995). Tại mỗi điểm nút lưới $(i, j)$ trên từng mực mô hình, một trường vô hướng bất kỳ $f_{ij}$ được phân rã thành:
$$f_{ij} = f_{ij}^E + f_{ij}^V = (f_{ij}^{EL} + f_{ij}^{ES}) + (f_{ij}^{VS} + f_{ij}^{VA})$$
Trong đó:
- $f_{ij}^E$ là trường môi trường; $f_{ij}^V$ là trường xoáy.
- $f_{ij}^{EL}$ là trường môi trường quy mô lớn; $f_{ij}^{ES}$ là trường môi trường quy mô nhỏ.
- $f_{ij}^{VS}$ là thành phần xoáy đối xứng (sóng số 0); $f_{ij}^{VA}$ là thành phần xoáy phi đối xứng (sóng số 1).
Quy trình bóc tách trường môi trường quy mô lớn $f_{ij}^{EL}$ được thực hiện qua các bước:
- Làm trơn sơ bộ bằng trung bình 4 điểm xung quanh trên một miền có kích thước $10^\circ \times 10^\circ$ kinh vĩ có tâm xấp xỉ trùng tâm bão quan trắc:
$$\bar{f}{i,j} = \frac{1}{4}(f{i,j-1} + f_{i,j+1} + f_{i-1,j} + f_{i+1,j})$$
Trường nhiễu động sơ bộ được xác định: $f'{i,j} = f{i,j} - \bar{f}_{i,j}$.
- Sử dụng sơ đồ lặp Barnes (1964) kết hợp bộ lọc tần thấp thông qua phép biến đổi Fourier nhanh một chiều (1D-FFT) để triệt tiêu các sóng ngắn hơn ngưỡng xác định, thu được trường xấp xỉ bước đầu tiên $f_{ij}^{EL(1)}$.
- Thiết lập trường phần dư $f_{ij}^{D(1)} = \bar{f}{i,j} - f{ij}^{EL(1)}$ và lặp làm trơn cho đến khi độ lệch giữa hai bước lặp liên tiếp đạt ngưỡng hội tụ.
- Xác định tâm xoáy phân tích (dựa vào cực đại độ xoáy, nhiệt độ hoặc cực tiểu khí áp), sau đó phân tích chuỗi Fourier theo góc phương vị xung quanh tâm để tách sóng số 0 ($f_{ij}^{VS}$) và sóng số 1 ($f_{ij}^{VA}$). Phần còn lại sau khi loại bỏ xoáy chính là trường môi trường quy mô nhỏ $f_{ij}^{ES}$. Quy trình áp dụng cho hai thành phần gió ngang $(u, v)$ và nhiệt độ ($T$).
Lý thuyết xây dựng xoáy ba chiều cân bằng
Để tạo xoáy nhân tạo, luận án ứng dụng lý thuyết xoáy cân bằng ba chiều của Smith (2005). Trong hệ tọa độ trụ đối xứng trục $(r, z)$, các trường khí áp $p(r, z)$, mật độ $\rho(r, z)$ và nhiệt độ $T(r, z)$ được thiết lập cân bằng động lực học với trường gió tiếp tuyến $V(r, z)$ thông qua hệ phương trình:
- Phương trình cân bằng gió gradient:
$$\frac{1}{\rho}\frac{\partial p}{\partial r} = \frac{V^2}{r} + fV$$
- Phương trình cân bằng thủy tĩnh:
$$\frac{1}{\rho}\frac{\partial p}{\partial z} = -g$$
- Phương trình trạng thái khí lý tưởng: $p = \rho R T$.
Phân bố gió tiếp tuyến đối xứng được mô tả bằng hàm giải tích kinh nghiệm của DeMaria (1987):
$$V(r) = V_m \left(\frac{r}{r_m}\right) \exp\left{\frac{1}{b}\left[1 - \left(\frac{r}{r_m}\right)^b\right]\right}$$
Trong đó $V_m$ là vận tốc gió tiếp tuyến cực đại, $r_m$ là bán kính gió cực đại ($R_{max}$), và $b$ là tham số hình dạng điều khiển mức độ suy giảm gió ra phía ngoài, được tính toán thông qua bán kính gió 15 m/s ($R_{15}$). Cấu trúc đứng của xoáy được kiểm soát bằng hàm trọng số thẳng đứng $W(z)$ hoặc $W(p)$.
Công cụ mô hình số trị và dữ liệu
- Mô hình WRF (Weather Research and Forecasting): Sử dụng phiên bản ARW với lưới tọa độ thẳng đứng thủy tĩnh $\eta$, thiết lập các kịch bản lý tưởng đối xứng để kiểm chứng sự cân bằng động lực học và quá trình tiến triển của xoáy ba chiều trong 24 giờ tích phân.
- Mô hình HRM (High-Resolution Regional Model): Mô hình khu vực hạn chế thủy tĩnh, sử dụng hệ phương trình nguyên thủy khí quyển, lưới ngang xen kẽ Arakawa C, tọa độ thẳng đứng lai (hybrid vertical coordinate) gồm 20 hoặc 30 mực. Luận án can thiệp mã nguồn để xây dựng phiên bản HRM_TC tích hợp mô-đun ban đầu hóa xoáy ba chiều.
- Dữ liệu đầu vào: Trường phân tích khí quyển toàn cầu GME (độ phân giải lưới tam giác xấp xỉ 40 km, 31 mực thẳng đứng) cung cấp điều kiện ban đầu và điều kiện biên. Thông tin chỉ thị bão (tâm bão quan trắc, $P_c$, $V_{max}$, $R_{15}$, $R_{max}$) trích xuất từ các bản tin của các trung tâm dự báo bão khu vực và quốc tế.
- Chỉ tiêu đánh giá định lượng: Sai số vị trí tâm bão (Position Error — PE), sai số dọc quỹ đạo (Along-Track Error — ATE), sai số ngang quỹ đạo (Cross-Track Error — CTE), độ lệch chuẩn sai số (Standard Deviation of Absolute Error — SDA), và chỉ số kỹ năng dự báo tương đối (Skill) so với phương án đối chứng:
$$\text{Skill} = \frac{\text{PE}{\text{control}} - \text{PE}{\text{bogus}}}{\text{PE}_{\text{control}}} \times 100%$$
Nội dung chính theo từng chương
Chương 1: Tổng quan các nghiên cứu về bão và ban đầu hóa xoáy bão
Chương 1 hệ thống hóa các thành tựu lý thuyết về cơ chế chuyển động của bão (dòng dẫn, hiệu ứng tương tác beta, tương tác hoàn lưu với độ đứt gió thẳng đứng). Tác giả so sánh chi tiết các dạng profile gió tiếp tuyến lý thuyết: Rankine cổ điển, Rankine sửa đổi, profile Fujita và profile DeMaria (1987). Văn bản chỉ ra rằng profile Rankine tuy đơn giản nhưng có điểm gãy khúc không khả vi tại $r_m$, trong khi profile DeMaria khắc phục được nhược điểm này và dễ dàng kiểm soát cấu trúc hoàn lưu ngoài thông qua tham số $b$.
Chương này cũng phân tích sâu các kỹ thuật cài xoáy và đồng hóa số liệu (3DVar, 4DVar). Đối với điều kiện tính toán tại Việt Nam thời điểm năm 2008, phương pháp biến phân 4 chiều đòi hỏi tài nguyên máy tính rất lớn và số liệu quan trắc biển nghèo nàn dễ dẫn đến sai số nền. Do đó, việc xây dựng một sơ đồ ban đầu hóa xoáy bằng phương pháp giải tích cân bằng kết hợp lọc trường Fourier là giải pháp thực tế và khả thi nhất. Tác giả kết luận việc lựa chọn mô hình HRM để nâng cấp thành HRM_TC là hướng đi đúng đắn nhằm khắc phục các nhược điểm của mô hình chính áp WBAR và nâng cao độ chính xác dự báo bão tại Biển Đông.
Chương 2: Nghiên cứu phát triển sơ đồ ban đầu hóa xoáy ba chiều
Chương 2 mô tả chi tiết thuật toán bóc tách trường xoáy và thuật toán xây dựng trường khối lượng cân bằng từ trường gió tiếp tuyến. Thuật toán phân tích xoáy cho phép loại bỏ xoáy sai lệch trong trường phân tích toàn cầu GME mà không làm biến dạng trường môi trường quy mô lớn xung quanh.
Để kiểm chứng tính vững và khả thi của sơ đồ xoáy cân bằng Smith (2005), tác giả thiết kế các thử nghiệm lý tưởng 3 chiều trên mô hình WRF:
- Kịch bản TH1: Xoáy khô cân bằng trong môi trường không ma sát.
- Kịch bản TH2: Xoáy khô cân bằng có tính đến ma sát lớp biên.
- Kịch bản TH3: Xoáy cân bằng được bổ sung ban đầu hóa ẩm (độ ẩm tương đối được gán giá trị cao ở vùng tâm bão).
Kết quả mô phỏng tích phân trong 24 giờ cho thấy:
- Trường gió tiếp tuyến và trường áp suất mực biển duy trì cấu trúc đối xứng ổn định; biến đổi khí áp tại tâm diễn ra êm thuận, không xuất hiện các dao động sóng trọng trường giả mạo.
- Cấu trúc lõi ấm (warm core) tại tầng đối lưu trên và dị thường áp suất thấp ở bề mặt được tái lập cân bằng hoàn hảo.
- Ở kịch bản ẩm (TH3), sau 24 giờ tích phân, mô hình tự hình thành hoàn lưu thứ cấp chân thực: dòng hội tụ mạnh ở tầng biên thấp (mực 1), dòng thăng mạnh quanh bán kính gió cực đại kèm theo ngưng kết ẩm tạo hàm lượng nước trong mây, và dòng phân kỳ tỏa ra ở tầng cao (mực 11). Kết quả này khẳng định sơ đồ xoáy cân bằng hoàn toàn tương thích và ổn định khi áp dụng vào các mô hình số trị ba chiều nguyên thủy.
Chương 3: Áp dụng sơ đồ ban đầu hóa xoáy ba chiều dự báo quĩ đạo bão
Chương 3 trình bày việc tích hợp sơ đồ vào mô hình HRM, xây dựng quy trình công nghệ cho phiên bản HRM_TC. Tác giả thực hiện một chuỗi các thí nghiệm độ nhạy có hệ thống trên các ca bão thực tế (Imbudo 2003, Conson 2004, Kaitak 2005, Chanchu 2006) để xác lập bộ tham số khả dụng:
- Thí nghiệm bán kính gió cực đại ($R_{max}$): Khảo sát 3 phương án RM1 ($R_{max} = 60\text{ km}$), RM2 ($R_{max} = 90\text{ km}$), và RM3 ($R_{max} = 120\text{ km}$). Kết quả cho thấy phương án RM2 ($R_{max} = 90\text{ km}$) mang lại sai số vị trí trung bình (MPE) nhỏ nhất và độ lệch chuẩn sai số (SDA) thấp nhất trên các hạn dự báo.
- Thí nghiệm bán kính gió 15 m/s ($R_{15}$): Khảo sát 4 phương án kích thước hoàn lưu ngoài: S1 (200 km), S2 (250 km), S3 (300 km) và S4 (400 km). Kết quả chỉ ra rằng cấu trúc hoàn lưu gió ngoài có ảnh hưởng rất lớn đến tốc độ và hướng di chuyển của bão; phương án S2–S3 (250–300 km) phản ánh sát thực nhất kích thước bão khu vực Biển Đông, cải thiện rõ rệt kỹ năng dự báo so với phương án đối chứng (không cài xoáy).
- Thí nghiệm hàm trọng số thẳng đứng ($W$): So sánh 4 dạng hàm trọng số (W1 đến W4). Hàm trọng số duy trì hoàn lưu xoáy mạnh từ mặt biển lên đến 700 hPa và suy giảm dần lên 200 hPa cho kết quả sai số vị trí thấp nhất và phân bố trường áp suất mực biển đối xứng mịn nhất.
- Thí nghiệm kết hợp trường gió tiếp tuyến: So sánh giữa phương án chỉ dùng gió xoáy nhân tạo thuần túy (M1) và phương án kết hợp gió nhân tạo ở vùng trung tâm với gió xoáy phân tích ở vùng ngoài (M2). Phương án M2 cho thấy sự chuyển tiếp trường động lực trơn tru hơn, triệt tiêu sự bất liên tục tại biên ngoài của xoáy và nâng cao kỹ năng dự báo tổng thể.
Chương 4: Thử nghiệm sơ đồ ban đầu hóa xoáy mới
Chương 4 kiểm chứng hiệu quả của mô hình HRM_TC với bộ tham số đã xác lập ở Chương 3 ($R_{max} = 90\text{ km}$, kết hợp gió M2, hàm trọng số chuẩn) thông qua các trường hợp bão dự báo thử nghiệm độc lập.
Hai ca bão điển hình được phân tích chi tiết:
- Bão Chanthu (00Z ngày 11/6/2004): Trường hợp bão di chuyển nhanh và sắp đổ bộ vào bờ biển Trung Bộ. Trường phân tích ban đầu GME thể hiện xoáy bão rất mờ nhạt và tâm bão lệch đáng kể so với vị trí thực tế. Sau khi áp dụng sơ đồ ban đầu hóa xoáy, trường áp suất mực biển tái hiện hoàn lưu bão sâu sắc và tâm bão trùng khớp với quan trắc. Phương án cài xoáy (bogus) dự báo chính xác khu vực đổ bộ với sai số quỹ đạo giảm mạnh so với phương án không cài xoáy (nobogus).
- Bão Koni (00Z ngày 20/7/2003): Trường hợp bão di chuyển phức tạp trên Biển Đông trong bối cảnh có sự xuất hiện đồng thời của cơn bão rất mạnh Imbudo ở rìa phía đông miền tính. HRM_TC đã bóc tách chính xác xoáy bão Koni và tái lập trường dòng dẫn quy mô lớn hợp lý, giúp mô phỏng chuyển động của bão sát với quỹ đạo quan trắc thực tế.
Bảng 2: Tổng hợp sai số dự báo trung bình và kỹ năng của sơ đồ ban đầu hóa xoáy (trích từ văn bản Chương 4)
| Đại lượng đánh giá |
Hạn 12h |
Hạn 24h |
Hạn 36h |
Hạn 48h |
| Sai số vị trí PE (km): |
|
|
|
|
| - Không ban đầu hóa (nobogus) |
124,5 |
198,2 |
275,6 |
362,4 |
| - Có ban đầu hóa xoáy (bogus) |
78,3 |
122,6 |
176,4 |
231,5 |
| Kỹ năng dự báo tương đối (%): |
|
|
|
|
| - Kỹ năng theo PE |
+37,1% |
+38,1% |
+36,0% |
+36,1% |
| - Kỹ năng theo sai số dọc (ATE) |
Dương |
Dương |
Dương |
Dương |
| - Kỹ năng theo sai số ngang (CTE) |
Dương |
Dương |
Dương |
Dương |
Kết quả tổng hợp cho thấy sơ đồ ban đầu hóa xoáy mới giúp giảm sai số vị trí trung bình từ 36% đến 38% ở tất cả các hạn dự báo từ 12h đến 48h. Cả hai thành phần sai số dọc quỹ đạo (ATE — phản ánh sai số tốc độ di chuyển) và sai số vuông góc quỹ đạo (CTE — phản ánh sai số lệch hướng) đều được cải thiện vượt trội so với phương án đối chứng.
Kết quả và những đóng góp mới
Đóng góp mới về mặt khoa học và phương pháp
- Luận án đã nghiên cứu, phát triển và xây dựng thành công một sơ đồ ban đầu hóa xoáy ba chiều mới dựa trên lý thuyết xoáy cân bằng giải tích (Smith, 2005) kết hợp với kỹ thuật phân tích lọc trường Fourier và lặp Barnes (Weber và Smith, 1995). Sơ đồ có tính tự nhất quán cao về mặt thủy nhiệt động lực học, có thể ứng dụng linh hoạt cho cả các mô phỏng lý tưởng lẫn các mô hình dự báo thời tiết số trị với số liệu quan trắc thực tế.
- Đã áp dụng thành công sơ đồ ban đầu hóa xoáy ba chiều mới vào mô hình khu vực phân giải cao HRM, phát triển mô hình này thành một phiên bản chuyên dụng mới mang tên HRM_TC. Phiên bản này tích hợp hoàn chỉnh quy trình xử lý dữ liệu chỉ thị bão, tách xoáy phân tích, kiến tạo xoáy cân bằng và kết hợp trường động lực, giúp mô hình vừa duy trì năng lực dự báo thời tiết chung vừa nâng cao năng lực dự báo chuyển động bão.
- Luận án đã khảo sát hệ thống và xác lập được bộ tham số tối ưu cho sơ đồ ban đầu hóa xoáy của mô hình HRM_TC (bao gồm bán kính gió cực đại $R_{max} = 90\text{ km}$, phân bố kích thước hoàn lưu ngoài $R_{15}$, hàm trọng số thẳng đứng phù hợp và phương thức kết hợp gió tiếp tuyến M2), phản ánh chân thực đặc thù cấu trúc xoáy thuận nhiệt đới hoạt động trên khu vực Biển Đông.
Đóng góp thực tiễn và giải pháp đề xuất
- Cung cấp một công cụ mô hình số trị tà áp ba chiều (HRM_TC) có kỹ năng dự báo quỹ đạo bão ổn định tới hạn 48 giờ, góp phần chuyển dịch quy trình dự báo bão tại Việt Nam từ chỗ phụ thuộc chủ yếu vào synop và mô hình thống kê sang ứng dụng thực chất sản phẩm mô hình số trị phân giải cao.
- Đề xuất đưa sơ đồ HRM_TC vào thử nghiệm nghiệp vụ tại Trung tâm Dự báo Khí tượng Thủy văn Trung ương, kết hợp việc khai thác kịp thời các thông tin chỉ thị bão (từ ảnh mây vệ tinh và quan trắc thám sát) để phục vụ công tác cảnh báo bão sớm phòng chống thiên tai.
Hạn chế và hướng nghiên cứu tiếp
Dựa trên các nội dung kỹ thuật mà tác giả trình bày trong luận án:
- Chưa ban đầu hóa trường độ ẩm: Sơ đồ ban đầu hóa xoáy trong mô hình HRM_TC hiện tại chỉ tập trung vào việc cân bằng trường gió ngang $(u, v)$ và trường nhiệt độ $(T)$. Trường độ ẩm chưa được ban đầu hóa tương thích (tác giả giả định độ ẩm ít ảnh hưởng trực tiếp đến quỹ đạo bão trong ngắn hạn), điều này có thể hạn chế khả năng mô phỏng quá trình biến đổi cường độ bão và giáng thủy chi tiết.
- Trường độ cao địa thế vị phụ thuộc vào cân bằng thủy tĩnh: Sơ đồ không can thiệp trực tiếp vào trường độ cao địa thế vị mà để mô hình tự hiệu chỉnh thông qua quan hệ tích phân thủy tĩnh trong quá trình giải số.
- Hiện tượng điều chỉnh ban đầu (spin-up): Dù xoáy nhân tạo được thiết kế cân bằng giải tích nhưng khi đưa vào trường nền tà áp thực tế của mô hình HRM, hệ thống vẫn cần một khoảng thời gian tích phân đầu tiên để các trường tương thích hoàn toàn.
- Dữ liệu thử nghiệm: Do hạn chế nguồn số liệu phân tích lưu trữ, các thử nghiệm dự báo số trị mới được thực hiện trên một số lượng trường hợp bão chọn lọc trong các mùa bão 2003–2006, cần tiếp tục mở rộng quy mô mẫu trên nhiều năm để hoàn thiện bộ tham số thống kê.
Hướng nghiên cứu tiếp theo được văn bản định hướng gồm:
- Nghiên cứu thuật toán ban đầu hóa trường ẩm và tham số hóa vi vật lý mây đi kèm xoáy nhân tạo.
- Tích hợp số liệu viễn thám (gió vệ tinh, radar thời tiết) vào sơ đồ ban đầu hóa thông qua các kỹ thuật đồng hóa dữ liệu tiên tiến hơn (như đồng hóa biến phân 3DVar/4DVar) khi điều kiện tài nguyên tính toán đáp ứng được.
Giá trị tham khảo
Luận án là tài liệu học thuật có giá trị ứng dụng và tham khảo rõ rệt trong chuyên ngành Khí tượng học và Dự báo thời tiết số trị:
- Đối với nghiên cứu sinh, học viên cao học và sinh viên ngành Khí tượng Thủy văn: Luận án cung cấp cơ sở lý thuyết chuẩn mực về động lực học xoáy thuận nhiệt đới, phương pháp giải tích xoáy cân bằng ba chiều của Smith (2005), và quy trình toán học chi tiết của phương pháp lọc Barnes kết hợp biến đổi Fourier nhanh (1D-FFT) để phân tách quy mô trường khí quyển.
- Đối với các nhà nghiên cứu phát triển mô hình: Cung cấp giải pháp mã nguồn và phương pháp luận can thiệp điều kiện ban đầu cho mô hình số trị tà áp (cụ thể là mô hình HRM), phương pháp thiết kế các kịch bản thử nghiệm độ nhạy tham số động lực học ($R_{max}$, $R_{15}$, hàm trọng số thẳng đứng).
- Đối với cán bộ dự báo khí tượng nghiệp vụ: Các kết quả thực nghiệm ở Chương 3 và Chương 4 cung cấp bộ tham số thực tế về cấu trúc bão trên Biển Đông, làm tài liệu tham chiếu tin cậy khi cài đặt xoáy nhân tạo hoặc hiệu chỉnh sai số dự báo quỹ đạo bão tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia.
Câu hỏi thường gặp
1. Tại sao cần phải thực hiện ban đầu hóa xoáy bão trong mô hình số trị dự báo thời tiết?
Do bão nhiệt đới hình thành và di chuyển trên các vùng biển rộng lớn nơi mạng lưới quan trắc synop truyền thống rất thưa thớt, các trường phân tích toàn cầu (như GME) thường biểu diễn xoáy bão rất yếu, mờ nhạt hoặc sai lệch vị trí tâm. Nếu giữ nguyên trường ban đầu này để tích phân, sai số vị trí sẽ tích lũy và tăng nhanh theo thời gian dự báo. Kỹ thuật ban đầu hóa xoáy giúp bóc tách xoáy phân tích sai lệch ra khỏi trường môi trường và thay thế bằng một xoáy nhân tạo có cấu trúc, cường độ và vị trí sát thực tế quan trắc, từ đó nâng cao độ chính xác dự báo quỹ đạo.
2. Thuật toán phân tích xoáy của luận án phân chia một trường vô hướng thành những thành phần nào?
Dựa trên phương pháp của Smith và Weber (1995), sơ đồ phân chia một trường vô hướng bất kỳ $f_{ij}$ thành 4 thành phần riêng biệt:
- Trường môi trường quy mô lớn ($f_{ij}^{EL}$): tách bằng phương pháp làm trơn sơ bộ kết hợp thuật toán lặp Barnes và bộ lọc tần thấp 1D-FFT.
- Trường môi trường quy mô nhỏ ($f_{ij}^{ES}$): phần còn lại của trường môi trường sau khi loại trừ xoáy.
- Trường xoáy đối xứng ($f_{ij}^{VS}$): xác định là sóng số 0 trong phân tích Fourier theo góc phương vị quanh tâm bão.
- Trường xoáy phi đối xứng ($f_{ij}^{VA}$): xác định là sóng số 1 trong phân tích Fourier theo góc phương vị (phản ánh thành phần xoắn beta).
3. Lý thuyết xoáy cân bằng của Smith (2005) được ứng dụng để tạo xoáy nhân tạo như thế nào?
Từ một phân bố gió tiếp tuyến đối xứng được chỉ định theo bán kính bằng hàm giải tích DeMaria (1987) và hàm trọng số thẳng đứng, tác giả tích phân hệ phương trình cân bằng gradient và phương trình cân bằng thủy tĩnh trong hệ tọa độ trụ $(r, z)$. Thuật toán này cho phép suy dẫn trực tiếp và đồng thời các trường khối lượng (mật độ $\rho$, áp suất $p$, và nhiệt độ $T$ thể hiện cấu trúc lõi ấm) ở trạng thái cân bằng động lực hoàn chỉnh mà không cần phải tích phân tốn thời gian một mô hình đối xứng trục phụ trợ.
4. Vì sao tác giả phát triển sơ đồ trên mô hình HRM thay vì mô hình WBAR?
Mô hình WBAR là mô hình chính áp một mực (sử dụng gió trung bình lớp sâu), có giả thiết đơn giản hóa cao về khí quyển nên chỉ dự báo tốt các cơn bão mạnh có quỹ đạo thẳng; khi bão gặp các hệ thống thời tiết tà áp phức tạp hoặc đổi hướng đột ngột, WBAR cho sai số rất lớn. Mô hình HRM là mô hình khu vực tà áp ba chiều phân giải cao đầy đủ, có khả năng nắm bắt các tương tác nhiệt động lực phức tạp. Việc tích hợp sơ đồ ban đầu hóa xoáy vào HRM cho phép tận dụng ưu thế vật lý của mô hình tà áp để dự báo chính xác quỹ đạo bão ở các hạn dài hơn tới 48 giờ.
5. Bộ tham số tối ưu nào đã được xác lập cho mô hình HRM_TC trên khu vực Biển Đông?
Qua các thí nghiệm độ nhạy ở Chương 3, luận án đã xác lập bộ tham số khả dụng tối ưu gồm:
- Bán kính gió cực đại: $R_{max} = 90\text{ km}$ (cho sai số vị trí trung bình nhỏ nhất so với các mức 60 km hay 120 km).
- Bán kính gió 15 m/s: $R_{15}$ trong khoảng 250–300 km.
- Hàm trọng số thẳng đứng: dạng suy giảm từ tầng đối lưu dưới lên tầng đối lưu trên phù hợp với số liệu thám sát bão thực tế.
- Phương án kết hợp gió: phương án M2 (kết hợp gió xoáy nhân tạo ở vùng trung tâm với gió xoáy phân tích ở vùng ngoài hoàn lưu).
Kết luận
Luận án tiến sĩ của Bùi Hoàng Hải đã giải quyết trọn vẹn bài toán xây dựng và ứng dụng sơ đồ ban đầu hóa xoáy ba chiều cho mô hình dự báo số trị khu vực HRM tại Việt Nam. Bằng việc kết hợp chặt chẽ giữa thuật toán phân tích lọc trường Fourier và lý thuyết xoáy ba chiều cân bằng giải tích, công trình đã phát triển thành công phiên bản mô hình chuyên dụng HRM_TC và xác lập được bộ tham số cấu trúc xoáy tối ưu cho bão trên Biển Đông. Kết quả thử nghiệm dự báo số trị độc lập khẳng định mô hình HRM_TC giúp giảm sai số vị trí tâm bão từ 36% đến 38% ở các hạn dự báo từ 12 đến 48 giờ so với phương án không ban đầu hóa xoáy, mang lại đóng góp thiết thực cho khoa học khí tượng và công tác nghiệp vụ dự báo bão tại Việt Nam.
