Tổng quan luận án

Bão và áp thấp nhiệt đới (gọi chung là xoáy thuận nhiệt đới - XTNĐ) là hiện tượng thời tiết nguy hiểm gây mưa to, gió lớn, lũ lụt diện rộng và nước dâng ven biển. Theo Tổ chức Khí tượng Thế giới (WMO), Việt Nam nằm ở ổ bão khu vực Tây Bắc Thái Bình Dương với trung bình 30 cơn bão mỗi năm, chiếm khoảng 38% tổng số bão toàn cầu. Hàng năm, nước ta chịu nhiều thiệt hại về người và vật chất do bão gây ra. Việc ứng dụng các mô hình số trị để nâng cao độ chính xác và kéo dài thời hạn dự báo quỹ đạo bão là yêu cầu cấp thiết nhằm phục vụ công tác phòng chống thiên tai.

Tại Việt Nam, trước thời điểm nghiên cứu này, chưa có công trình nào đặt vấn đề dự báo quỹ đạo bão hạn đến 5 ngày (120 giờ), trong khi các trung tâm khí tượng quốc tế lớn như Nhật Bản, Hoa Kỳ hay Trung Quốc đã đưa hạn dự báo này vào nghiệp vụ. Nguyên nhân chính khiến chất lượng dự báo bão của các mô hình số khu vực chưa cao là do sai số vị trí và cường độ bão trong các trường phân tích và dự báo toàn cầu còn lớn, cùng với các xấp xỉ vật lý và thuật toán tồn tại trong từng mô hình. Để cải thiện trường ban đầu, kỹ thuật cài xoáy giả ba chiều trước đây tuy có hiệu quả về mặt định vị tâm bão ban đầu nhưng chưa mô tả được cấu trúc bất đối xứng của bão thực tế. Do đó, khoảng trống nghiên cứu được luận án xác định là: cần kết hợp kỹ thuật nuôi nhiễu cho khu vực có bão (nhiễu xoáy) nhằm mô tả cấu trúc bão với kỹ thuật nuôi nhiễu khu vực bên ngoài bão (nhiễu môi trường), đồng thời kết hợp thay đổi các sơ đồ đối lưu trong mô hình khu vực để xây dựng một hệ thống dự báo tổ hợp quỹ đạo bão hạn 5 ngày phù hợp với điều kiện Việt Nam.

Luận án xác định 2 mục đích nghiên cứu chính (được đánh số theo văn bản):

  1. Cải tiến và áp dụng thành công kỹ thuật nuôi nhiễu phát triển trên mô hình RAMS khi có bão trên khu vực Tây Bắc Thái Bình Dương.
  2. Xây dựng được quy trình dự báo tổ hợp quỹ đạo bão ảnh hưởng đến Việt Nam hạn 5 ngày bằng phương pháp tổ hợp dựa trên phương pháp nuôi nhiễu phát triển.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Quỹ đạo di chuyển của các cơn bão (xoáy thuận nhiệt đới) ảnh hưởng đến Việt Nam.
  • Phạm vi không gian: Khu vực Tây Bắc Thái Bình Dương và Biển Đông, tập trung vào các cơn bão đổ bộ hoặc tác động trực tiếp đến đất liền và vùng biển Việt Nam.
  • Phạm vi thời gian và tư liệu: Chuỗi số liệu thống kê khí hậu bão thời kỳ 1959–2008 và 1961–2010; tập số liệu thử nghiệm gồm 199 trường hợp dự báo bão trong giai đoạn 5 năm (2009–2013) ở các hạn từ 24 đến 120 giờ, với các trường hợp điển hình như bão Megi (2010), Nesat (2011), Nalgae (2011), Washi (2011), Gaemi (2012) và Bopha (2012).

Tổng quan tài liệu và vị trí của luận án

