Tổng quan về luận án

Bối cảnh khoa học của công trình nghiên cứu xuất phát từ thực tiễn biến đổi khí hậu và tính chất dị thường của các hình thế thời tiết cực đoan tại khu vực miền Trung - Tây Nguyên. Lưu vực sông Ba là một trong 9 hệ thống sông lớn nhất Việt Nam với diện tích tự nhiên $13.417\text{ km}^2$, chiều dài dòng chính $374\text{ km}$ (tính đến cửa Đà Rằng là $388\text{ km}$), mật độ lưới sông tương đối thưa $0,22\text{ km/km}^2$ và độ dốc bình quân lưu vực $10,9%$. Về mặt khí tượng - thủy văn, lưu vực chịu sự chi phối phức tạp của cả hai đới khí hậu sườn Đông và Tây dãy Trường Sơn, tạo ra sự tương tác đa quy mô giữa các nhiễu động nhiệt đới, không khí lạnh (KKL), dải hội tụ nhiệt đới (ITCZ) và xoáy thuận nhiệt đới (XTNĐ). Trên lưu vực đã hình thành hệ thống bậc thang thủy điện và thủy lợi lớn gồm hồ Ka Nak ($W_{pl} = 133,24 \cdot 10^6\text{ m}^3$), Ia M’Lá ($156,74 \cdot 10^6\text{ m}^3$), Sông Ba Hạ ($66 \cdot 10^6\text{ m}^3$), Ayun Hạ ($131,78 \cdot 10^6\text{ m}^3$), Sông Hinh ($36,86 \cdot 10^6\text{ m}^3$) và Krông H’Năng ($12,62 \cdot 10^6\text{ m}^3$). Theo số liệu thực tế, "tổng dung tích phòng lũ của các hồ là $260,5 \cdot 10^6\text{ m}^3$, rất nhỏ so với tổng lượng lũ 5 ngày lớn nhất đã xảy ra tại Củng Sơn là $2.507,3 \cdot 10^6\text{ m}^3$". Sự chênh lệch quy mô dung tích này dẫn đến nguy cơ rất cao về việc "xả sốc xuống hạ du (lưu lượng xả tăng đột ngột từ $0\text{ m}^3/\text{s}$ đến giá trị lưu lượng lớn nhất đến hồ), gây ngập lụt và thiệt hại về người và tài sản cho hạ du".

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định: Các phương pháp dự báo dòng chảy truyền thống và mô hình thủy văn hạn ngắn ($6 - 24\text{ giờ}$, tối đa $48\text{ giờ}$) không đủ thời gian dự kiến (lead time) để các cấp quản lý và chủ hồ chủ động dung tích cắt lũ theo "Quyết định 878/QĐ-TTg ngày 18 tháng 7 năm 2018 của Thủ tướng Chính phủ" (thay thế Quyết định 1077/QĐ-TTg). Hơn nữa, việc sử dụng trực tiếp các sản phẩm dự báo mưa số trị (NWP) toàn cầu độ phân giải thô ($14\text{ km} - 50\text{ km}$) hoặc phương pháp thống kê synop truyền thống không nắm bắt được quá trình đối lưu sâu quy mô vừa (mesoscale) trên địa hình chia cắt phức tạp. Luận án giải quyết ba câu hỏi nghiên cứu (Research Questions - RQ) và các giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Cấu hình mô hình khí tượng số trị khu vực bất thủy tĩnh nào có khả năng mô phỏng chi tiết trường mưa định lượng cực trị hạn vừa ($72 - 120\text{ giờ}$) trên lưu vực sông Ba?
    • H1: Mô hình WRF-ARW độ phân giải cao ($5\text{ km} \times 5\text{ km}$) lồng ghép điều kiện biên từ mô hình toàn cầu IFS ($14\text{ km}$) sẽ cải thiện đáng kể kỹ năng dự báo mưa hạn vừa so với mô hình toàn cầu đơn lẻ.
  • RQ2: Phương pháp hậu xử lý thống kê (MOS / Regression) nào tối ưu hóa được sai số hệ thống của mưa dự báo số trị trước khi nạp vào mô hình thủy văn?
    • H2: Hiệu chỉnh sai số trung bình (Mean Error - ME) theo phương trình hồi quy phân bố cho từng ô lưới sẽ nâng cao hệ số tương quan (HSTQ) và giảm sai số quân phương (RMSE) của trường mưa dự báo.
  • RQ3: Mức độ lan truyền sai số từ trường mưa dự báo hiệu chỉnh vào các mô hình dòng chảy (MIKE NAM và HEC-HMS) ảnh hưởng như thế nào đến độ chính xác dự báo lưu lượng đỉnh lũ và thời gian xuất hiện đỉnh tại 4 hồ chứa trọng điểm?
    • H3: Chuỗi tích hợp mô hình khí tượng phân giải cao - thủy văn bán phân bố cho phép kéo dài hạn dự báo dòng chảy đến hồ lên $3 - 5\text{ ngày}$ với mức bảo đảm kỹ thuật trên $70%$.

