Tổng quan nghiên cứu

Khu vực miền Trung và Tây Nguyên là địa bàn trọng điểm thường xuyên hứng chịu các đợt mưa lớn cực đoan, gây thiệt hại nghiêm trọng về người và tài sản với hàng chục đợt thiên tai mỗi năm. Địa hình hiểm trở kéo dài theo dãy Trường Sơn kết hợp với 5 dạng hình thế thời tiết synop điển hình (xoáy thuận nhiệt đới, dải hội tụ nhiệt đới ITCZ, không khí lạnh tăng cường, nhiễu động gió Đông và hội tụ gió trên cao) tạo nên chế độ mưa vô cùng phức tạp. Đặc biệt, trường mưa mang bản chất phi tuyến, không cố định theo không gian và thời gian, khiến công tác dự báo định lượng mưa luôn là bài toán thách thức hàng đầu trong ngành khí tượng học.

Nhằm giải quyết bài toán cấp thiết này, luận văn tập trung đánh giá định lượng và xác thực kỹ năng dự báo mưa lớn tại các hạn dự báo 24 giờ, 48 giờ và 72 giờ của 3 mô hình dự báo số trị khu vực đang vận hành nghiệp vụ gồm HRM-gfs, WRFARW-gfs và WRFNMM-gfs. Nghiên cứu được triển khai trên phạm vi toàn diện gồm 4 tiểu vùng khí hậu: Bắc Trung Bộ, Trung Trung Bộ, Nam Trung Bộ và Tây Nguyên dựa trên chuỗi số liệu quan trắc liên tục trong giai đoạn 3 năm (2008 - 2010) tại 64 trạm khí tượng bề mặt.

Ý nghĩa khoa học và thực tiễn của nghiên cứu thể hiện qua việc cung cấp bức tranh định lượng chi tiết về sai số hệ thống, xác suất phát hiện và độ tin cậy của từng mô hình theo từng ngưỡng mưa (mưa vừa 16-50 mm/24h, mưa to 50-100 mm/24h và mưa rất to trên 100 mm/24h). Kết quả thẩm tra đóng vai trò làm căn cứ chuẩn mực giúp các dự báo viên tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia lựa chọn mô hình tối ưu cho từng khu vực, nâng cao xác suất phát hiện mưa lớn lên khoảng 10-15% và giảm thiểu tỷ lệ báo động sai, phục vụ đắc lực cho công tác vận hành hồ chứa và phòng chống lũ quét, sạt lở đất.

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 dựa trên nền tảng lý thuyết thủy động lực học khí quyển và hệ phương trình nguyên thủy mô tả các quá trình đối lưu, nhiệt động lực học và truyền bức xạ trong khí quyển. Hai lý thuyết cốt lõi được áp dụng gồm: lý thuyết mô hình hóa số trị thời tiết quy mô vừa (Mesoscale Numerical Weather Prediction) và lý thuyết thẩm tra đánh giá dự báo khí tượng (Meteorological Forecast Verification) theo tiêu chuẩn của Tổ chức Khí tượng Thế giới (WMO).

Mô hình nghiên cứu so sánh 3 cấu hình động lực học tiên tiến:

  • Mô hình HRM (phiên bản 8) sử dụng hệ phương trình thủy tĩnh, sơ đồ tham số hóa đối lưu Tiedtke, bức xạ Ritter/Geleyn và lưới sai phân Arakawa C với độ phân giải ngang 0.125 độ (khoảng 14 km), 40 mực thẳng đứng, tích phân theo bước thời gian 60 giây.
  • Mô hình WRF-ARW (phiên bản 3.3) tích hợp sơ đồ tham số hóa đối lưu Kain-Fritsch, vi vật lý mây WSM-3, bức xạ RRTM/Dudhia, hệ tọa độ theo áp suất với độ phân giải 13.75 km và 40 mực khí quyển.
  • Mô hình WRF-NMM vận hành với độ phân giải ngang 13.75 km trên lưới Arakawa E/C cải tiến, 40 mực thẳng đứng và bước tích phân 60 giây.

Các khái niệm chuyên ngành then chốt bao gồm: dự báo biến liên tục (lượng mưa tích lũy mm/24h), dự báo pha nhị phân (có/không xuất hiện mưa theo ngưỡng), hạn dự báo (lead time 24h, 48h, 72h), điều kiện biên toàn cầu (GFS 0.5 độ cập nhật 6 giờ một lần) và bảng liên hợp thống kê 2x2.

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

