Tổng quan nghiên cứu

Dự báo định lượng mưa lớn luôn là một trong những bài toán phức tạp và thách thức nhất trong nghiệp vụ khí tượng thủy văn tại Việt Nam. Vùng duyên hải miền Trung với địa hình hẹp, dốc, một bên là dãy Trường Sơn hùng vĩ, một bên là Biển Đông, 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. Điển hình như đợt mưa lịch sử từ ngày 13 đến 17/10/2010, tổng lượng mưa phổ biến đạt từ 500 đến 800 mm, cá biệt tại trạm Hòa Duyệt (Hà Tĩnh) đạt tới 903 mm, gây ngập lụt diện rộng tại các tỉnh từ Nghệ An đến Quảng Bình.

Vấn đề cốt lõi của các mô hình dự báo số trị hiện nay là sự nhạy cảm cao với điều kiện ban đầu; những sai số nhỏ ở trường khởi tạo có thể dẫn đến sự sai lệch lớn trong dự báo sau 24 đến 72 giờ. Mục tiêu cụ thể của luận văn là nghiên cứu, thiết kế và thử nghiệm hệ thống đồng hóa số liệu lọc Kalman tổ hợp địa phương (LETKF) tích hợp vào mô hình thời tiết quy mô vừa phi thủy tĩnh WRF (phiên bản ARW V3.2). Nghiên cứu tập trung đánh giá khả năng mô phỏng diện mưa và lượng mưa tích lũy cho đợt mưa lớn tại miền Trung trong hạn dự báo 72 giờ.

Ý nghĩa của công trình thể hiện ở việc nâng cao độ chính xác dự báo định lượng mưa, giảm sai số vị trí tâm mưa lớn từ 50 đến 100 km so với mô hình tham chiếu không đồng hóa. Đồng thời, nghiên cứu mở ra hướng ứng dụng công nghệ đồng hóa tổ hợp tiên tiến nhằm hiện đại hóa công tác cảnh báo thiên tai tại Việt Nam.

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 ứng dụng lý thuyết đồng hóa số liệu chuỗi và phương pháp lọc Kalman tổ hợp địa phương (LETKF - Local Ensemble Transform Kalman Filter), kết hợp với hệ thống phương pháp số trị của mô hình dự báo khí tượng WRF-ARW. Về mặt toán học và động lực học khí quyển, hệ thống vận hành dựa trên 4 khái niệm và cấu trúc lý thuyết nền tảng:

  • Ma trận sai số hiệp biến nền ($P^f$ hoặc $B$): Đại diện cho độ bất định của trường dự báo tại thời điểm phân tích, được tính toán động lực từ một tập hợp các thành phần tổ hợp thay vì giả thiết cố định như trong phương pháp nội suy tối ưu (OI) hay đồng hóa biến phân 3 chiều (3DVAR).
  • Toán tử quan trắc ($H$): Thực hiện phép chiếu không gian từ các biến trạng thái của mô hình sang không gian quan trắc thực tế, cho phép tính toán độ lệch gia số quan trắc (innovation vector).
  • Không gian con tổ hợp địa phương: LETKF thực hiện tính toán biến đổi độc lập tại từng điểm nút lưới thông qua vùng lân cận được xác định trước, giúp chuyển bài toán nghịch đảo ma trận kích thước lớn sang không gian số chiều nhỏ tương ứng với số thành phần tổ hợp ($K$).
  • Hệ phương trình thủy nhiệt động phi thủy tĩnh WRF: Mô phỏng sự tiến triển của các trường khí tượng như độ cao địa thế vị, nhiệt độ, độ ẩm tương đối, tốc độ và hướng gió theo không gian 3 chiều.

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