Luận án đã hệ thống hóa các cơ sở lý thuyết và công trình nghiên cứu về dự báo tổ hợp (EPS) từ quốc tế đến trong nước:

  • Lý thuyết về độ bất định khí quyển: Khởi nguồn từ Lorenz (1963, 1969) với hiệu ứng cánh bướm, chỉ ra rằng các sai số nhỏ ở trường ban đầu sẽ khuếch đại theo thời gian do tính phi tuyến của khí quyển. Khả năng dự báo bị giới hạn bởi sai số phân tích ban đầu và sai số mô hình. Tiếp đó, Epstein (1969) đề xuất phương pháp ngẫu nhiên động lực; Leith đề xuất phương pháp ước lượng tối ưu.
  • Các phương pháp tạo nhiễu trường ban đầu: Hoffman và Kalnay (1983) phát triển phương pháp dự báo trung bình trễ (LAF); Ebisuzaki và Kalnay (1991) cải tiến thành phương pháp dự báo trung bình trễ có chuẩn hóa (SLAF). Toth và Kalnay (1993, 1997) xây dựng kỹ thuật nuôi những dao động phát triển (Breeding of Growing Modes - BGM/BV) qua 2 phiên bản (sử dụng 1 nhiễu đơn lẻ và sử dụng cặp nhiễu dương - âm). Kalnay (2003) và Du (2007) phân tích các đặc tính của dao động phát triển (tính sát thực, tính phân kỳ và tính trực giao).
  • Các biến thể cải tiến kỹ thuật nuôi nhiễu và véc tơ kỳ dị: Wang và Bishop (2003) phát triển kỹ thuật nuôi nhiễu ẩn (Masked Breeding); Lee (2006) tại Cục Khí tượng Hàn Quốc (KMA) áp dụng kỹ thuật nuôi nhiễu kết hợp quay nhân tố (RBV); Wei và cộng sự (2007) áp dụng kỹ thuật chuyển vị tổ hợp (ET-Breeding) tại NCEP. Buizza (1994), Palmer và cộng sự (1998) nghiên cứu kỹ thuật véc tơ kỳ dị (SV) dựa trên mô hình tiếp tuyến tuyến tính (TLM) và phương pháp liên hợp tại ECMWF; Barkmeijer và cộng sự (1998), Magnusson và cộng sự (2008a, 2008b) nghiên cứu công thức tổng năng lượng và kỹ thuật BV-EOF; Wang và Kann (2005) ứng dụng trong hệ thống ALADIN-LAEF.
  • Phương pháp đồng hóa và lọc Kalman: Houtekamer (1995) với phương pháp nhiễu thám sát (PO); Evensen (1994) đề xuất lọc Kalman tổ hợp (EnKF); Anderson (2001), Bishop (2001), Whitaker và Hamill (2002) phát triển các bộ lọc căn bậc hai; Wang và Bishop (2003) đưa ra kỹ thuật ETKF. Các nghiên cứu so sánh của Hamill và cộng sự (2000b), Bowler (2006), Descamps và Talagrand (2007), Magnusson và cộng sự (2008b) làm rõ ưu nhược điểm của từng phương pháp.
  • Nghiên cứu về điều kiện biên ngang (LBC) và đa cơ chế vật lý: Warner và cộng sự (1997), Du và Tracton (1999), Nutter và Xue (2004a, 2004b) làm rõ tác động của LBC đến độ tán tổ hợp. Houtekamer và cộng sự (1995), Stensrud và cộng sự (2000), Chen và cộng sự (2003) nghiên cứu sai số mô hình từ việc thay đổi tham số vật lý. Điển hình là hệ thống SREF của NCEP (Du, 2011) kết hợp 3 mô hình (ETA, RSM, WRF) cùng các sơ đồ đối lưu khác nhau.
  • Các hệ thống dự báo bão tổ hợp khu vực: Queenie và cộng sự (2005) tại Đài Khí tượng Hồng Kông (HKO); Choo (2006) với hệ thống TMEPS tại KMA; Chen và cộng sự (2009) tại Trung tâm Khí tượng Quốc gia Trung Quốc (NMC/Bắc Kinh) sử dụng T213L31 kết hợp WRF; Huang và cộng sự (2007) với hệ thống GRAPES_TCM tại Thượng Hải; Yamaguchi (2012) tại Cục Khí tượng Nhật Bản (JMA).
  • Tình hình nghiên cứu trong nước: Đinh Văn Ưu (2009), Nguyễn Văn Thắng (2010) cung cấp số liệu khí hậu bão; các nghiên cứu trước đây tại Việt Nam chủ yếu áp dụng cài xoáy giả đối xứng trên mô hình đơn lẻ với hạn dự báo ngắn (24–72 giờ).

Khoảng trống luận án lựa chọn giải quyết: Các phương pháp nuôi nhiễu thông thường trên thế giới chủ yếu áp dụng cho quy mô vừa hoặc quy mô lớn mà chưa xử lý chuyên biệt cấu trúc xoáy khi có bão trong mô hình khu vực hạn chế. Luận án giải quyết vấn đề này bằng cách kết hợp đồng thời phương án nuôi nhiễu môi trường và nhiễu xoáy (có cài xoáy giả 3 chiều), thực hiện chuẩn hóa và trực giao hóa, đồng thời kết hợp 3 sơ đồ đối lưu trong mô hình RAMS để thiết lập hệ thống dự báo tổ hợp quỹ đạo bão hạn 5 ngày tại Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án dựa trên nền tảng lý thuyết động lực học khí quyển phi tuyến và lý thuyết dự báo tổ hợp (Lorenz, 1963; Kalnay, 2003). Nguyên lý của phương pháp nuôi những dao động phát triển (BV) là tận dụng quá trình tích phân mô hình phi tuyến để lọc ra các dạng nhiễu có tốc độ tăng trưởng nhanh nhất tương ứng với sai số phân tích ban đầu ("sai số ngày").