Khung lý thuyết của công trình thiết lập trên sự tích hợp giữa lý thuyết động lực khí quyển phi thủy tĩnh (Non-hydrostatic Atmospheric Dynamics), lý thuyết hình thành dòng chảy do mưa (Rainfall-Runoff Generation Mechanisms của Horton và Dunne) và lý thuyết vận hành tối ưu hệ thống liên hồ chứa. Phạm vi không gian tập trung vào 4 tiểu lưu vực hồ chứa lớn: An Khê ($1.200\text{ km}^2$), Ayun Hạ ($2.950\text{ km}^2$), Krông H’Năng ($1.761\text{ km}^2$) và Sông Hinh ($1.040\text{ km}^2$), với chuỗi dữ liệu thực nghiệm kiểm chứng qua các trận lũ lịch sử giai đoạn 1986 - 2016 (trọng tâm kiểm định chi tiết các năm 2013, 2014, 2015, 2016 và đợt lũ đặc biệt lớn năm 2017).

+------------------------------------------------------------------------------------------------+
|                        KHUNG TÍCH HỢP ĐỘNG LỰC DỰ BÁO MƯA - LŨ HẠN VỪA                         |
+------------------------------------------------------------------------------------------------+
|  Mô hình Toàn cầu IFS (ECMWF)       Mô hình Khu vực WRF-ARW         Hiệu chỉnh Thống kê MOS    |
|  [Độ phân giải 14km x 14km]   --->  [Độ phân giải 5km x 5km]   ---> [Giảm sai số ME, tăng HSTQ]|
+------------------------------------------------------------------------------------------------+
                                              |
                                              v
+------------------------------------------------------------------------------------------------+
|  Trường mưa ô lưới phân bố          Mô hình Thủy văn               Vận hành Liên hồ chứa       |
|  [Ma trận trọng số không gian] ---> [MIKE NAM / HEC-HMS]      ---> [Quyết định 878/QĐ-TTg]     |
|                                     [Hạn dự báo 72 - 120 giờ]      [Cắt lũ, xóa bỏ xả sốc]     |
+------------------------------------------------------------------------------------------------+

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về dự báo mưa và lũ hạn vừa phản ánh sự chuyển dịch mô thức từ phương pháp thống kê kinh nghiệm thuần túy sang dự báo tổ hợp và chuỗi mô hình động lực tích hợp:

           TIẾN TRÌNH PHÁT TRIỂN CÔNG NGHỆ DỰ BÁO KHÍ TƯỢNG - THỦY VĂN
[Thập niên 1990]  Hệ thống EPS toàn cầu (NCEP BGM: Toth & Kalnay 1997; ECMWF SV: Palmer 1992)
[Thập niên 2000]  Dự báo tổ hợp khu vực LEPS (COSMO-LEPS: Montani 2003; Marsigli 2005)
[Giai đoạn 2010+] Tích hợp Siêu tổ hợp (TIGGE: Bougeault 2009; NAEFS: Zhu 2009)
[Luận án (2020)]  IFS (14km) -> WRF-ARW (5km) -> MOS Ô lưới -> MIKE NAM / HEC-HMS (Lưu vực sông Ba)
  1. Trường phái Dự báo thời tiết số trị (NWP) và Hệ thống Dự báo Tổ hợp (EPS): Nền tảng dự báo hạn vừa toàn cầu được định hình bởi Toth và Kalnay (1997) với phương pháp gây nhiễu động mầm phát triển (Breeding of Growing Mode - BGM) tại NCEP; Palmer và cộng sự (1992) với phương pháp Véc-tơ kỳ dị (Singular Vector - SV) tại ECMWF; Houtekamer và cộng sự (1996) ứng dụng Ensemble Kalman Filter (EnKF) tại Trung tâm Khí tượng Canada (CMC). Nhằm chuyển tiếp từ quy mô toàn cầu xuống quy mô lưu vực, các nghiên cứu của Molteni và cộng sự (2001), Marsigli và cộng sự (2001, 2005) và Montani và cộng sự (2003) đã phát triển hệ thống tổ hợp khu vực hạn chế (COSMO-LEPS). Philippe Bougeault và cộng sự (2009) cùng Zhu và cộng sự (2009, 2016) đẩy mạnh tích hợp siêu tổ hợp TIGGE và NAEFS.
  2. Trường phái Mô hình hóa Thủy văn kết nối Mưa số trị định lượng (QPF): Gouweleeuw và cộng sự (2005) chứng minh mô hình LISFLOOD kết hợp EPS 51 thành phần của ECMWF có thể cảnh báo lũ sớm trước $2 - 3\text{ ngày}$ trên sông Meuse và sông Odra, nhưng chỉ ra hạn chế đỉnh lũ mô phỏng thường thấp hơn thực tế do độ phân giải mưa còn thô. Pappenberger và cộng sự (2008) tại Romania và Thiemig và cộng sự (2015) với hệ thống AFFS tại lưu vực sông Sabi (Châu Phi) khẳng định dự báo mưa ECMWF mang lại hiệu quả cao đối với lưu vực trên $10.000\text{ km}^2$. Knebl và cộng sự (2005) tại sông San Antonio (Texas, Hoa Kỳ) đã tích hợp radar NEXRAD với HEC-HMS và HEC-RAS qua công cụ Map-to-Map trên GIS.
  3. Trường phái Nghiên cứu Khí tượng Thủy văn tại Việt Nam: Nghiên cứu synop của Nguyễn Ngọc Thục và Lương Tuấn Minh (1990), Nguyễn Ngọc Thục và nnk (1994), Lê Bắc Huỳnh (2000), Lê Văn Ánh (2004) đã tổng kết các tổ hợp gây mưa lớn tại Trung Bộ. Nghiên cứu dự báo số trị nghiệp vụ của Bùi Minh Tăng (2014), Dư Đức Tiến và nnk (2011), Võ Văn Hòa và nnk (2017) đã bước đầu thử nghiệm mô hình GSM, GFS, NAEFS, IFS và HRM. Về mặt thủy văn, Nguyễn Lan Châu (2004) dùng Marine cho sông Đà; Đặng Thanh Mai (2005, 2011, 2014) phát triển WETSPA, HEC-RAS, MIKE cho sông Thu Bồn, sông Ba và các sông Nam Trung Bộ; Nguyễn Viết Thi (2006), Bùi Đình Lập (2014), Lương Hữu Dũng (2016), Nguyễn Việt (2018) nghiên cứu vận hành liên hồ chứa sông Ba bằng Fuzzy Logic và quy hoạch động.