Nguồn dữ liệu đầu vào và điều kiện biên của mô hình được trích xuất từ tập số liệu phân tích toàn cầu GFS (NCEP, Hoa Kỳ) với độ phân giải không gian 1.0° x 1.0°. Nguồn số liệu đồng hóa bao gồm mạng lưới trạm thám không vô tuyến cao không khu vực châu Á và tập số liệu thám không độ phân giải cao tại 5 trạm nghiệp vụ trọng điểm của Việt Nam gồm Hà Nội, Điện Biên, Vinh, Đà Nẵng và Tân Sơn Hòa, với độ cao quan trắc trung bình đạt trên 28 km và tần suất quan trắc 2 lần/ngày (00Z và 12Z).

Mô hình WRF được thiết lập với cấu hình 2 miền lưới lồng tương tác hai chiều: miền ngoài có kích thước 91 x 91 điểm lưới với bước lưới 54 km (tâm tại 16°N; 108°E), miền trong bao phủ trọn vẹn khu vực miền Trung với 155 x 155 điểm lưới và bước lưới 18 km. Cỡ mẫu thử nghiệm bao gồm 7 chu trình dự báo liên tiếp (mỗi chu trình cách nhau 12 giờ) và quy mô tổ hợp gồm 5 thành phần ($K = 5$), tích phân thời gian 72 giờ từ 00Z ngày 14/10/2010 đến 00Z ngày 17/10/2010 với 25 bước thời gian trích xuất lượng mưa. Phương pháp chọn mẫu tập trung vào đợt thiên tai điển hình có sự tương tác mạnh mẽ giữa không khí lạnh tăng cường và dải hội tụ nhiệt đới. Lý do lựa chọn thuật toán LETKF vì phương pháp này cho phép xử lý tính toán song song độc lập tại từng điểm lưới, tối ưu hóa bộ nhớ và khắc phục triệt để giả định sai số tĩnh của các phương pháp đồng hóa cổ điển.

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

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

Quá trình thử nghiệm đối chứng giữa mô hình WRF chuẩn không đồng hóa (CTL) và hệ thống WRF tích hợp đồng hóa lọc Kalman (WRF-KF) đã mang lại các kết quả định lượng rõ rệt:

  1. Cải thiện đáng kể cường độ và phân bố không gian lượng mưa: Tại các thời điểm dự báo 12Z ngày 14/10 và 12Z ngày 15/10/2010, trường mưa dự báo từ mô hình WRF-KF tái hiện rất sát với thực tế quan trắc trạm. Mô hình WRF chuẩn (CTL) dự báo lượng mưa thấp hơn thực tế từ 40% đến 60%, trong khi WRF-KF đã nắm bắt được các tâm mưa lớn tại dải ven biển Hà Tĩnh và Quảng Bình.
  2. Khả năng phân loại và bắt các vùng không mưa: Mô hình WRF-KF thể hiện tính ưu việt khi phân định chính xác ranh giới các trạm có lượng mưa nhỏ hoặc không mưa ở sườn phía tây và phía bắc (như Quế Phong, Cửa Rào, Đô Lương với lượng mưa quan trắc dưới 25 mm/ngày), hạn chế tình trạng báo động giả thường gặp ở mô hình WRF đơn lẻ.
  3. Hiệu quả vượt trội của tập số liệu cao không nội địa: Khi đồng hóa bổ sung bộ số liệu thám không cao không chi tiết của 5 trạm Việt Nam, chất lượng dự báo lượng mưa tích lũy 24 giờ tại các trạm trọng điểm tăng thêm khoảng 25% so với việc chỉ sử dụng số liệu thám không khu vực châu Á thông thường.
  4. Giới hạn trong việc mô phỏng các đỉnh mưa cực đoan: Dù WRF-KF nâng mức dự báo lượng mưa ngày lên mức 150 đến 200 mm/ngày tại các tâm mưa lớn (Hương Khê, Linh Cảm, Chu Lễ), hệ thống vẫn chưa đạt tới các mức đo cực đoan thực tế (vốn dao động từ 456 đến 548 mm/ngày trong ngày 16/10/2010).

Thảo luận kết quả

