Tổng quan nghiên cứu

Khu vực Đồng bằng Bắc Bộ (ĐBBB) có lượng mưa trung bình năm dao động từ 1.600 mm đến 1.800 mm với khoảng 130 đến 140 ngày mưa mỗi năm. Mùa mưa kéo dài từ tháng 5 đến tháng 10, tập trung tới 85% tổng lượng giáng thủy toàn năm, trong đó lượng mưa tích lũy của ba tháng cao điểm (tháng 7, 8, 9) chiếm trên 49%. Các hiện tượng mưa lớn cực đoan thường gây ra ngập lụt nghiêm trọng tại các đô thị lớn như Hà Nội, Nam Định và lũ quét tại các khu vực bán sơn địa như Ba Vì, Ninh Bình, điển hình như trận mưa lịch sử năm 2008 ghi nhận lượng mưa ngày tại trạm Hà Đông lên tới 514 mm.

Dự báo mưa lớn định lượng luôn là một trong những thách thức lớn nhất của ngành khí tượng thủy văn do tính bất định cao của các quá trình đối lưu quy mô vừa. Hệ thống dự báo tổ hợp hạn ngắn (SREPS) được vận hành nghiệp vụ tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cung cấp nguồn thông tin đa dạng nhưng các dự báo viên còn gặp nhiều khó khăn do thiếu các đánh giá định lượng chi tiết về sai số theo không gian, thời gian và các ngưỡng mưa.

Mục tiêu cốt lõi của nghiên cứu là đánh giá có hệ thống kỹ năng dự báo mưa lớn của 20 thành phần dự báo thuộc hệ thống SREPS và thử nghiệm 3 phương án tổ hợp tối ưu hóa kết quả dự báo. Nghiên cứu được triển khai trên tập số liệu gồm 155 đợt mưa lớn diện rộng trong 8 năm (2010–2017) tại 14 trạm quan trắc khí tượng bề mặt trên toàn khu vực ĐBBB cho các hạn dự báo 24 giờ, 48 giờ và 72 giờ. Kết quả đánh giá cung cấp cơ sở khoa học định lượng giúp các dự báo viên lựa chọn sản phẩm mô hình phù hợp, góp phần nâng cao độ tin cậy của các bản tin cảnh báo thiên tai và giảm thiểu rủi ro kinh tế - xã 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 lý thuyết dự báo thời tiết bằng mô hình số trị (NWP) và phương pháp dự báo tổ hợp đa mô hình - đa phân tích. Hệ thống SREPS tích hợp 4 mô hình động lực khu vực quy mô vừa (HRM, WRF-ARW, WRF-NMM, BoLAM) kết hợp với điều kiện biên và trường ban đầu từ 5 mô hình toàn cầu lớn (GFS của Mỹ, GME của Đức, NOGAPS của Hải quân Mỹ, GEM của Canada, GSM của Nhật Bản) để tạo ra 20 thành phần dự báo độc lập.

Các khái niệm chuyên ngành cốt lõi được áp dụng xuyên suốt bao gồm:

  • Dự báo giáng thủy định lượng (QPF): Ước lượng giá trị mưa tích lũy trong khoảng thời gian xác định (R24 từ 19 giờ hôm trước đến 19 giờ hôm sau).
  • Dự báo tất định từ tập tổ hợp: Kỹ thuật tổng hợp từ 20 trường dự báo thành phần để đưa ra một trường dự báo đại diện tối ưu duy nhất.
  • Bảng phân loại nhị phân 2x2: Cấu trúc ma trận đánh giá hiện tượng xuất hiện (Hit, Miss, False Alarm, Correct Negative) cho các ngưỡng mưa định danh.

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

Tập dữ liệu nghiên cứu bao gồm 155 ngày mưa lớn diện rộng được lựa chọn nghiêm ngặt theo tiêu chí nghiệp vụ của Trung tâm Dự báo Khí tượng Thủy văn Quốc gia (ngày có từ 50% số trạm trở lên trong khu vực đạt lượng mưa R24 từ 16 mm trở lên) trong giai đoạn từ năm 2010 đến 2017. Phân bố các ngày mưa lớn theo từng năm gồm: 2010 (13 ngày), 2011 (21 ngày), 2012 (23 ngày), 2013 (28 ngày), 2014 (15 ngày), 2015 (13 ngày), 2016 (17 ngày) và 2017 (25 ngày). Dữ liệu thực đo được thu thập đồng bộ tại 14 trạm khí tượng bề mặt đại diện cho toàn bộ địa hình ĐBBB.