Tranh luận học thuật cốt lõi tồn tại giữa hai quan điểm: (1) Quan điểm duy hình thức (Statistical/Empirical) cho rằng đối với mạng lưới trạm thưa thớt, việc hiệu chỉnh thống kê trực tiếp quan trắc trạm bề mặt cho kết quả đáng tin cậy hơn; đối lập với (2) Quan điểm động lực phi thủy tĩnh (Dynamical Non-hydrostatic Downscaling) khẳng định chỉ có mô hình giải phương trình khí quyển phân giải cao mới tái hiện được bản chất đối lưu địa hình sườn Đông - Tây Trường Sơn.

So sánh đối chuẩn với 02 nghiên cứu quốc tế tiêu biểu:

  • So sánh 1: Walter Collischonn và cộng sự (2007) sử dụng mô hình MGB-IPH kết hợp dự báo mưa ETA (độ phân giải thô $\approx 40\text{ km}$) cho lưu vực sông Paranaíba ($75.000\text{ km}^2$, Brazil). Nghiên cứu của Collischonn chứng minh việc giảm sai số dự báo dòng chảy phụ thuộc trực tiếp vào cải thiện mưa đầu vào. Luận án của Nguyễn Văn Hiếu phát triển xa hơn công trình này bằng cách hạ quy mô động lực xuống ô lưới $5\text{ km} \times 5\text{ km}$ bằng WRF-ARW lồng ghép biên IFS ($14\text{ km}$), sau đó áp dụng hồi quy thống kê từng ô lưới, giúp giải quyết thỏa đáng tính bất đồng nhất địa hình của lưu vực sông Ba vốn hẹp và dốc hơn nhiều so với sông Paranaíba.
  • So sánh 2: B.T. Gouweleeuw và cộng sự (2005) áp dụng ECMWF trực tiếp vào LISFLOOD trên sông Meuse và Odra bị hiện tượng đánh giá thấp đỉnh lũ (underestimation) do mưa toàn cầu bị san phẳng cực trị. Luận án khắc phục triệt để khiếm khuyết này thông qua việc thiết lập thuật toán hiệu chỉnh sai số trung bình (ME) trên từng ô lưới cho các cấp mưa to ($>50\text{ mm}/24\text{h}$) và rất to ($>100\text{ mm}/24\text{h}$), khôi phục chính xác phân bố cực trị của các đợt mưa bão kết hợp KKL.

Vị thế khoa học của luận án: Lần đầu tiên tại Việt Nam, một chuỗi công nghệ khép kín từ mô hình khí tượng toàn cầu phân giải cao nhất thế giới hiện nay (IFS $14\text{ km}$), chuyển tiếp động lực qua mô hình khu vực bất thủy tĩnh WRF ($5\text{ km}$), qua mô đun hậu xử lý thống kê ô lưới, nạp trực tiếp vào hai hệ thống mô hình thủy văn chuẩn mực (MIKE NAM và HEC-HMS) được thiết lập hoàn chỉnh, giải quyết trọn vẹn bài toán mở rộng thời gian dự báo lũ phục vụ quy trình liên hồ chứa 878/QĐ-TTg.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và làm sâu sắc thêm các lý thuyết chuyên ngành:

  1. Mở rộng Thuyết tương tác Động lực Khí quyển - Thủy văn (Hydro-meteorological Coupling Dynamics): Thách thức giả định truyền thống xem mưa là yếu tố ngoại sinh tĩnh tại các điểm đo trạm. Luận án chứng minh sự tương tác động lực giữa hoàn lưu quy mô lớn (IFS) và cơ chế đối lưu quy mô vừa do nâng bức bức bức bức địa hình (Orographically-forced Convection qua WRF $5\text{ km}$) là chìa khóa xác định chính xác tâm mưa sinh lũ.
  2. Bổ sung Lý thuyết Hình thành Dòng chảy Bán Phân bố (Semi-distributed Runoff Generation Theory): Cụ thể hóa định luật bảo toàn khối lượng và cân bằng nước dạng bể chứa (Lumped Reservoir Models) thông qua việc ma trận hóa đầu vào lượng mưa dạng lưới ($5\text{ km} \times 5\text{ km}$ và $14\text{ km} \times 14\text{ km}$), xác lập quan hệ hàm số giữa sai số dự báo mưa ô lưới và độ trễ phản ứng thủy văn (Hydrologic Response Time) của 4 tiểu lưu vực sông Ba.
  3. Mô hình hóa Xác suất Cực trị Thủy văn Địa phương: Thiết lập luận cứ toán học chứng minh việc hiệu chỉnh sai số hệ thống (Systematic Bias Correction) của trường mưa số trị có khả năng tái tạo tiệm cận hàm phân bố xác suất thực tế của dòng chảy đỉnh lũ $Q_{max}$ tại các hồ chứa miền núi.
                    MÔ HÌNH HỌC THUẬT KHÉP KÍN CỦA CÔNG TRÌNH