Nguồn dữ liệu thực nghiệm được thu thập đồng bộ từ 64 trạm quan trắc khí tượng bề mặt phân bố khắp khu vực miền Trung và Tây Nguyên, bao gồm: 20 trạm tại Bắc Trung Bộ, 15 trạm tại Trung Trung Bộ, 12 trạm tại Nam Trung Bộ và 17 trạm tại Tây Nguyên. Toàn bộ tập dữ liệu dự báo dạng lưới từ 3 mô hình được trích xuất và nội suy song tuyến tính (Bilinear Interpolation) về đúng tọa độ địa lý của 64 trạm quan trắc tương ứng.

Phương pháp chọn mẫu là chọn mẫu toàn bộ theo chuỗi thời gian liên tục 36 tháng (giai đoạn 2008 - 2010), lọc ra tất cả các đợt mưa lớn diện rộng và cục bộ điển hình do bão, áp thấp nhiệt đới, không khí lạnh và dải hội tụ nhiệt đới gây ra. Dung lượng mẫu gồm hàng nghìn cặp số liệu dự báo - quan trắc được phân lớp nghiêm ngặt theo 3 cấp độ mưa chuẩn quốc tế: mưa vừa (16 - 50 mm/24h), mưa to (51 - 100 mm/24h) và mưa rất to (trên 100 mm/24h).

Lý do lựa chọn kết hợp hai nhóm phương pháp phân tích định lượng và phân hạng là do tính chất đặc thù của mưa: vừa có giá trị định lượng liên tục vừa có tính gián đoạn không gian cao:

  • Phương pháp đánh giá biến liên tục: Áp dụng các chỉ số sai số trung bình (ME), sai số tuyệt đối trung bình (MAE), sai số bình phương trung bình quân phương (RMSE) và hệ số tương quan Pearson (r) để lượng hóa độ lệch biên độ và độ phân tán của mô hình.
  • Phương pháp đánh giá theo ngưỡng (Contingency Matrix): Sử dụng tỷ số vùng dự báo trên quan trắc (BIAS/FBI), xác suất phát hiện hiện tượng đúng (POD), tỷ lệ báo động sai (FAR) và điểm số thành công hợp lý (ETS) nhằm đo lường khả năng bắt trúng các hiện tượng mưa cực đoan mà không bị thiên lệch bởi các ngày không mưa.

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 HRM-gfs thể hiện kỹ năng dự báo tổng thể vượt trội nhất trên toàn bộ khu vực miền Trung và Tây Nguyên. Ở hạn dự báo 24 giờ và 48 giờ, hệ số tương quan tuyến tính (r) của mô hình HRM đạt mức 0.40 - 0.45, cao hơn đáng kể so với WRFARW-gfs (đạt khoảng 0.35 - 0.38) và WRFNMM-gfs (đạt dưới 0.30). Đặc biệt tại khu vực Nam Trung Bộ, hệ số tương quan của mô hình đạt mức đỉnh xấp xỉ 0.50.

Thứ hai, nghiên cứu đã phát hiện rõ nét quy luật sai số hệ thống đối lập giữa các mô hình:

  • Mô hình HRM-gfs và WRFARW-gfs luôn cho giá trị sai số trung bình dương (ME > 0) trên hầu hết các hạn dự báo, phản ánh xu hướng dự báo thiên cao (dự báo lượng mưa và diện mưa lớn hơn thực tế từ 10% đến 25%).
  • Ngược lại, mô hình WRFNMM-gfs liên tục ghi nhận giá trị ME âm (ME < 0) trên cả 3 hạn dự báo 24h, 48h, 72h, cho thấy mô hình này có thiên hướng dự báo thiên thấp, thường xuyên đánh giá thấp lượng mưa thực tế đo được tại các trạm.