Sai số dự báo sau một chu kỳ nuôi từ thời điểm $t_i$ đến $t_j$ được xác định qua hiệu số giữa các tích phân mô hình: $$\delta x_i^f(t_j) = x_i^{fp}(t_j) - x_i^{fn}(t_j)$$ Trong đó $x_i^{fp}$ và $x_i^{fn}$ là các trường dự báo tích hợp từ cặp nhiễu dương và âm. Nhiễu phân tích được xác định lại qua việc co tỷ lệ (scaling down) bằng một hằng số chuẩn hóa $c$: $$\delta x^a(t_j) = c \cdot \delta x_i^f(t_j)$$

Phương pháp nghiên cứu thực nghiệm và công cụ cụ thể bao gồm:

  1. Mô hình số trị: Sử dụng mô hình khí quyển khu vực RAMS (Regional Atmospheric Modeling System). Cấu hình miền tính bao trùm khu vực Tây Bắc Thái Bình Dương và Biển Đông với hệ thống phương trình động lực đầy đủ, điều kiện biên ngang (LBC) cập nhật từ trường dự báo GFS.
  2. Thiết lập đa cơ chế vật lý: Luận án khảo sát và sử dụng 3 sơ đồ tham số hóa đối lưu trong mô hình RAMS:
    • Sơ đồ KUO
    • Sơ đồ Kain-Fritsch (KF)
    • Sơ đồ Kain-Fritsch cải tiến bởi Nguyễn Minh Trường (KFCT)
  3. Quy trình nuôi nhiễu kép:
    • Tạo 6 nhân ban đầu (ký hiệu từ D1 đến D6) từ phân bố chuẩn Gauss với độ lớn nhiễu tốc độ gió là 3 m/s và nhiệt độ là 1°C.
    • Cài xoáy giả ba chiều (BOGUS) vào trường phân tích toàn cầu GFS để xác định cấu trúc xoáy ban đầu.
    • Tách biệt quy trình thành 2 thành phần: nuôi nhiễu xoáy (trong phạm vi bán kính bão) và nuôi nhiễu môi trường (khu vực bên ngoài hoàn lưu bão). Chu kỳ nuôi thực hiện trong 24 giờ trước thời điểm bắt đầu dự báo ($T_{00}$), chia làm 4 bước tích phân 6 giờ: từ -24h đến -18h, -12h, -06h và 00h.
    • Chuẩn hóa biên độ nhiễu môi trường và nhiễu xoáy dựa trên hàm năng lượng.
    • Trực giao hóa các véc tơ nhiễu để đảm bảo tính độc lập thông tin giữa các thành phần tổ hợp.
  4. Kỹ thuật tổ hợp kết quả:
    • Phương pháp trung bình đơn giản (Ensemble Mean) số học của các thành phần dự báo.
    • Phương pháp Siêu tổ hợp (Superensemble) dựa trên lý thuyết của Krishnamurti và Jordan (2005), sử dụng kỹ thuật hồi quy tuyến tính đa biến để tính toán hệ số trọng số tối ưu cho kinh độ và vĩ độ tâm bão ở từng hạn dự báo (24, 48, 72, 96, 120 giờ).
    • Dự báo xác suất quỹ đạo bão dựa trên bán kính vòng tròn độ tán tổ hợp của các thành phần.
  5. Phương pháp đánh giá sai số: Tính toán 3 loại sai số theo Best Track của JMA:
    • Sai số khoảng cách vị trí tâm bão (km)
    • Sai số dọc theo hướng di chuyển (along-track error, km)
    • Sai số ngang vuông góc với hướng di chuyển (cross-track error, km)

Nguồn tư liệu và dữ liệu:

  • Dữ liệu khí quyển toàn cầu: Các trường phân tích và dự báo GFS (Global Forecast System) của Trung tâm Quốc gia Dự báo Môi trường Hoa Kỳ (NCEP).
  • Dữ liệu quỹ đạo bão thực tế: Số liệu Best Track và các bản đồ dự báo của Cơ quan Khí tượng Nhật Bản (JMA).
  • Số liệu thống kê bão lịch sử: Bộ số liệu bão hoạt động trên Biển Đông và đổ bộ Việt Nam thời kỳ 1959–2008 (Đinh Văn Ưu, 2009; Bộ Tài nguyên và Môi trường, 2012) và 1961–2010 (Nguyễn Văn Thắng, 2010).