Nguyên nhân chính giúp mô hình WRF-KF vượt trội hơn WRF-CTL là nhờ cơ chế cập nhật ma trận sai số hiệp biến phân tích liên tục qua từng chu trình 12 giờ. Quá trình này giúp trường ban đầu phản ánh chính xác các nhiễu động ẩm tầng thấp và sự tương tác giữa đới gió đông bắc lạnh với địa hình đón gió miền Trung. Dữ liệu đối chứng có thể được trực quan hóa thông qua bảng so sánh chỉ số đánh giá phân loại mưa (chỉ số đe dọa TS, chỉ số độ lệch BIAS) và biểu đồ phân tán giữa lượng mưa dự báo tổ hợp trung bình với chuỗi số liệu thực đo tại 12 trạm khí tượng thủy văn miền Trung.

Tuy nhiên, việc mô hình chưa nắm bắt trọn vẹn đỉnh mưa trên 400 mm/ngày xuất phát từ 2 yếu tố: bước lưới 18 km chưa đủ mịn để giải đoán các ổ đối lưu sâu quy mô cực nhỏ, và số lượng thành phần tổ hợp ($K = 5$) còn hạn chế do tài nguyên tính toán. So sánh với các nghiên cứu sử dụng 3DVAR trước đây của các nhà khoa học trong nước, phương pháp LETKF cho thấy sự ổn định hơn trong việc duy trì cấu trúc động lực của các hình thế gây mưa lớn.

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

Nhằm hoàn thiện và đưa công nghệ đồng hóa lọc Kalman tổ hợp vào ứng dụng nghiệp vụ thực tế, luận văn đưa ra 4 khuyến nghị then chốt:

  • Nâng cao độ phân giải miền tính: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần thiết lập lưới lồng chi tiết thứ ba với độ phân giải ngang từ 2 đến 3 km cho khu vực trọng điểm miền Trung trong giai đoạn 2024–2026, giúp mô hình tự giải đoán trực tiếp các quá trình đối lưu mây tích mà không cần sơ đồ tham số hóa thô.
  • Mở rộng quy mô thành phần tổ hợp: Tăng số lượng thành phần tổ hợp từ 5 lên 30 đến 50 thành phần trong vòng 12 đến 18 tháng tới. Mục tiêu đạt độ tin cậy phương sai sai số nền trên 85%, loại bỏ các tương quan giả trên quy mô không gian lớn.
  • Tích hợp đồng hóa đa nguồn dữ liệu thời gian thực: Các đơn vị nghiên cứu và dự báo cần phát triển mô-đun đồng hóa bổ sung dữ liệu phản xạ Radar Doppler ven biển (Đồng Hới, Tam Kỳ, Vinh) và dữ liệu gió vệ tinh độ phân giải cao, hướng tới chu kỳ cập nhật dữ liệu từ 1 đến 3 giờ/lần.
  • Tự động hóa và song song hóa hệ thống: Cục Khí tượng Thủy văn cần đầu tư nâng cấp hạ tầng tính toán hiệu năng cao (HPC), tối ưu hóa mã nguồn LETKF để quy trình từ tiếp nhận dữ liệu quan trắc đến xuất bản tin dự báo hoàn tất trong thời gian dưới 45 phút cho mỗi ốp chạy.

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

Tài liệu là công trình nghiên cứu chuyên sâu, mang lại giá trị thiết thực cho nhiều nhóm chuyên môn:

  • Dự báo viên khí tượng thủy văn: Nắm bắt quy trình thiết lập, vận hành hệ thống WRF-LETKF và hiểu rõ cơ chế lan truyền thông tin quan trắc để hiệu chỉnh các bản tin cảnh báo mưa lũ nguy hiểm.
  • Học viên cao học và nghiên cứu sinh ngành Khoa học Khí quyển: Sử dụng công trình như một tài liệu tham khảo chuẩn mực về toán học đồng hóa số liệu, các công thức dẫn xuất ma trận hiệp biến và phương pháp tích phân ngẫu nhiên Monte Carlo.
  • Kỹ sư phát triển mô hình số trị: Khai thác quy trình tự động hóa các khâu tiền xử lý (WPS), tích phân động lực (WRFV3), ghép nối module WRFDA và xử lý đồ họa hậu kỳ bằng GrADS.
  • Cán bộ quản lý rủi ro thiên tai: Tiếp cận các đánh giá khoa học về độ bất định của dự báo thời tiết, từ đó xây dựng các kịch bản sơ tán dân cư và vận hành hồ chứa thủy điện hợp lý hơn.

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