Dữ liệu dự báo số trị được trích xuất từ phiên chạy 12 UTC (19 giờ Hà Nội) của hệ thống SREPS trên lưới tính toán 0.15° (khoảng 15 km) gồm 162 x 135 nút lưới, bao phủ tọa độ 99.9°E–124.05°E và 9.9°N–25.0°N. Nghiên cứu áp dụng phương pháp nội suy điểm gần nhất kết hợp mặt nạ đất liền - biển để gán giá trị dự báo từ nút lưới về điểm trạm quan trắc, nhằm loại bỏ hiện tượng sai lệch do lấy nhầm các điểm lưới trên mặt biển.

Để đánh giá toàn diện, hai nhóm chỉ số thống kê toán học được áp dụng:

  1. Nhóm chỉ số định lượng cho biến liên tục: 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 tuyến tính Pearson (r). Lý do chọn nhóm chỉ số này là nhằm xác định chính xác độ lệch hệ thống và biên độ dao động sai số của mô hình.
  2. Nhóm chỉ số dự báo pha cho biến nhị phân: Tỷ số thiên lệch tần suất (FBI), xác suất phát hiện (POD) và tỷ lệ báo động sai (FAR) phân theo 3 cấp độ mưa: mưa vừa (16–50 mm/24h), mưa to (51–100 mm/24h) và mưa rất to (> 100 mm/24h).

Ba phương án tính toán dự báo tổ hợp tuyến tính được xây dựng và so sánh:

  • Phương án 1: Phương pháp trung bình số học đơn giản với trọng số đều $w_i = 1/20$.
  • Phương án 2: Phương pháp trọng số nghịch đảo theo sai số RMSE của từng mô hình thành phần.
  • Phương án 3: Phương pháp gán trọng số tối ưu dựa trên mô hình hồi quy tuyến tính đa biến được ước lượng trên ngôn ngữ lập trình thống kê R.

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

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

Thứ nhất, đối với hạn dự báo 24 giờ, toàn bộ 20 mô hình thành phần của SREPS đều thể hiện xu hướng dự báo thấp hơn thực tế với ME âm dao động từ -5 mm đến -15 mm. Sai số tuyệt đối trung bình MAE nằm trong khoảng 25 mm đến 27 mm và RMSE dao động từ 33 mm đến 36 mm. Hệ số tương quan Pearson đạt mức 0.2 đến 0.5. Thành phần M05 (mô hình HRM chạy với điều kiện biên GSM) cho kết quả dự báo định lượng tốt nhất với RMSE thấp nhất và hệ số tương quan cao nhất. Nhóm mô hình BoLAM (thành phần M16 đến M20) duy trì sai số ổn định nhất trên toàn bộ 14 trạm quan trắc.

Thứ hai, khi kéo dài hạn dự báo lên 48 giờ và 72 giờ, chất lượng dự báo định lượng suy giảm rõ rệt. Giá trị ME chuyển từ âm sang dương (ME từ +4 mm đến +8 mm ở hạn 48 giờ và từ +10 mm đến +15 mm ở hạn 72 giờ), cho thấy mô hình chuyển sang xu hướng dự báo thiên cao. Hệ số tương quan ở hạn 48 giờ giảm xuống dưới 0.1 và phần lớn các thành phần đều cho hệ số tương quan âm ở hạn dự báo 72 giờ.

Thứ ba, kỹ năng dự báo pha phụ thuộc lớn vào ngưỡng mưa. Ở cấp mưa vừa (16–50 mm/24h), chỉ số FBI duy trì quanh giá trị 1.0 đến 1.2, xác suất phát hiện POD đạt từ 40% đến 60%, tỷ lệ báo động sai FAR dao động từ 35% đến 50%. Tuy nhiên, ở các cấp mưa to (51–100 mm/24h) và mưa rất to (> 100 mm/24h), chỉ số FBI giảm sâu dưới 1.0, POD giảm mạnh xuống chỉ còn 0% đến 20%, trong khi FAR tăng vọt lên 80% đến 100% ở hạn dự báo 72 giờ.

Thứ tư, so sánh 3 phương án tổ hợp cho thấy Phương án 3 (hồi quy đa biến) đạt hiệu quả dự báo định lượng vượt trội nhất với MAE và RMSE thấp nhất, đồng thời duy trì hệ số tương quan cao nhất ở cả 3 hạn dự báo. Phương án 1 (trung bình đơn giản) cho kết quả ổn định thứ hai. Phương án 2 (trọng số nghịch đảo RMSE) kiểm soát tốt chỉ số FBI ở ngưỡng mưa vừa nhưng lại bộc lộ hạn chế lớn khi không bắt được các hiện tượng mưa cực đoan trên 100 mm.

Thảo luận kết quả