Nội dung chính theo từng chương

CHƯƠNG 1: Tổng quan về dự báo tổ hợp và dự báo tổ hợp quỹ đạo bão

Chương này trình bày cơ sở lý luận của phương pháp dự báo tổ hợp, lịch sử phát triển các kỹ thuật tạo nhiễu trường ban đầu và điều kiện biên, cùng tình hình hoạt động của bão tại khu vực Tây Bắc Thái Bình Dương, Biển Đông và Việt Nam. Tác giả phân tích cấu trúc của các hệ thống dự báo tổ hợp nghiệp vụ trên thế giới, tiêu biểu như hệ thống SREF của NCEP (Hoa Kỳ) gồm 21 thành phần được tổng hợp trong văn bản:

Mô hình Điều kiện ban đầu / Nhiễu ban đầu Số thành phần Độ phân giải (km) Hạn dự báo (giờ) Nhiễu biên (LBC)
Eta_BMJ Ndas / nhiễu nuôi từ mô hình khu vực 3 (ctl1, n1, p1) 32 87 GFS/GEFS
Eta_KF Ndas / nhiễu nuôi từ mô hình khu vực 3 (ctl1, n2, p2) 32 87 GFS/GEFS
RSM_SAS GFS / nhiễu nuôi từ mô hình khu vực 3 (ctl, n1, p1) 32 87 GFS/GEFS
RSM_RAS GFS / nhiễu nuôi từ mô hình khu vực 2 (n2, p2) 32 87 GFS/GEFS
WRF_NMM GFS / nhiễu ET từ mô hình toàn cầu 5 (ctl, n1, p1, n2, p2) 32 87 GFS/GEFS
WRF_ARW GFS / nhiễu ET từ mô hình toàn cầu 5 (ctl, n1, p1, n2, p2) 35 87 GFS/GEFS

Về đặc điểm khí hậu bão, tác giả tổng hợp các bảng số liệu lịch sử phản ánh tần suất hoạt động và quy luật di chuyển của bão:

Bảng 1: Số lượng bão trung bình theo từng thập kỷ (Đinh Văn Ưu, 2009)

Thời gian 1959–1968 1969–1978 1979–1988 1989–1998 1999–2008 Trung bình
Tây Bắc Thái Bình Dương (cơn/năm) 35,1 28,0 27,3 34,2 29,5 30,8
Biển Đông (cơn/năm) 12,4 11,6 10,6 13,5 11,4 11,9

Trong số 11,9 cơn/năm trên Biển Đông, có 7,6 cơn hình thành từ Tây Bắc Thái Bình Dương di chuyển vào và 4,3 cơn hình thành trực tiếp trên Biển Đông.

Bảng 2: Phân bố bão trên Biển Đông và bão đổ bộ Việt Nam thời kỳ 1961–2010 (Nguyễn Văn Thắng, 2010)

Chỉ tiêu / Tháng I II III IV V VI VII VIII IX X XI XII TB năm
Bão trên Biển Đông - - - - 0,62 1,08 1,88 2,02 2,22 1,90 1,48 0,66 -
Bão đổ bộ Việt Nam 0,12 0,02 0,08 0,06 0,26 0,60 0,88 1,08 1,60 1,36 1,16 0,40 7,62

Văn bản ghi nhận quy luật dịch chuyển vùng đổ bộ: từ vĩ độ 20°N trở lên có xu hướng giảm, từ 15°N–20°N ít thay đổi, từ 15°N trở xuống có xu hướng gia tăng, cho thấy hoạt động bão có xu thế dịch chuyển dần về phía Nam.

Chương 1 cũng tổng kết kinh nghiệm ứng dụng EPS tại HKO (giảm 10–15% sai số hạn 24–72 giờ), hệ thống TMEPS của Hàn Quốc (12 nguồn số liệu), hệ thống tổ hợp bão NMC Bắc Kinh (15 thành phần hạn 120 giờ) và GRAPES_TCM Thượng Hải (9 thành phần từ WRF).

CHƯƠNG 2: Phương pháp nuôi nhiễu khi có bão và lựa chọn các tham số cho mô hình RAMS để dự báo tổ hợp quỹ đạo bão ở Việt Nam