Thứ ba, kỹ năng dự báo có sự phân hóa sâu sắc theo cấp độ cường độ mưa:

  • Ở ngưỡng mưa vừa (16 - 50 mm/24h): HRM-gfs cho tỷ số diện mưa lý tưởng nhất với BIAS xấp xỉ 1.0; trong khi WRFARW dự báo khống diện tích với BIAS lên tới 1.4 và WRFNMM dự báo hụt với BIAS khoảng 0.8.
  • Ở ngưỡng mưa to (50 - 100 mm/24h): WRFARW-gfs lại chiếm ưu thế rõ rệt với chỉ số BIAS đạt sát ngưỡng tối ưu trong khoảng 0.80 - 0.90, trong khi HRM giảm xuống 0.60 - 0.70 và WRFNMM chỉ đạt xấp xỉ 0.50.
  • Ở ngưỡng mưa rất to (trên 100 mm/24h): Cả 3 mô hình đều suy giảm kỹ năng rõ rệt, diện tích dự báo nhỏ hơn nhiều so với thực tế (BIAS của WRFNMM tụt xuống xấp xỉ 0.20). Tuy nhiên, HRM-gfs vẫn duy trì xác suất phát hiện POD cao nhất và tỷ lệ báo động sai FAR thấp nhất (khoảng 0.45 - 0.50) ở hạn 24 giờ.

Thứ tư, sai số mô hình tăng mạnh theo hạn dự báo. Khi tăng thời gian dự báo từ 24h lên 48h và 72h, chỉ số sai số tuyệt đối MAE và sai số quân phương RMSE của cả 3 mô hình đều tăng trung bình từ 20% đến 35%, đồng thời xác suất phát hiện POD giảm đều từ 15% đến 25%.

Thảo luận kết quả

Nguyên nhân dẫn đến sự khác biệt về kỹ năng giữa các mô hình bắt nguồn chủ yếu từ cách thức tham số hóa vật lý đối lưu và xử lý điều kiện biên địa hình. Sơ đồ đối lưu Tiedtke trong mô hình HRM chứng minh khả năng mô phỏng cân bằng ẩm đối lưu nhiệt đới ổn định hơn, giúp hạn chế các sai số bùng nổ cục bộ. Trong khi đó, sơ đồ Kain-Fritsch trong WRF-ARW kích hoạt cơ chế đối lưu mạnh mẽ hơn, dẫn đến việc mở rộng diện tích mưa vượt mức thực tế ở các ngưỡng mưa vừa và nhỏ nhưng lại bắt trúng tâm mưa cường độ lớn (50-100 mm/24h) tốt hơn khi có sự tương tác của bão và không khí lạnh.

Đối với khu vực địa hình đèo dốc hiểm trở như Bắc Trung Bộ và Trung Trung Bộ, hiệu ứng nâng địa hình chắn gió phía Đông dãy Trường Sơn gây ra sai số vị trí tâm mưa cục bộ. Dữ liệu nghiên cứu có thể được minh họa trực quan thông qua các toán đồ tụ điểm phân tán 45 độ kết hợp biểu đồ cột so sánh bộ chỉ số đánh giá theo từng dải cường độ mưa (16-30, 30-50, 50-70, 70-100 và trên 100 mm/24h). Cách trình bày này làm nổi bật mối liên hệ giữa bước thời gian tích phân và độ phân tán sai số, khẳng định tính tất yếu của việc hiệu chỉnh địa hình cục bộ trước khi đưa sản phẩm số vào bản tin dự báo tác động.

So sánh với các nghiên cứu trước đây về mô hình MM5 tại Việt Nam, kết quả của luận văn khẳng định bước tiến vượt bậc của thế hệ mô hình HRM và WRF với hệ số tương quan nâng từ mức 0.30 lên 0.45 - 0.50, cung cấp nguồn dữ liệu đầu vào đáng tin cậy hơn cho các mô hình thủy văn tính toán dòng chảy lũ.

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