Kết quả sai số định lượng và dự báo pha được trực quan hóa qua chuỗi biểu đồ phân tán, biểu đồ hộp sai số ME/RMSE theo các hạn dự báo và bảng phân loại nhị phân tổng hợp. Nguyên nhân chính khiến hệ thống SREPS dự báo thiếu hụt lượng mưa ở hạn 24 giờ là do độ tán của tập tổ hợp còn tương đối hẹp, khiến mô hình chưa nắm bắt đầy đủ các nhiễu động đối lưu cục bộ có cường độ mạnh. Hiện tượng mô hình dự báo thiên cao ở hạn 48 giờ và 72 giờ phản ánh xu thế khuếch tán diện rộng của các sơ đồ tham số hóa đối lưu (như Tiedtke, Kain-Fritsch, Grell-Devenyi) khi tích phân phương trình theo thời gian dài.

So sánh với các nghiên cứu quốc tế, phát hiện này hoàn toàn tương đồng với kết quả của Trung tâm Khí tượng Anh (Met Office) trên hệ thống MOGREPS và nghiên cứu của McBride & Ebert tại Úc, khi các mô hình số trị quy mô khu vực đều gặp khó khăn lớn trong việc dự báo chính xác vị trí và độ lớn của tâm mưa cực trị vượt ngưỡng 50 mm/24h. Sự vượt trội của phương án hồi quy đa biến khẳng định rằng việc tối ưu hóa trọng số dựa trên chuỗi số liệu nền lịch sử có khả năng hiệu chỉnh sai số hệ thống tốt hơn nhiều so với việc chỉ gán trọng số cơ học theo sai số RMSE đơn lẻ.

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

Dựa trên các kết luận rút ra từ thực nghiệm, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm tối ưu hóa công tác dự báo mưa lớn:

  1. Ứng dụng thuật toán tổ hợp hồi quy tuyến tính đa biến (Phương án 3): Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần tích hợp module hồi quy đa biến vào quy trình vận hành tự động của SREPS, đặt mục tiêu giảm sai số RMSE xuống dưới 28 mm cho hạn dự báo 24 giờ và hoàn thành chuyển giao nghiệp vụ trong vòng 6 đến 12 tháng.
  2. Tối ưu hóa việc lựa chọn mô hình thành phần: Đội ngũ dự báo viên tại Đài Khí tượng Thủy văn khu vực Đồng bằng Bắc Bộ cần ưu tiên tham khảo các sản phẩm từ mô hình BoLAM và mô hình HRM kết hợp điều kiện biên GSM trong các tình huống thời tiết nguy hiểm, nhằm nâng xác suất phát hiện (POD) các đợt mưa diện rộng lên trên 65% trong giai đoạn 2024–2025.
  3. Nâng cấp độ phân giải mô hình và đồng hóa dữ liệu: Tổng cục Khí tượng Thủy văn cần đầu tư nâng cấp độ phân giải ngang của hệ thống từ 0.15° xuống dưới 0.05° (khoảng 3–5 km), tích hợp đồng hóa dữ liệu viễn thám từ mạng lưới radar thời tiết và vệ tinh địa tĩnh, hướng tới mục tiêu kéo giảm tỷ lệ báo sai (FAR) ở ngưỡng mưa to xuống dưới 40% trước năm 2026.
  4. Xây dựng bảng tra cứu trọng số động theo hình thế thời tiết: Các nhóm nghiên cứu chuyên trách cần phân loại trọng số mô hình riêng biệt cho từng dạng hình thế điển hình (rãnh thấp, xoáy thuận nhiệt đới, không khí lạnh kết hợp gió đông) để đưa chỉ số thiên lệch tần suất FBI tiệm cận mức lý tưởng 1.0 trong lộ trình 2024–2026.

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

  1. Dự báo viên khí tượng thủy văn tại các đài trung ương và địa phương: Nắm vững đặc trưng sai số của từng mô hình thành phần và áp dụng trực tiếp các phương án tổ hợp trọng số để hiệu chỉnh bản tin dự báo mưa lớn hàng ngày.
  2. Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Khí tượng học: Sử dụng tài liệu như một cẩm nang phương pháp luận hoàn chỉnh về quy trình hậu xử lý số liệu mô hình số trị và kỹ thuật kiểm định thống kê giáng thủy định lượng.
  3. Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các tỉnh ĐBBB: Hiểu rõ giới hạn tin cậy của các hạn dự báo 24 giờ, 48 giờ và 72 giờ để chủ động xây dựng phương án tiêu thoát nước đô thị và cảnh báo sơ tán dân cư sát với thực tế.
  4. Kỹ sư phát triển phần mềm và công nghệ môi trường: Khai thác các thuật toán nội suy điểm trạm, mô hình tổ hợp và cấu trúc xử lý dữ liệu khí tượng để tích hợp vào các hệ thống cảnh báo sớm thông minh phục vụ nông nghiệp và đô thị hóa.

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