Chương 2 tập trung vào giải pháp kỹ thuật cốt lõi do tác giả đề xuất: phương án nuôi nhiễu môi trường và nhiễu xoáy riêng biệt khi xuất hiện bão.

Nội dung kỹ thuật gồm các bước cụ thể:

  1. Tạo nhân ban đầu: Sử dụng phân bố chuẩn Gauss để tạo 6 trường nhân ban đầu (D1 đến D6) với độ lệch vận tốc gió 3 m/s và nhiệt độ 1°C tại các mực khí áp.
  2. Cài xoáy giả cho trường GFS: Áp dụng phương pháp cài xoáy giả 3 chiều (Bogus) nhằm khắc phục nhược điểm định vị tâm và cường độ thô của trường phân tích toàn cầu.
  3. Quy trình nuôi nhiễu trong 24 giờ: Tiến hành nuôi nhiễu trước thời điểm làm dự báo $T_{00}$ qua các mốc -24h, -18h, -12h, -06h và 00h. Quá trình nuôi được phân tách độc lập giữa nhiễu xoáy (khu vực hoàn lưu bão) và nhiễu môi trường (khu vực trường nền bên ngoài).
  4. Chuẩn hóa và trực giao hóa:
    • Chuẩn hóa biên độ: Đưa biên độ của các véc tơ nuôi về mức sai số phân tích thực tế thông qua tích phân năng lượng toàn phần (gió, nhiệt độ thế vị, độ ẩm, khí áp bề mặt).
    • Trực giao hóa: Áp dụng phương pháp trực giao hàm trực giao tự nhiên (EOF) hoặc ma trận quay để loại bỏ tính tương quan giữa các véc tơ nuôi, đảm bảo các thành phần nhiễu độc lập tối đa và bao quát được không gian bất định.
  5. Cấu hình mô hình RAMS: Thiết lập miền tính bao trùm khu vực Tây Bắc Thái Bình Dương và Biển Đông; tích hợp 3 sơ đồ tham số hóa đối lưu (KUO, KF, KFCT). Kết hợp 12 thành phần nhiễu nuôi với 3 sơ đồ đối lưu và thành phần kiểm chứng tạo ra tập hợp 39 thành phần dự báo.
  6. Phương pháp xác định tâm bão và tổ hợp: Xây dựng công cụ dò tìm tâm bão tự động dựa trên cực tiểu khí áp và hoàn lưu gió xoáy. Thiết lập 2 thuật toán tổ hợp:
    • Trung bình đơn giản số học của các thành phần.
    • Siêu tổ hợp (Superensemble): Ước lượng các hệ số hồi quy đa biến cho tọa độ kinh độ và vĩ độ bão theo từng hạn dự báo 24, 48, 72, 96 và 120 giờ.
    • Dự báo xác suất: Xác định bán kính vòng tròn độ tán tổ hợp để thể hiện xác suất xuất hiện tâm bão.

CHƯƠNG 3: Dự báo tổ hợp quỹ đạo bão ở Việt Nam dựa trên phương pháp nuôi nhiễu phát triển