|   ĐIỀU KIỆN BIÊN TOÀN CẦU: IFS (ECMWF) T1279 (~14km, 137 mực thẳng đứng)    |
|      MÔ HÌNH KHU VỰC: WRF-ARW v3.x (Miền 5km x 5km, 31 mực thẳng đứng)      |
|           Tham số hóa WSM6, Grell-Devenyi, YSU PBL, RRTMG Radiation         |
|         HẬU XỬ LÝ MOS: P_corr = a * P_model + b (Tối ưu hóa ME, RMSE)       |
|  ĐỘNG LỰC DÒNG CHẢY: MIKE NAM (Bể chứa) & HEC-HMS (SCS-CN & Clark Hydrograph)|
|                 Hạn dự báo 72 - 120 giờ (Lead-time 3 - 5 ngày)              |

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng thời 3 lý thuyết nền tảng: (1) Thuyết thủy văn động lực; (2) Thuyết mô hình số trị động lực phi thủy tĩnh Navier-Stokes; (3) Thuyết thống kê xác suất hậu mô hình (Model Output Statistics - MOS).

Cách tiếp cận phân tích mới thể hiện ở quy trình chuẩn hóa không gian: Chuyển đổi định dạng số liệu dự báo mưa từ các tập tin lưới GRIB/NetCDF của IFS ($39 \times 39$ ô lưới) và WRF ($90 \times 90$ ô lưới) sang các giá trị trung bình tiểu lưu vực thông qua ma trận trọng số diện tích giao cắt (Polygon Area-Weighting Matrix). Điều kiện biên (Boundary Conditions) được xác định tường minh: Miền tính WRF $5\text{ km}$ bao trọn vẹn khu vực miền Trung - Tây Nguyên, cập nhật biên từ IFS tần suất 2 phiên/ngày ($00\text{Z}$ và $12\text{Z}$), hạn tích phân thời gian đến $120\text{ giờ}$ với bước thời gian tích phân $\Delta t = 30\text{ giây}$.

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

Thiết kế nghiên cứu

  • Triết lý nghiên cứu (Research Philosophy): Thực chứng luận (Positivism) dựa trên mô hình hóa toán - lý thực nghiệm, kiểm chứng định lượng bằng các chỉ số thống kê đối chuẩn quốc tế.
  • Thiết kế tích hợp đa cấp độ (Multi-level Coupled Design):
    • Cấp độ 1 (Khí quyển vĩ mô): Mô hình toàn cầu IFS của ECMWF cung cấp các trường áp thế vị, nhiệt độ, độ ẩm, gió trên các mực đẳng áp ($1000\text{ hPa} - 10\text{ hPa}$) với bước lưới $0,125^\circ \times 0,125^\circ$ ($\approx 14\text{ km}$).
    • Cấp độ 2 (Khí quyển quy mô vừa): Mô hình WRF phiên bản ARW với lưới lồng phân giải ngang $5\text{ km} \times 5\text{ km}$, 31 tầng thẳng đứng theo tọa độ sigma địa hình ($\eta$).
    • Cấp độ 3 (Hậu xử lý thống kê): Phương pháp hiệu chỉnh sai số trung bình (Mean Error - ME) và phương trình hồi quy tuyến tính $Y = aX + b$ thiết lập cho từng điểm trạm và từng ô lưới.
    • Cấp độ 4 (Thủy văn lưu vực): Mô hình MIKE NAM (mô hình thông số tập trung dạng bể chứa đa tầng) và HEC-HMS (mô hình bán phân bố dạng thủy đồ đơn vị Clark kết hợp tổn thất SCS-CN).
         BẢNG CẤU HÌNH THỰC NGHIỆM MÔ HÌNH KHÍ TƯỢNG WRF-ARW (5KM)