Nhằm tối ưu hóa hiệu quả ứng dụng các mô hình số trị trong công tác dự báo thời tiết nguy hiểm tại miền Trung và Tây Nguyên, 4 giải pháp trọng tâm được đề xuất:

  • Thiết lập quy trình khai thác tổ hợp đa mô hình phân vùng: Ban hành hướng dẫn nghiệp vụ yêu cầu các dự báo viên ưu tiên tham khảo mô hình HRM-gfs đối với các đợt dự báo mưa diện rộng và mưa cực lớn trên 100 mm/24h; đồng thời kết hợp đối chiếu mô hình WRFARW-gfs khi nhận định các tâm mưa cường độ 50 - 100 mm/24h tại Trung Trung Bộ. Mục tiêu nâng xác suất phát hiện POD tổng thể thêm 12 - 18% trong giai đoạn 2024 - 2026, do Trung tâm Dự báo Khí tượng Thủy văn Quốc gia chủ trì thực hiện.
  • Triển khai công nghệ đồng hóa dữ liệu quan trắc địa phương độ phân giải cao: Tích hợp hệ thống đồng hóa 3DVAR/EnKF thu nhận dữ liệu từ mạng lưới radar thời tiết ven biển và ảnh vệ tinh đa phổ vào mô hình WRF nhằm chuẩn hóa trường ban đầu. Mục tiêu giảm sai số tuyệt đối MAE từ mức 25 - 30 mm xuống dưới 18 mm/24h trong lộ trình 24 tháng, do Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu phối hợp thực hiện.
  • Xây dựng mô đun hậu xử lý thống kê khử sai số hệ thống: Áp dụng phương pháp hiệu chỉnh sai số MOS (Model Output Statistics) kết hợp thuật toán học máy để tự động triệt tiêu độ lệch dương (ME > 0) của HRM và độ lệch âm (ME < 0) của WRFNMM tại từng trạm quan trắc. Đích đến là nâng hệ số tương quan r tại Bắc Trung Bộ từ mức 0.20 hiện tại lên trên 0.55 trước năm 2026, do Phòng Dự báo Số trị đảm trách.
  • Tăng cường độ phân giải lưới tính toán cho khu vực sườn đón gió: Nâng cấp độ phân giải mô hình WRF từ 13.75 km xuống dưới 3 km bao phủ dải ven biển và vùng núi cao Tây Nguyên, nhằm mô phỏng chính xác hoàn lưu quy mô vừa do đứt gãy địa hình gây ra, cắt giảm tỷ lệ báo động sai FAR xuống dưới 0.35 vào năm 2027.

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

Nội dung luận văn mang giá trị học thuật sâu sắc và tính ứng dụng thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  • Dự báo viên khí tượng tại Trung tâm Dự báo KTTV Quốc gia và các Đài KTTV khu vực: Luận văn cung cấp cẩm nang tra cứu định lượng chính xác về ưu/nhược điểm, sai số bias của từng mô hình trên 64 trạm cụ thể, giúp rút ngắn 30% thời gian phân tích synop và nâng cao độ chính xác khi ban hành bản tin cảnh báo thiên tai sớm 24-72 giờ.
  • Các chuyên gia nghiên cứu và phát triển mô hình số trị (NWP Modelers): Cung cấp bộ dữ liệu kiểm chứng thực nghiệm chi tiết về tác động của các sơ đồ tham số hóa đối lưu Tiedtke, Kain-Fritsch đối với khí hậu nhiệt đới gió mùa, làm cơ sở tinh chỉnh hệ số vật lý cho 40 tầng khí quyển.
  • Ban Chỉ huy Phòng thủ dân sự, Phòng chống thiên tai và Tìm kiếm cứu nạn các cấp: Giúp các nhà quản lý thấu hiểu bản chất độ tin cậy và biên độ sai số của các bản tin dự báo mưa hạn 1 đến 3 ngày, từ đó tối ưu hóa phương án điều tiết xả lũ an toàn cho các hồ chứa thủy điện/thủy lợi và chủ động sơ tán dân cư trước các tình huống khẩn cấp.
  • Học viên cao học, nghiên cứu sinh và sinh viên ngành Khí tượng, Thủy văn, Biến đổi khí hậu: Đóng vai trò là tài liệu tham khảo chuẩn mực về phương pháp luận thẩm tra mô hình số trị (Verification Metrics), quy trình xử lý dữ liệu lưới NetCDF và kỹ thuật nội suy thống kê trong nghiên cứu khoa học khí quyển.

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

Mô hình nào cho kết quả dự báo mưa lớn tổng thể tốt nhất cho miền Trung và Tây Nguyên? Mô hình HRM-gfs đạt kỹ năng toàn diện nhất với sai số trung bình ME thấp nhất, hệ số tương quan r duy trì ổn định ở mức 0.40 - 0.45 ở hạn dự báo 24 giờ trên chuỗi 64 trạm, đồng thời kiểm soát diện mưa vừa và mưa rất to chuẩn xác nhất.