Chương 3 trình bày các kết quả thử nghiệm mô phỏng và đánh giá độ chính xác của hệ thống:

  1. Thử nghiệm trên cơn bão Washi (tháng 12/2011) hạn 120 giờ:
    • Phân tích quá trình phát triển của nhiễu môi trường, nhiễu xoáy và kết hợp cả hai trong 24 giờ nuôi. Kết quả cho thấy các trường nhiễu phát triển rõ nét về cấu trúc thẳng đứng và phương ngang.
    • Đánh giá sai số khoảng cách của 12 thành phần dự báo: Trường hợp kết hợp đồng thời nhiễu môi trường và nhiễu xoáy cho sai số khoảng cách trung bình tổ hợp nhỏ hơn rõ rệt so với trường hợp chỉ nuôi riêng rẽ nhiễu môi trường hoặc chỉ nuôi riêng rẽ nhiễu xoáy ở tất cả các hạn từ 24 đến 120 giờ.
    • Vai trò của đa cơ chế vật lý: Khi kết hợp 12 thành phần nhiễu với từng sơ đồ đối lưu đơn lẻ (KUO, KF, KFCT), sai số có sự khác biệt; việc kết hợp cả 3 sơ đồ đối lưu (tạo thành 39 thành phần) mang lại độ tán tổ hợp tốt hơn, bao trùm được quỹ đạo thực tế của bão.
  2. Thử nghiệm mở rộng và xây dựng phương trình siêu tổ hợp:
    • Thử nghiệm trên 199 trường hợp bão giai đoạn 2009–2013 ở các hạn dự báo.
    • Tuyển chọn thành phần: Tác giả phân tích tần suất xuất hiện của các thành phần trong 40 phương trình hồi quy dự báo kinh độ và vĩ độ tâm bão. Đồ thị tốc độ suy giảm sai số cho thấy khi tăng số lượng thành phần tham gia, sai số khoảng cách giảm nhanh từ 1 đến khoảng 20–25 thành phần; từ thành phần thứ 25 trở đi, tốc độ suy giảm sai số bắt đầu bão hòa. Do đó, luận án lựa chọn 25 thành phần tối ưu để xây dựng phương trình dự báo siêu tổ hợp.
    • Bảng phương trình hồi quy siêu tổ hợp được xác lập riêng biệt cho các hạn 24, 48, 72, 96 và 120 giờ.
  3. Đánh giá sai số trên bộ số liệu độc lập năm 2012 và 2013:
    • So sánh giữa dự báo siêu tổ hợp và trung bình tổ hợp 25 thành phần: Phương pháp siêu tổ hợp đạt độ chính xác cao hơn, làm giảm rõ rệt sai số vị trí tâm bão, sai số dọc và sai số ngang ở hầu hết các hạn dự báo từ 24 đến 120 giờ.
    • Thử nghiệm trên các cơn bão có tính chất quỹ đạo điển hình:
      • Bão Megi (18/10/2010): Siêu tổ hợp mô phỏng sát quỹ đạo thực khi bão đổi hướng.
      • Bão đôi Nesat và Nalgae (27/10/2011): Dự báo mô phỏng thành công sự tương tác giữa hai cơn bão gần nhau.
      • Bão Gaemi (30/09/2012): Xử lý tốt trường hợp bão có hướng di chuyển phức tạp, dị thường trên Biển Đông.
      • Bão Bopha (29/11/2012): Cơn bão mạnh hình thành ở vĩ độ thấp.
  4. Dự báo xác suất và chỉ số tin cậy:
    • Phân tích mối quan hệ giữa độ tán tổ hợp tích lũy 6 giờ và sai số vị trí thực tế.
    • Xây dựng bảng độ tin cậy của dự báo theo 3 mức A, B, C: Mức A tương ứng với độ tán tổ hợp nhỏ và sai số vị trí thấp (độ tin cậy cao); Mức B là trung bình; Mức C tương ứng với độ tán lớn và sai số vị trí cao (độ tin cậy thấp).
  5. Đề xuất quy trình công nghệ: Luận án sơ đồ hóa quy trình dự báo bão hạn 5 ngày hoàn chỉnh từ khâu tiếp nhận số liệu GFS, cài xoáy giả, chạy chu trình nuôi nhiễu 24 giờ, tích hợp mô hình RAMS đa cơ chế, chạy hồi quy siêu tổ hợp, trích xuất sai số và xuất bản đồ xác suất.

Kết quả và những đóng góp mới

Theo công bố của tác giả, luận án đạt được các kết quả và đóng góp mới sau:

Về mặt lý luận và phương pháp:

  • Áp dụng và cải tiến thành công phương pháp nuôi nhiễu tạo trường ban đầu trong điều kiện có bão cho khu vực Tây Bắc Thái Bình Dương và Biển Đông. Điểm cốt lõi là việc phân tách, chuẩn hóa và trực giao hóa độc lập giữa nhiễu xoáy (mô tả cấu trúc bão có kết hợp cài xoáy giả 3 chiều) và nhiễu môi trường (mô tả bất định trường nền).
  • Chứng minh được tính hiệu quả của việc kết hợp đa cơ chế tham số hóa đối lưu (KUO, KF, KFCT) với các thành phần nuôi nhiễu trong mô hình RAMS nhằm tạo độ tán tổ hợp tối ưu.
  • Xây dựng thành công phương pháp siêu tổ hợp (Superensemble) gồm 25 thành phần được tuyển chọn qua phân tích tốc độ suy giảm sai số và tần suất xuất hiện trong các phương trình hồi quy đa biến cho kinh vĩ độ tâm bão.

Về mặt thực tiễn và ứng dụng:

  • Xây dựng thành công hệ thống và quy trình công nghệ dự báo tổ hợp quỹ đạo bão ảnh hưởng đến Việt Nam với hạn dự báo kéo dài đến 5 ngày (120 giờ) dựa trên phương pháp nuôi nhiễu phát triển.
  • Cung cấp đồng thời sản phẩm dự báo tất định (quỹ đạo từ siêu tổ hợp và trung bình tổ hợp) cùng sản phẩm dự báo xác suất (bán kính vòng tròn độ tán tổ hợp và chỉ số tin cậy A, B, C), giúp lượng hóa độ bất định của dự báo bão.
  • Kết quả kiểm chứng trên 199 trường hợp trong 5 năm (2009–2013) khẳng định tính khả thi và độ tin cậy của quy trình dự báo trong điều kiện thực tế tại Việt Nam.