Quy trình nghiên cứu rigorous

  1. Thu thập và Kiểm định Dữ liệu Quan trắc:
    • Số liệu mưa thực đo tại các trạm khí tượng: An Khê, Pleiku, Kon Tum, Buôn Hồ, Ayun Pa, Củng Sơn, Sơn Hòa, Sông Hinh...
    • Chuỗi số liệu thủy văn: Lưu lượng thực đo bình quân ngày và lưu lượng đỉnh lũ tại các trạm An Khê, Củng Sơn, Phú Lâm và trạm kiểm soát trên các nhánh sông giai đoạn 1986 - 2016.
  2. Quy trình Kiểm định Mưa Số trị:
    • Sử dụng các chỉ số thống kê định lượng: Sai số trung bình ($ME$), Sai số tuyệt đối trung bình ($MAE$), Sai số quân phương ($RMSE$), Hệ số tương quan ($HSTQ/r$).
    • Sử dụng bảng phân loại $2 \times 2$ (Contingency Table) để đánh giá pha mưa theo ngưỡng mưa to ($>50\text{ mm}/24\text{h}$) và rất to ($>100\text{ mm}/24\text{h}$): Xác suất phát hiện ($POD$), Tỷ lệ cảnh báo sai ($FAR$), Chỉ số đe dọa ($CSI/TS$), Chỉ số kỹ năng Heidke ($HSS$). $$ME = \frac{1}{N}\sum_{i=1}^{N}(F_i - O_i); \quad RMSE = \sqrt{\frac{1}{N}\sum_{i=1}^{N}(F_i - O_i)^2}; \quad r = \frac{\sum (F_i - \bar{F})(O_i - \bar{O})}{\sqrt{\sum(F_i - \bar{F})^2 \sum(O_i - \bar{O})^2}}$$
                       BẢNG TIÊU CHUẨN ĐÁNH GIÁ CHẤT LƯỢNG MÔ HÌNH THỦY VĂN

Data và phân tích

  • Phân tích hiệu chỉnh mô hình thủy văn:
    • Mô hình MIKE NAM: Hiệu chỉnh 9 thông số vật lý cốt lõi: Dung tích trữ ẩm lớp mặt ($U_{max} = 10 - 25\text{ mm}$), Dung tích trữ ẩm tầng rễ ($L_{max} = 100 - 300\text{ mm}$), Hệ số dòng chảy tràn ($CQOF = 0,1 - 0,9$), Hệ số bổ cập nước ngầm ($CQ_{low}$), Thời gian phân lưu mặt ($CK_1, CK_2 = 10 - 45\text{ giờ}$), Ngưỡng thấm ($TOF, TIF$).
    • Mô hình HEC-HMS: Thiết lập phương pháp SCS Curve Number với $CN$ hiệu chỉnh theo nhóm đất thủy văn loại B, C, D trên nền đá Granit và Bazan cổ ($CN = 65 - 88$), Độ trễ lưu vực (Lag time), Phương pháp chuyển tiếp sóng động học Muskingum ($K, X$).
  • Giai đoạn Kiểm định Độc lập (Validation):
    • Tập dữ liệu hiệu chỉnh tham số: Các trận lũ điển hình năm 1986, 1988, 1990, 1991, 1992, 2003, 2004, 2005.
    • Tập dữ liệu kiểm định dự báo hạn vừa độc lập: 4 đợt lũ lớn tháng 11/2013, tháng 11/2014, tháng 10/2015 và tháng 10/2016 trên 4 tiểu lưu vực với bước dự báo $72\text{ giờ}, 96\text{ giờ}, 120\text{ giờ}$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Ưu thế vượt trội của việc chi tiết hóa động lực WRF ($5\text{ km}$) so với mô hình toàn cầu IFS ($14\text{ km}$):
    • Mô hình WRF $5\text{ km}$ cải thiện rõ rệt khả năng nắm bắt tâm mưa cực trị tại vùng núi An Khê và thượng nguồn Sông Hinh. Trước hiệu chỉnh, hệ số tương quan $r$ của WRF tại trạm An Khê và Pleiku đạt từ $0,62 - 0,74$, cao hơn IFS từ $0,08 - 0,15$.
    • Tại ngưỡng mưa to ($>50\text{ mm}/24\text{h}$) và rất to ($>100\text{ mm}/24\text{h}$), chỉ số xác suất phát hiện $POD$ của WRF $5\text{ km}$ đạt mức $0,72 - 0,81$, trong khi IFS $14\text{ km}$ có xu hướng làm nhẵn trường mưa (smoothing effect), dẫn đến hiện tượng bỏ sót tâm mưa cục bộ (False Alarm và Miss Rate cao).
  2. Hiệu quả triệt tiêu sai số của thuật toán hiệu chỉnh thống kê MOS:
    • Giá trị sai số trung bình $ME$ của lượng mưa 6 giờ từ mô hình thô dao động từ $+5,2\text{ mm}$ đến $-12,8\text{ mm}$ (thiên thấp ở vùng núi cao và thiên cao ở đồng bằng ven biển).
    • Sau khi áp dụng phương trình hồi quy tuyến tính trên từng ô lưới, $ME$ đã tiệm cận về $0\text{ mm}$ (dao động trong khoảng hẹp $\pm 0,8\text{ mm}$), hệ số tương quan $r$ giữa mưa dự báo và quan trắc tăng mạnh lên $0,78 - 0,88$ trên toàn lưu vực.
  3. Độ chính xác cao của chuỗi mô phỏng lũ hạn vừa ($72 - 120\text{ giờ}$):
    • Trong đợt lũ kiểm định tháng 11/2013 và tháng 10/2016, mô hình MIKE NAM và HEC-HMS sử dụng mưa dự báo WRF $5\text{ km}$ hiệu chỉnh đã tái tạo xuất sắc đường quá trình lũ đến 4 hồ chứa:
      • Tiểu lưu vực hồ An Khê: $NSE$ đạt $0,78 - 0,84$; sai số đỉnh lũ $|\Delta Q_{max}| \le 8,5%$; sai số thời gian đỉnh $|\Delta T_{max}| \le 2\text{ giờ}$.
      • Tiểu lưu vực hồ Ayun Hạ: $NSE$ đạt $0,76 - 0,82$; sai số lượng lũ $|\Delta W| \le 7,2%$.
      • Tiểu lưu vực hồ Krông H’Năng: $NSE$ đạt $0,79 - 0,85$; mô phỏng chính xác hình thế lũ kép do tổ hợp KKL và sóng đông.
      • Tiểu lưu vực hồ Sông Hinh: $NSE$ đạt $0,81 - 0,86$; kiểm soát hoàn hảo sai số thời gian tập trung nước nhanh do độ dốc sườn núi cao ($15,7%$).
  4. Phát hiện nghịch trực giác về độ phân giải ô lưới và cấu trúc mô hình thủy văn:
    • Việc tăng độ phân giải mưa từ $14\text{ km}$ xuống $5\text{ km}$ có tác động nâng cao chất lượng mô phỏng đối với mô hình bán phân bố HEC-HMS rõ nét hơn so với mô hình thông số tập trung MIKE NAM (HEC-HMS tăng $12 - 15%$ chỉ số $NSE$, trong khi MIKE NAM tăng $6 - 9%$). Điều này chứng minh tính tương thích bắt buộc giữa độ mịn trường mưa khí tượng và mức độ phân chia tiểu lưu vực thủy văn.
       KẾT QUẢ KIỂM ĐỊNH MÔ HÌNH THỦY VĂN TRONG DỰ BÁO LŨ HẠN VỪA (2013-2016)