1. Phương pháp LETKF có điểm gì vượt trội hơn so với phương pháp 3DVAR truyền thống? Phương pháp LETKF tính toán ma trận sai số hiệp biến nền biến đổi theo thời gian và phụ thuộc vào dòng chảy thực tế thông qua các thành phần tổ hợp. Trong khi đó, 3DVAR sử dụng ma trận sai số tĩnh, không phản ánh được sự biến động nhanh của các nhiễu động đối lưu gây mưa lớn.

2. Tại sao cần đồng hóa riêng bộ số liệu thám không cao không của Việt Nam? Số liệu cao không quốc tế thường có độ phân giải thưa và độ trễ cao. Năm trạm thám không vô tuyến của Việt Nam cung cấp thông tin khí áp, nhiệt độ, độ ẩm và gió rất chi tiết theo phương thẳng đứng lên tới trên 28 km, giúp mô hình bắt chính xác hoàn lưu gió mùa và độ ẩm tầng thấp.

3. Vì sao mô hình chưa dự báo được các điểm mưa cực đoan trên 400 mm/ngày? Hạn chế này do bước lưới mô hình ở mức 18 km và số thành phần tổ hợp chỉ dừng lại ở 5 thành phần. Ở độ phân giải này, mô hình phải sử dụng sơ đồ tham số hóa đối lưu, dẫn đến hiện tượng làm mịn trường mưa và giảm đỉnh mưa cục bộ.

4. Chu kỳ đồng hóa số liệu 12 giờ được thực hiện như thế nào trong mô hình? Tại mỗi mốc thời gian (00Z và 12Z), dự báo hạn rất ngắn 12 giờ từ chu trình trước được lấy làm trường nền. Sau đó, hệ thống kết hợp trường nền này với dữ liệu quan trắc mới qua bộ lọc LETKF để tạo ra trường phân tích tối ưu cho chu trình kế tiếp.

5. Cần cấu hình phần cứng như thế nào để vận hành hệ thống WRF-LETKF? Để chạy ổn định hệ thống WRF-LETKF với tổ hợp 20 đến 30 thành phần và 2 miền lưới lồng, hệ thống đòi hỏi cụm máy chủ tính toán hiệu năng cao (HPC) tối thiểu từ 64 đến 128 core CPU, bộ nhớ RAM trên 256 GB và hệ thống lưu trữ tốc độ cao.

Kết luận

  • Luận văn đã xây dựng và thử nghiệm thành công quy trình đồng hóa số liệu lọc Kalman tổ hợp địa phương (WRF-LETKF) phục vụ dự báo mưa lớn miền Trung.
  • Kết quả thử nghiệm trên đợt mưa lịch sử tháng 10/2010 chứng minh WRF-KF vượt trội hơn mô hình WRF thông thường về khả năng mô phỏng diện mưa và lượng mưa tích lũy.
  • Việc tích hợp chuỗi số liệu thám không độ phân giải cao tại 5 trạm cao không Việt Nam tạo ra bước chuyển biến rõ rệt trong việc định vị tâm mưa lớn.
  • Kế hoạch nghiên cứu tiếp theo cần tập trung tăng quy mô tổ hợp lên 30–50 thành phần và giảm bước lưới xuống dưới 3 km trong giai đoạn 2024–2026.
  • Độc giả và các đơn vị nghiên cứu quan tâm có thể khai thác toàn văn tài liệu để ứng dụng trực tiếp khung thuật toán vào bài toán nghiệp vụ dự báo thời tiết nguy hiểm.