Kiến nghị và giải pháp đề xuất:

  • Đề xuất đưa quy trình dự báo tổ hợp bão hạn 5 ngày vào vận hành thử nghiệm nghiệp vụ tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia.
  • Kết hợp sử dụng thông tin độ tán tổ hợp và chỉ số tin cậy (A, B, C) trong việc ban hành các bản tin cảnh báo bão sớm phục vụ Ban Chỉ đạo Quốc gia về Phòng chống thiên tai.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu:

  • Điều kiện biên ngang (LBC) của mô hình RAMS trong nghiên cứu vẫn sử dụng từ một nguồn dự báo toàn cầu duy nhất (mô hình GFS của NCEP) chứ chưa tạo được nhiễu trên biên từ một hệ thống tổ hợp toàn cầu chuyên biệt hoặc đa mô hình.
  • Mẫu số liệu dùng để xây dựng phương trình hồi quy siêu tổ hợp được thực hiện trên 199 trường hợp bão giai đoạn 2009–2013; chuỗi số liệu này cần tiếp tục được kéo dài để tăng độ ổn định của các trọng số hồi quy.
  • Luận án mới tập trung giải quyết bài toán dự báo quỹ đạo di chuyển của tâm bão, chưa giải quyết sâu bài toán dự báo tổ hợp cường độ bão, cấu trúc trường gió mạnh và lượng mưa lớn do bão gây ra.

Hướng nghiên cứu tiếp theo:

  • Nghiên cứu áp dụng kỹ thuật gây nhiễu trên điều kiện biên ngang (LBC) từ các hệ thống tổ hợp toàn cầu để hạn chế sự suy giảm độ tán tổ hợp ở các hạn dự báo dài.
  • Mở rộng hệ thống dự báo tổ hợp bằng phương pháp nuôi nhiễu để dự báo cường độ bão, bán kính gió mạnh và phân bố lượng mưa bão.
  • Tiếp tục cập nhật chuỗi số liệu các cơn bão mới để hiệu chỉnh, cập nhật bộ hệ số phương trình siêu tổ hợp.

Giá trị tham khảo

Luận án có giá trị tham khảo chuyên môn đối với các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và cán bộ nghiên cứu chuyên ngành Khí tượng, Khí hậu học, Dự báo thời tiết số trị: Luận án là tài liệu tham khảo chi tiết về thuật toán nuôi nhiễu (BV), kỹ thuật chuẩn hóa năng lượng, trực giao hóa véc tơ nuôi, kỹ thuật cài xoáy giả 3 chiều và phương pháp hồi quy siêu tổ hợp (Krishnamurti). Phần đáng chú ý nhất là phương pháp luận phân tách nhiễu xoáy - nhiễu môi trường ở Chương 2 và các phương trình hồi quy siêu tổ hợp ở Chương 3.
  • Giảng viên đại học: Sử dụng làm tài liệu giảng dạy cho các học phần chuyên đề về mô hình hóa khí quyển, dự báo thời tiết hạn vừa và lý thuyết dự báo tổ hợp (EPS).
  • Các cơ quan dự báo khí tượng nghiệp vụ và quản lý thiên tai: Bản mô tả quy trình công nghệ dự báo 5 ngày, thuật toán xác định tâm bão tự động và cách phân cấp chỉ số tin cậy (A, B, C) cung cấp cơ sở kỹ thuật thực tế để ứng dụng vào công tác cảnh báo bão sớm tại Việt Nam.

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

1. Điểm khác biệt cơ bản giữa kỹ thuật nuôi nhiễu của luận án so với kỹ thuật nuôi nhiễu (BV) truyền thống của Toth và Kalnay là gì?
Kỹ thuật nuôi nhiễu truyền thống của Toth và Kalnay (1993, 1997) được áp dụng trên toàn bộ trường khí quyển quy mô lớn mà không phân biệt khu vực có bão. Luận án của tác giả Công Thanh đã cải tiến bằng cách phân tách độc lập thành 2 thành phần: nhiễu xoáy (được kết hợp cài xoáy giả 3 chiều để mô tả cấu trúc bất đối xứng của bão) và nhiễu môi trường (mô tả bất định của trường nền bên ngoài). Hai dạng nhiễu này được nuôi, chuẩn hóa năng lượng và trực giao hóa độc lập trước khi tích hợp vào trường ban đầu.