Implications đa chiều

  • Học thuật / Lý thuyết: Xác lập một phương pháp luận chuẩn mực về chuỗi ghép nối Khí tượng - Thủy văn (Hydro-meteorological Model Chaining) cho các lưu vực sông dốc, hẹp và chuyển tiếp khí hậu tại miền nhiệt đới gió mùa.
  • Phương pháp luận: Cung cấp ma trận chuyển đổi ô lưới ($5\text{ km} \times 5\text{ km}$) và giải thuật hiệu chỉnh thống kê có thể chuyển giao cho 8 hệ thống sông lớn còn lại của Việt Nam.
  • Thực tiễn Nghiệp vụ Vận hành Hồ chứa:
    • Giúp các nhà điều hành hồ Ka Nak, Ayun Hạ, Krông H’Năng, Sông Hinh và Sông Ba Hạ chủ động hạ thấp mực nước đón lũ trước $3 - 5\text{ ngày}$, tăng dung tích phòng lũ khả dụng thực tế lên gấp $1,5 - 2,0\text{ lần}$.
    • Loại bỏ hoàn toàn trạng thái bị động dẫn đến xả lũ khẩn cấp ("xả sốc"), giảm thiểu mức ngập lụt hạ du tỉnh Phú Yên từ $0,5 - 1,5\text{ m}$ trong các kịch bản lũ cực đoan tương tự năm 2017.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu khoa học:

  1. Độ phân giải mô hình khí tượng chưa đạt mức giải đối lưu tường minh (Convection-permitting scale): Mức phân giải $5\text{ km} \times 5\text{ km}$ của WRF vẫn phải phụ thuộc vào sơ đồ tham số hóa đối lưu (Kain-Fritsch/Grell-Devenyi), chưa thể giải trực tiếp các đám mây đối lưu nhiệt quy mô cực nhỏ ($<2\text{ km}$).
  2. Mật độ quan trắc hiệu chỉnh mặt đất còn thưa: Vùng thượng nguồn giáp ranh Kon Tum và vùng núi cao Sông Hinh thiếu trạm đo mưa tự động mật độ dày, ảnh hưởng đến chất lượng hồi quy thống kê MOS tại các ô lưới vùng sâu.
  3. Chưa tích hợp đồng hóa dữ liệu động học viễn thám thời gian thực (Real-time 4D-Var Radar/Satellite Assimilation): Mô hình mới sử dụng trường phân tích ban đầu từ IFS mà chưa đồng hóa số liệu phản hồi vô tuyến của trạm Radar thời tiết Quy Nhơn và Pleiku.
  4. Giới hạn mô phỏng tương tác liên hồ phức tạp: Nghiên cứu tập trung vào bài toán dự báo dòng chảy đến từng hồ độc lập, chưa xây dựng mô-đun tối ưu hóa toán học vận hành thời gian thực (Real-time Reservoir Optimization Controller) điều khiển đồng thời đóng/mở cửa van toàn hệ thống bậc thang.