Tại sao hệ thống SREPS lại có xu hướng dự báo lượng mưa thấp hơn thực tế ở hạn 24 giờ?

Ở hạn dự báo 24 giờ, toàn bộ 20 mô hình thành phần đều có sai số ME âm từ -5 mm đến -15 mm. Nguyên nhân là do độ phân giải 0.15° của mô hình chưa bao quát hết các tâm mưa đối lưu quy mô nhỏ, đồng thời độ tán tổ hợp của hệ thống còn hẹp, dẫn đến việc ước lượng thấp cường độ của các tâm mưa cực đoan có lượng mưa trên 50 mm đến 100 mm/24h.

Trong 20 thành phần dự báo của SREPS, mô hình nào cho kết quả định lượng tốt nhất?

Mô hình HRM kết hợp điều kiện biên GSM (thành phần M05) và dòng mô hình BoLAM kết hợp các trường biên toàn cầu GME, GSM, NOGAPS (từ M17 đến M20) cho kết quả tốt nhất. Nhóm này đạt sai số RMSE thấp nhất (khoảng 33 mm) và duy trì hệ số tương quan r cao nhất (từ 0.35 đến 0.45), thể hiện độ ổn định cao qua 155 đợt mưa lớn.

Phương pháp tổ hợp nào mang lại độ chính xác cao nhất cho khu vực Đồng bằng Bắc Bộ?

Phương pháp tổ hợp tính trọng số bằng hồi quy tuyến tính đa biến (Phương án 3) mang lại hiệu quả cao nhất. Phương pháp này đạt giá trị ME nhỏ nhất, hạ thấp MAE và RMSE, đồng thời duy trì hệ số tương quan tuyến tính cao hơn đáng kể so với phương pháp trung bình số học đơn giản và phương pháp nghịch đảo RMSE.

Khả năng dự báo mưa to và rất to của mô hình thay đổi ra sao khi thời hạn dự báo tăng lên 72 giờ?

Kỹ năng dự báo suy giảm rất nhanh theo thời gian. Ở hạn 72 giờ với ngưỡng mưa to (> 50 mm) và rất to (> 100 mm), xác suất phát hiện POD giảm sâu xuống dưới 10%, trong khi tỷ lệ báo động sai FAR tăng lên mức 90% đến 100%, đồng thời hệ số tương quan giữa dự báo và quan trắc tiệm cận mức 0 hoặc mang giá trị âm.

Thuật toán nào được sử dụng để chuyển đổi dữ liệu từ lưới mô hình về trạm quan trắc?

Luận văn áp dụng thuật toán nội suy điểm gần nhất kết hợp bộ lọc mặt nạ đất liền - biển. Phương pháp này tính toán khoảng cách hình học từ trạm quan trắc đến các nút lưới xung quanh và chọn nút lưới đất liền gần nhất, giúp bảo toàn tính chất giáng thủy cục bộ và tránh sai số do ảnh hưởng của các điểm nút nằm ngoài biển.

Kết luận

  • Luận văn đã thực hiện đánh giá toàn diện và hệ thống kỹ năng dự báo mưa lớn của 20 thành phần thuộc hệ thống SREPS và 3 phương pháp tổ hợp qua 155 đợt mưa diện rộng giai đoạn 2010–2017 trên 14 trạm khí tượng ĐBBB.
  • Xác định rõ quy luật sai số: Mô hình có xu hướng dự báo thiên thấp ở hạn 24 giờ (ME âm) nhưng chuyển dịch sang dự báo thiên cao ở hạn 48 giờ và 72 giờ (ME dương).
  • Khẳng định tính ổn định vượt trội của mô hình BoLAM và chứng minh phương pháp tổ hợp hồi quy tuyến tính đa biến mang lại chất lượng dự báo định lượng tối ưu nhất.
  • Chỉ rõ giới hạn của hệ thống đối với các ngưỡng mưa cực đoan (> 50 mm và > 100 mm/24h), khi chỉ số POD suy giảm mạnh và FAR tăng cao ở các hạn dự báo xa.
  • Đóng góp bộ cơ sở dữ liệu kiểm định khoa học và giải pháp thực tiễn có giá trị ứng dụng cao cho ngành khí tượng thủy văn Việt Nam.

Kế hoạch tiếp theo cần tập trung chuẩn hóa thuật toán hồi quy vào hệ thống dự báo nghiệp vụ trong năm 2024, đồng thời mở rộng nghiên cứu thử nghiệm trên chuỗi số liệu cập nhật giai đoạn 2018–2025. Các chuyên gia, dự báo viên và nhà nghiên cứu quan tâm được khuyến khích khai thác toàn văn tài liệu để ứng dụng hiệu quả vào công tác cảnh báo thiên tai và nghiên cứu chuyên sâu.