2. Hệ thống dự báo tổ hợp trong luận án gồm bao nhiêu thành phần và tại sao tác giả lựa chọn 25 thành phần cho phương án siêu tổ hợp?
Hệ thống ban đầu tạo ra 39 thành phần dự báo (kết hợp từ 12 thành phần nhiễu nuôi và 1 thành phần kiểm chứng qua 3 sơ đồ đối lưu KUO, KF, KFCT: $13 \times 3 = 39$). Khi phân tích đồ thị tốc độ suy giảm sai số khoảng cách theo số lượng thành phần tham gia, tác giả nhận thấy sai số giảm rất nhanh từ 1 đến khoảng 20–25 thành phần; từ thành phần thứ 25 trở đi, tốc độ suy giảm bắt đầu bão hòa. Đồng thời qua phân tích tần suất xuất hiện trong 40 phương trình hồi quy dự báo kinh vĩ độ, tác giả đã chọn ra 25 thành phần tối ưu nhất để xây dựng phương trình siêu tổ hợp.

3. Mô hình RAMS trong luận án đã sử dụng những sơ đồ tham số hóa đối lưu nào và vai trò của chúng ra sao?
Mô hình RAMS sử dụng 3 sơ đồ tham số hóa đối lưu: KUO, Kain-Fritsch (KF), và Kain-Fritsch cải tiến bởi Nguyễn Minh Trường (KFCT). Việc thay đổi đồng thời 3 sơ đồ đối lưu này đóng vai trò mô tả sự bất định về mặt vật lý của mô hình, giúp các thành phần tổ hợp phân nhánh tốt hơn và tăng độ tán tổ hợp so với việc chỉ sử dụng một sơ đồ đối lưu duy nhất.

4. Chỉ số độ tin cậy (A, B, C) trong dự báo xác suất quỹ đạo bão của luận án được xác định như thế nào?
Chỉ số độ tin cậy được xác định dựa trên mối quan hệ định lượng giữa độ tán tổ hợp tích lũy (tổng độ tán tổ hợp 6 giờ từ hạn 6 đến 120 giờ) và sai số vị trí thực tế của bão. Cụ thể: Mức A biểu thị trường hợp độ tán tổ hợp nhỏ tương ứng với sai số vị trí dự báo thấp (độ tin cậy cao); Mức B biểu thị mức độ tán và sai số trung bình; Mức C biểu thị trường hợp độ tán tổ hợp lớn tương ứng với sai số vị trí cao (độ tin cậy thấp).

5. Kết quả thử nghiệm trên cơn bão Washi (2011) đã chứng minh điều gì về vai trò của việc kết hợp nhiễu xoáy và nhiễu môi trường?
Kết quả mô phỏng cơn bão Washi ở hạn dự báo 120 giờ chỉ ra rằng sai số khoảng cách trung bình tổ hợp của trường hợp kết hợp đồng thời nhiễu môi trường và nhiễu xoáy luôn thấp hơn đáng kể so với hai trường hợp chỉ nuôi riêng rẽ nhiễu môi trường hoặc chỉ nuôi riêng rẽ nhiễu xoáy, khẳng định tính đúng đắn của phương án nuôi nhiễu kép do tác giả đề xuất.

Kết luận

Luận án tiến sĩ của tác giả Công Thanh đã giải quyết bài toán khoa học về kéo dài thời hạn dự báo quỹ đạo bão ảnh hưởng đến Việt Nam lên hạn 5 ngày bằng phương pháp dự báo tổ hợp trên mô hình khu vực RAMS. Đóng góp nổi bật của công trình là cải tiến kỹ thuật nuôi nhiễu thông qua việc phân tách, chuẩn hóa và trực giao hóa đồng thời nhiễu xoáy (kết hợp cài xoáy giả ba chiều) và nhiễu môi trường, kết hợp đa cơ chế đối lưu (KUO, KF, KFCT). Thử nghiệm trên 199 trường hợp bão giai đoạn 2009–2013 cho thấy phương pháp siêu tổ hợp 25 thành phần cùng dự báo xác suất phân cấp theo độ tin cậy (A, B, C) giúp giảm thiểu rõ rệt sai số vị trí tâm bão so với phương pháp dự báo đơn lẻ. Công trình đã đề xuất được một quy trình dự báo hoàn chỉnh, đáp ứng yêu cầu nâng cao chất lượng cảnh báo sớm thiên tai bão tại Việt Nam.