Tại sao mô hình WRFARW-gfs thường dự báo diện mưa lớn hơn thực tế đo đạc? Nguyên nhân xuất phát từ sơ đồ tham số hóa đối lưu Kain-Fritsch kích hoạt độ bất ổn định khí quyển mạnh mẽ, khiến mô hình có xu hướng phóng đại diện tích vùng mưa vừa (16-50 mm/24h) với chỉ số BIAS trung bình đạt xấp xỉ 1.4 trên toàn khu vực.

Kỹ năng dự báo của các mô hình biến đổi như thế nào khi kéo dài hạn dự báo từ 24h lên 72h? Kỹ năng dự báo suy giảm tuyến tính rõ rệt: ở hạn 72 giờ, sai số bình phương trung bình quân phương RMSE tăng khoảng 25-35%, hệ số tương quan r giảm xuống dưới mức 0.30 và xác suất phát hiện đúng hiện tượng POD giảm từ 15% đến 25% so với hạn 24 giờ.

Tại sao cần kết hợp cả chỉ số sai số liên tục và chỉ số phân hạng theo ngưỡng khi đánh giá mưa? Do trường mưa mang tính phi liên tục cao, các chỉ số liên tục (ME, RMSE) chỉ đo lường độ lớn sai số lượng nước (mm), trong khi bảng liên hợp (POD, FAR, ETS) giúp thẩm tra chuẩn xác xác suất bắt trúng sự kiện cực đoan mà không bị sai lệch bởi các ngày không mưa.

Kết quả nghiên cứu của luận văn có thể ứng dụng trực tiếp cho vận hành hồ chứa thủy điện không? Hoàn toàn khả thi. Dữ liệu dự báo định lượng 24h-72h được hiệu chỉnh sai số theo từng lưu vực trạm giúp các kỹ sư hồ chứa dự tính chính xác lưu lượng nước về hồ trước 1 đến 2 ngày, giảm thiểu 20-30% nguy cơ xả lũ đột ngột gây ngập lụt hạ du.

Kết luận

  • Thẩm tra toàn diện và có hệ thống kỹ năng dự báo mưa lớn của 3 mô hình số trị khu vực (HRM-gfs, WRFARW-gfs, WRFNMM-gfs) trên chuỗi số liệu thực nghiệm 3 năm (2008 - 2010) tại mạng lưới 64 trạm khí tượng miền Trung và Tây Nguyên.
  • Khẳng định HRM-gfs là công cụ có độ tin cậy và tính ổn định cao nhất toàn vùng (r ~ 0.45, ME tối ưu), trong khi WRFARW-gfs chiếm ưu thế ở ngưỡng mưa to 50-100 mm/24h và WRFNMM-gfs luôn dự báo thiên thấp.
  • Phân tích chi tiết quy luật phân hóa sai số không gian trên 4 tiểu vùng địa hình, chỉ ra Nam Trung Bộ đạt tương quan dự báo cao nhất (r ~ 0.50) trong khi Bắc Trung Bộ và Trung Trung Bộ chịu ảnh hưởng mạnh bởi hiệu ứng chắn gió địa hình.
  • Làm sáng tỏ quy luật suy giảm chất lượng dự báo từ hạn 24h đến 72h, cung cấp bằng chứng định lượng phục vụ việc lựa chọn hạn dự báo an toàn trong cảnh báo thiên tai.
  • Đóng góp bộ cơ sở khoa học thực chứng vững chắc phục vụ công tác chuẩn hóa nghiệp vụ dự báo khí tượng và định hướng phát triển công nghệ đồng hóa dữ liệu quy mô vừa tại Việt Nam.

Lộ trình triển khai tiếp theo cần tập trung ứng dụng công nghệ đồng hóa dữ liệu radar 3DVAR và kỹ thuật học máy khử sai số MOS trong giai đoạn 12 - 24 tháng tới. Các cơ quan quản lý, trung tâm khí tượng và viện nghiên cứu chuyên ngành hãy chủ động khai thác toàn văn luận văn để tích hợp bộ trọng số mô hình tối ưu vào quy trình tác nghiệp thực tế, góp phần nâng cao năng lực chủ động ứng phó thiên tai và bảo vệ an toàn cộng đồng.