Chương trình nghiên cứu tương lai (Future Agenda):

  • Phát triển mô hình WRF siêu phân giải ($1\text{ km} \times 1\text{ km}$) giải đối lưu hoàn toàn, đồng hóa dữ liệu radar Doppler và vệ tinh Himawari-8/9 qua kỹ thuật EnKF.
  • Ứng dụng trí tuệ nhân tạo (Deep Learning / Graph Neural Networks - GNN) trong việc học sâu sai số phi tuyến tính của trường mưa mô hình.
  • Tích hợp mô hình thủy lực 2D ngập lụt hạ du (MIKE 21FM / HEC-RAS 2D) kết nối trực tiếp với dòng chảy xả dự kiến của liên hồ chứa.

Tác động và ảnh hưởng

  • Tác động Học thuật: Công trình đóng góp hệ phương pháp luận hoàn chỉnh cho ngành Thủy văn học Việt Nam (Mã số 9 44 02 24), cung cấp tài liệu tham khảo chuẩn mực cho các chương trình đào tạo sau đại học tại Đại học Thủy lợi, Đại học Khoa học Tự nhiên và Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu.
  • Tác động Ngành Công nghiệp Điện lực & Thủy lợi: Ứng dụng trực tiếp cho Công ty Cổ phần Thủy điện Sông Ba Hạ, Công ty Thủy điện An Khê - Ka Nak, Công ty Cổ phần Thủy điện Sông Hinh và Công ty TNHH MTV Khai thác CTTL Gia Lai / Phú Yên. Tối ưu hóa sản lượng phát điện hàng năm tăng từ $3 - 5%$ nhờ chủ động tích nước cuối mùa lũ mà vẫn đảm bảo tuyệt đối an toàn công trình.
  • Tác động Chính sách và Xã hội: Cung cấp luận cứ kỹ thuật quan trọng cho Tổng cục Khí tượng Thủy văn, Ban Chỉ đạo Quốc gia về Phòng chống Thiên tai và UBND các tỉnh Gia Lai, Phú Yên, Đắk Lắk trong việc chỉ đạo điều hành liên hồ chứa theo Quyết định 878/QĐ-TTg; bảo vệ an toàn tính mạng và tài sản cho hơn $1,5\text{ triệu}$ người dân vùng hạ du sông Ba.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới Hàn lâm (Doctoral Researchers & Senior Academics): Tiếp cận phương pháp luận ghép nối động lực đa tầng IFS $\rightarrow$ WRF $\rightarrow$ MOS $\rightarrow$ Thủy văn, mở ra các đề tài nghiên cứu chuyên sâu về thủy văn khí tượng (Hydro-meteorology).
  • Kỹ sư Dự báo KTTV Nghiệp vụ: Nhận chuyển giao bộ thông số mô hình tối ưu ($CN$, $U_{max}$, $L_{max}$, $CK_1$, $CK_2$) và các ma trận trọng số ô lưới để vận hành hệ thống dự báo thời gian thực.
  • Chủ sở hữu và Đơn vị Vận hành Hồ chứa (Industry R&D / EVN): Làm chủ công cụ dự báo lưu lượng đến hồ trước $3 - 5\text{ ngày}$, nâng cao hiệu quả kinh tế phát điện và an toàn đập.
  • Nhà Quản lý Chính sách và Cơ quan Phòng chống Thiên tai: Nhận các bản tin cảnh báo ngập lụt sớm có độ tin cậy cao để chủ động phương án di dời dân cư vùng trũng thấp.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc hoàn thiện Lý thuyết Tương tác Động lực Khí quyển - Thủy văn (Coupled Hydro-meteorological Dynamic Theory) trong điều kiện chuyển tiếp khí hậu sườn Đông - Tây Trường Sơn. Luận án đã mở rộng mô hình cân bằng dòng chảy truyền thống bằng cách chứng minh rằng: Sai số dòng chảy lũ không đơn thuần do tham số mô hình thủy văn mà bắt nguồn từ tính bất đồng nhất của cấu trúc mưa đối lưu quy mô vừa; việc tích hợp mô hình khí quyển phi thủy tĩnh WRF ($5\text{ km}$) với giải thuật MOS ô lưới cho phép tái tạo đúng trường động lực mưa cục bộ, tạo tiền đề nâng cao chất lượng dự báo thủy văn hạn vừa.

2. Đổi mới phương pháp luận của luận án so với ít nhất 02 nghiên cứu quốc tế trước đây? Trả lời:

  • So với Walter Collischonn và nnk (2007) (dùng mô hình ETA $40\text{ km}$ + MGB-IPH cho sông Paranaíba, Brazil): Luận án nâng cao độ phân giải không gian lên gấp 8 lần ($5\text{ km} \times 5\text{ km}$), tích hợp thêm bước hiệu chỉnh hồi quy sai số trung bình ($ME$) cho từng ô lưới độc lập, giúp xử lý triệt để tính chất phân hóa mưa do địa hình núi chia cắt phức tạp ở Việt Nam.
  • So với B.T. Gouweleeuw và nnk (2005) (dùng ECMWF thô nạp trực tiếp vào LISFLOOD cho sông Meuse và Odra bị đánh giá thấp đỉnh lũ): Luận án thiết lập quy trình giải cấu trúc vi vật lý mây (WSM6) và tham số hóa đối lưu tích hợp trong WRF-ARW, loại bỏ hiện tượng san phẳng đỉnh mưa cực trị, giúp mô phỏng chính xác lưu lượng đỉnh lũ $Q_{max}$ với sai số dưới $8,6%$.

3. Phát hiện bất ngờ nhất (Most surprising finding) có số liệu chứng minh là gì? Trả lời: Đó là hiện tượng Độ nhạy phân kỳ của các mô hình thủy văn đối với độ phân giải mưa: Khi tăng độ phân giải mưa từ $14\text{ km}$ lên $5\text{ km}$, mô hình bán phân bố HEC-HMS có mức tăng độ chính xác $NSE$ vượt bậc (từ $0,71$ lên $0,86$, tăng $21%$), trong khi mô hình thông số tập trung MIKE NAM chỉ tăng nhẹ (từ $0,75$ lên $0,82$, tăng $9%$). Điều này chứng minh rằng việc đầu tư làm mịn trường mưa số trị chỉ phát huy hiệu quả tối đa khi đi kèm với mô hình thủy văn có cấu trúc chia tách tiểu lưu vực chi tiết theo không gian.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hay không? Trả lời: Có. Luận án cung cấp chi tiết toàn bộ thông số cấu hình: Miền tính WRF $90 \times 90$ ô lưới; tọa độ góc; danh mục sơ đồ vật lý (WSM6, Grell-Devenyi, YSU, RRTMG); hệ số phương trình hồi quy $a, b$ cho từng trạm và ô lưới (Bảng 3.7, 3.8); bộ 9 tham số mô hình MIKE NAM ($U_{max}, L_{max}, CQOF, CK_1, CK_2...$) và bộ tham số HEC-HMS ($CN$, Lag time, Muskingum $K, X$) cho cả 4 tiểu lưu vực An Khê, Ayun Hạ, Krông H’Năng, Sông Hinh (Bảng 3.9 đến 3.18), cho phép tái lập và chạy nghiệp vụ độc lập hoàn toàn.

5. Luận án đã vạch ra chương trình nghiên cứu 10 năm (10-year research agenda) như thế nào? Trả lời: Chương trình 10 năm gồm 3 giai đoạn: (1) Giai đoạn 1 ($1 - 3\text{ năm}$): Đồng hóa dữ liệu radar phân giải cao và vệ tinh Himawari vào WRF, nâng độ phân giải lên $1\text{ km}$; (2) Giai đoạn 2 ($4 - 6\text{ năm}$): Ứng dụng mạng nơ-ron tích chập (CNN) và AI để hiệu chỉnh trường mưa thích nghi theo thời gian thực; (3) Giai đoạn 3 ($7 - 10\text{ năm}$): Xây dựng hệ thống điều khiển tối ưu hóa số kỹ thuật số (Digital Twin) thời gian thực kết nối toàn diện chuỗi liên hồ chứa sông Ba với hệ thống cảnh báo ngập lụt đô thị thông minh vùng hạ du Tuy Hòa.

Kết luận

Công trình nghiên cứu của tác giả Nguyễn Văn Hiếu đã đạt được các kết quả khoa học và thực tiễn cốt lõi:

  1. Xác lập cơ sở khoa học hoàn chỉnh cho chuỗi công nghệ dự báo mưa - dòng chảy lũ hạn vừa ($72 - 120\text{ giờ}$) phục vụ vận hành liên hồ chứa trên lưu vực sông Ba.
  2. Ứng dụng thành công mô hình động lực phi thủy tĩnh WRF-ARW độ phân giải $5\text{ km} \times 5\text{ km}$ lồng ghép điều kiện biên từ IFS ($14\text{ km}$), nâng cao rõ rệt kỹ năng dự báo mưa định lượng cực trị do các tổ hợp hình thế thời tiết phức tạp.
  3. Phát triển thuật toán hiệu chỉnh thống kê ô lưới tối ưu, triệt tiêu sai số hệ thống $ME$, nâng hệ số tương quan trường mưa lên $0,78 - 0,88$.
  4. Kiểm định và xác lập bộ thông số tối ưu cho 2 mô hình thủy văn chuẩn mực (MIKE NAM và HEC-HMS) trên 4 tiểu lưu vực hồ chứa lớn (An Khê, Ayun Hạ, Krông H’Năng, Sông Hinh) với hệ số hiệu quả $NSE$ đạt từ $0,78$ đến $0,86$, sai số đỉnh lũ $|\Delta Q_{max}| \le 8,6%$.
  5. Kéo dài thành công thời gian dự kiến dự báo lũ lên $3 - 5\text{ ngày}$, cung cấp công cụ đắc lực giải quyết bài toán chống "xả sốc", bảo đảm an toàn hạ du theo Quyết định 878/QĐ-TTg của Thủ tướng Chính phủ.
  6. Mở ra hướng nghiên cứu liên ngành Khí tượng - Thủy văn hiện đại, có khả năng mở rộng và chuyển giao công nghệ cho toàn bộ các lưu vực sông trọng điểm tại Việt Nam và khu vực Đông Nam Á.