Tổng quan về luận án

Mô hình hóa và dự báo dòng chảy lũ trên các lưu vực sông có vai trò sống còn trong quản lý tài nguyên nước, thiết kế công trình thủy lợi – thủy điện và giảm thiểu rủi ro thiên tai. Tuy nhiên, tại các lưu vực sông thuộc khu vực Nam Trung Bộ Việt Nam, mạng lưới quan trắc khí tượng thủy văn còn rất thưa thớt, thiếu hụt nghiêm trọng số liệu mặt cắt hình học lòng dẫn và số liệu mưa phân bố không gian có độ phân giải cao. Luận án tiến sĩ Thủy văn học của tác giả Bùi Văn Chanh với đề tài "Nghiên cứu cải tiến mô hình MARINE để mô phỏng và dự báo dòng chảy cho lưu vực sông thiếu số liệu - Áp dụng cho khu vực Nam Trung Bộ" (Mã số: 9440224, bảo vệ năm 2022 tại Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu, dưới sự hướng dẫn khoa học của PGS.TS Trần Ngọc Anh và PGS.TS Lương Tuấn Anh) là công trình nghiên cứu tiên phong giải quyết bài toán mô phỏng dòng chảy phân bố trong điều kiện thiếu dữ liệu thực đo.

Nghiên cứu tập trung giải quyết khoảng trống học thuật cốt lõi (research gap): Mô hình thủy văn thông số phân bố vật lý MARINE gốc (Estupina, 2000; HydroEco - IMFT) dù mô tả xuất sắc quá trình dòng chảy sườn dốc nhưng tồn tại ba hạn chế lớn khi ứng dụng thực tế: (1) chỉ diễn toán dòng chảy trên một trục sông chính thông qua thủ tục cộng dồn lưu lượng CalcApLat, xem dòng chảy sông nhánh như dòng chảy sườn dốc; (2) hoàn toàn khuyết thiếu mô đun điều tiết dòng chảy qua công trình hồ chứa; và (3) phụ thuộc vào radar/vệ tinh vốn có sai số lớn tại vùng núi hoặc sử dụng phép chia đa giác Thiessen làm triệt tiêu tính biến thiên không gian của mưa. Luận án xây dựng và kiểm chứng hai câu hỏi nghiên cứu cùng hai giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa chính xác quá trình truyền sóng lũ trong mạng lưới sông phân nhánh phức tạp và qua hệ thống hồ chứa khi không có số liệu mặt cắt ngang đo đạc chi tiết?
  • Câu hỏi nghiên cứu 2 (RQ2): Phương pháp nội suy mưa không gian kết hợp gradient độ cao nào tối ưu hóa được dữ liệu đầu vào cho mô hình phân bố trong điều kiện địa hình chia cắt mạnh?
  • Giả thuyết 1 (H1): Việc tích hợp phương trình sóng động học 1 chiều (1D Kinematic Wave) giải bằng sơ đồ sai phân ẩn kết hợp thuật toán lặp phi tuyến Newton - Raphson và mô đun hồ chứa Runge - Kutta bậc 3 sẽ khắc phục hoàn toàn sự méo dạng đỉnh lũ của mô hình MARINE gốc.
  • Giả thuyết 2 (H2): Phương pháp nghịch đảo khoảng cách theo trọng số (IDW) kết hợp mô hình hồi quy bội phi tuyến theo độ cao địa hình sẽ nâng cao vượt bậc độ chính xác của trường mưa phân bố đầu vào so với phương pháp đa giác truyền thống.

Khung lý thuyết của luận án được xây dựng dựa trên sự liên kết giữa lý thuyết thủy lực Saint - Venant 2 chiều xấp xỉ sóng động học cho dòng chảy tràn, mô hình thấm bão hòa Green - Ampt, lý thuyết sóng động học 1 chiều (Lighthill & Whitham, 1955) cho mạng lưới sông và phương trình bảo toàn khối lượng hồ chứa giải bằng Runge - Kutta bậc 3. Phạm vi nghiên cứu bao quát 3 lưu vực sông trọng điểm Nam Trung Bộ không chịu ảnh hưởng triều: Sông Cái Nha Trang ($2000\text{ km}^2$, chiều dài sông chính $79\text{ km}$, độ dốc $3,7‰$), Sông Dinh Ninh Hòa ($1000\text{ km}^2$, đỉnh Chư Mu cao $2000\text{ m}$) và Sông Cái Phan Rang ($3000\text{ km}^2$, chiều dài $119\text{ km}$, độ dốc lưu vực $17,7‰$). Kết quả đã nâng cao chỉ số Nash-Sutcliffe (NSE) từ mức dưới $0,65$ của mô hình gốc lên trên $0,80 - 0,91$ trong các trận lũ lịch sử (2009, 2010, 2013, 2016, 2017, 2018) và thiết lập thành công quy trình dự báo nghiệp vụ thời gian thực hạn $24\text{h} - 48\text{h}$ trong mùa lũ năm 2020.

Literature Review và Positioning

Lịch sử phát triển mô hình toán thủy văn ghi nhận bước chuyển dịch mang tính thế kỷ từ các mô hình tuyến tính thông số tập trung sang các hệ thống mô hình thông số phân bố dựa trên cơ sở vật lý thực nghiệm (Biswas, 1970; Todini, 1988; Singh, 1995). Các mô hình thông số tập trung truyền thống như Tank (Sugawara, 1961), NAM (DHI, 1972), SSARR (US Army Corps of Engineers) hay HEC-HMS (Feldman, 2000) đã bộc lộ những giới hạn nghiêm trọng khi ứng dụng cho các lưu vực có địa hình chia cắt phức tạp. Như nhận định của Gupta và cộng sự (2008) cũng như Beven (2012), mô hình tập trung đã "trung bình hóa, tham số hóa nhiều yếu tố tự nhiên của lưu vực, làm suy giảm tính liên kết vật lý giữa quá trình mưa - dòng chảy và không phản ánh được tính phi đồng nhất không gian của mặt đệm".

Để khắc phục nhược điểm này, các dòng mô hình bán phân bố và phân bố đã ra đời, phân chia thành hai trường phái chính:

  1. Mô hình phân bố tham số/khái niệm: Điển hình là TOPMODEL (Beven & Kirkby, 1979), SWAT (Arnold et al., 1998), IFAS (ICHARM, Nhật Bản), FRASC và DIMOSOP (Todini, 1996). Nhóm mô hình này mô phỏng dòng chảy sườn dốc thông qua các cụm bể chứa phân tầng (tank cascade) trên từng ô lưới hoặc đơn vị phản ứng nhóm GRU (Grouped Response Unit). Tuy nhiên, các mô hình này phụ thuộc nặng nề vào các thông số kinh nghiệm và gặp khó khăn khi mô tả tương tác vật lý trực tiếp giữa dòng mặt và mạng lưới sông hở.
  2. Mô hình phân bố dựa trên vật lý cơ bản (Physically-Based Distributed Models): Tiêu biểu là MIKE SHE (Abbott et al., 1986), WETSPA (Wang et al., 1996), CATHY (Camporese et al., 2010) và MARINE (HydroEco - IMFT; Estupina, 2000; Dartus & Labat, 2005). Các mô hình này thiết lập các phương trình vi phân đạo hàm riêng mô tả cơ chế thủy lực dòng chảy mặt và phương trình thấm Richards hoặc Green - Ampt.

Trong bức tranh tổng quan đó, Dự án So sánh các Mô hình Thủy văn Thông số Phân bố (DMIP - Distributed Model Intercomparison Project do Smith et al., 2004, 2012 chủ trì) đã khẳng định: "Biến đổi lượng mưa và đặc trưng vật lý của các yếu tố tự nhiên trên lưu vực là lý do sử dụng mô hình thủy văn thông số phân bố thay cho mô hình thủy văn thông số tập trung". Tuy nhiên, rào cản lớn nhất khi ứng dụng mô hình phân bố là tình trạng khan hiếm dữ liệu. Đối với các lưu vực thiếu số liệu (PUB - Predictions in Ungauged Basins), Tim Van Emmerik và cộng sự (2015) định nghĩa rằng: "Lưu vực thiếu số liệu quan trắc là lưu vực thiếu đáng kể dữ liệu thủy văn tại chỗ như lượng mưa, dòng chảy, bốc hơi theo thời gian".

Về mặt học thuật, tồn tại cuộc tranh luận lớn giữa hai quan điểm: (1) Sử dụng các mô hình thủy lực phức tạp 1D/2D (như HEC-RAS, MIKE 11, Telemac-2D) ghép nối với mô hình mưa - dòng chảy; và (2) Sử dụng các phương pháp diễn toán thủy văn tinh giản (như Muskingum, Muskingum - Cunge). Quan điểm thứ nhất thất bại ở các lưu vực vùng núi thiếu số liệu mặt cắt địa hình thực đo, trong khi quan điểm thứ hai bộc lộ khiếm khuyết vì giả định tuyến tính không phản ánh được sự biến dạng phi tuyến của sóng lũ lớn.

So sánh với các nghiên cứu quốc tế, Ruelland và cộng sự (2008) khi đánh giá các phương pháp nội suy không gian tại Tây Phi đã chứng minh IDW và Kriging vượt trội hơn hẳn phép chia đa giác Thiessen. Trong khi đó, nghiên cứu của M. Coulibaly (2010) tại lưu vực Charlie Creek (Hoa Kỳ) và Nam Phi chỉ ra rằng việc kết hợp hồi quy địa hình với nội suy không gian giúp giảm sai số ước lượng mưa vùng núi trên 35%.

Luận án của NCS. Bùi Văn Chanh định vị chính xác vị trí khoa học của mình: Không phát triển một mô hình hoàn toàn mới từ đầu mà tiến hành cải tiến cấp độ hợp nhất (integrated coupling) trên nền tảng mã nguồn mở của mô hình MARINE, bổ sung mắt xích động lực học lòng dẫn và điều tiết hồ chứa vốn bị bỏ ngỏ, tạo nên một công cụ hoàn chỉnh, độc lập và tối ưu cho điều kiện thủy văn thực tế tại Việt Nam.

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

+-----------------------------------------------------------------------------+
|                  KHUNG PHÂN TÍCH TÍCH HỢP MÔ HÌNH MARINE CẢI TIẾN           |
+-----------------------------------------------------------------------------+
                    [ DỮ LIỆU ĐẦU VÀO ĐỊA HÌNH & KHÍ TƯỢNG ]
  |  MÔ ĐUN 1: NỘI SUY MƯA KHÔNG GIAN (IDW + HỒI QUY PHI TUYẾN THEO ĐỘ CAO) |
  |  - Phân tích tương quan mưa - độ cao địa hình: P = f(Z)                 |
  |  - Trọng số nghịch đảo khoảng cách IDW có bán kính tìm kiếm & vùng chắn |
  |  - Tạo trường mưa phân bố chi tiết theo từng ô lưới (Grid-based Rain)   |
  |  MÔ ĐUN 2: THỦY VĂN SƯỜN DỐC (MÔ HÌNH MARINE GỐC CẢI TIẾN)              |
  |  - Tính thấm lớp đất mặt bão hòa / chưa bão hòa: Cơ chế Green-Ampt      |
  |  - Dòng chảy tràn sườn dốc 2D: Phương trình Saint-Venant xấp xỉ         |
  |  - Trích xuất dòng chảy sườn dốc đổ vào từng đoạn sông (Lateral Inflow) |
  |  MÔ ĐUN 3: DIỄN TOÁN MẠNG LƯỚI SÔNG (SÓNG ĐỘNG HỌC 1D PHI TUYẾN)        |
  |  - Phương trình liên tục & phương trình động lượng sóng động học 1D      |
  |  - Sơ đồ sai phân ẩn + Thuật toán giải phi tuyến Newton-Raphson         |
  |  - Cấu trúc mạng lưới sông phân cấp (không cần mặt cắt ngang đo đạc)    |
  | MÔ ĐUN 4: DIỄN TOÁN HỒ CHỨA      |    | DÒNG CHẢY HẠ DU TỰ DO            |
  | - Phương trình cân bằng nước    |    | (Non-reservoir reaches)          |
  | - Quan hệ lòng hồ VFZ & tràn xả |    |                                  |
  | - Tích phân Runge-Kutta bậc 3   |    |                                  |
  |                  LƯU LƯỢNG DÒNG CHẢY TẠI CÁC VỊ TRÍ KIỂM SOÁT           |
  |         (Đồng Trăng, Ninh Hòa, Tân Mỹ... & Dự báo nghiệp vụ 24h/48h)     |

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

Luận án mở rộng và làm sâu sắc thêm ba lý thuyết nền tảng trong thủy văn học và thủy lực học:

  1. Lý thuyết sóng động học (Kinematic Wave Theory - Lighthill & Whitham, 1955): Nghiên cứu mở rộng ứng dụng của lý thuyết sóng động học một chiều từ các kênh hở đơn lẻ sang toàn bộ mạng lưới sông phân nhánh tự nhiên có bậc sông phức tạp. Công trình chứng minh rằng với các con sông vùng núi Nam Trung Bộ có độ dốc đáy lớn ($S_0 > 0,002$), số Froude thỏa mãn điều kiện ổn định, các số hạng quán tính và gradient áp lực trong phương trình Saint-Venant đầy đủ là hoàn toàn có thể bỏ qua mà không làm suy giảm độ chính xác của đường quá trình lũ.
  2. Lý thuyết thấm Green - Ampt (1911): Tích hợp phân bố tham số thổ nhưỡng theo không gian (độ dẫn thủy lực bão hòa $K_s$, thế áp suất ẩm tại màng thấm $\psi$, độ xốp hiệu dụng $\theta_s - \theta_i$) với mô hình số độ cao DEM để mô tả chính xác quá trình bão hòa cục bộ và hình thành dòng chảy tràn theo cơ chế Dunne và Horton.
  3. Mệnh đề khoa học (Theoretical Proposition):
    • Mệnh đề 1 (P1): Dòng chảy sườn dốc quy mô ô lưới nhỏ hơn $1\text{ km}^2$ có thể mô hình hóa dưới dạng 2 chiều xấp xỉ, nhưng khi diện tích tích lũy vượt quá $1\text{ km}^2$, quá trình tập trung dòng chảy chuyển pha hoàn toàn sang động lực học lòng dẫn 1 chiều, đòi hỏi phải giải quyết tương tác mạng lưới thay vì lũy tích đại số.
    • Mệnh đề 2 (P2): Hiệu ứng cản và dung tích điều tiết của các hồ chứa bậc thang tạo ra sự phi tuyến tính cực đại trong đường truyền sóng lũ; việc tuyến tính hóa hồ chứa dẫn đến sai số pha đỉnh lũ lên tới $3 - 6\text{ giờ}$, điều này chỉ có thể khắc phục bằng sơ đồ tích phân phi tuyến cấp cao Runge - Kutta bậc 3.

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

Khung phân tích của luận án tích hợp đồng thời 4 trụ cột công nghệ và toán học:

  • Tích hợp mô đun nội suy mưa 2 thành phần: Kết hợp phép nghịch đảo khoảng cách IDW có giới hạn bán kính tìm kiếm và vùng chắn địa hình với phương trình hồi quy phi tuyến $P = a \cdot Z^2 + b \cdot Z + c$ (trong đó $P$ là lượng mưa giờ, $Z$ là cao trình địa hình từ DEM).
  • Thuật toán sóng động học phi tuyến mạng lưới: Sử dụng sơ đồ sai phân 4 điểm ẩn (Preissmann type), tuyến tính hóa hệ phương trình bằng thuật toán Newton - Raphson, cho phép giải đồng thời lưu lượng $Q$ và diện tích mặt cắt ướt $A$ trên toàn bộ các đoạn sông nhánh và dòng chính mà chỉ cần thông số độ dốc đáy và bề rộng lòng sông ước lượng qua ảnh viễn thám/DEM.
  • Mô đun diễn toán dòng chảy qua hồ chứa (Reservoir Module): Thiết lập phương trình vi phân cấp 1 về cân bằng thể tích: $$\frac{dV}{dt} = I(t) - O(H, t)$$ trong đó $I(t)$ là lưu lượng đến hồ (được cấp từ mô hình sóng động học nhánh thượng lưu), $O(H, t)$ là lưu lượng xả qua công trình tháo nước và tràn tự do, kết hợp đường đặc tính hình học lòng hồ $V = f(Z)$, $F = f(Z)$ (VFZ). Phương trình được giải chính xác bằng sơ đồ Runge - Kutta 3 điểm gia số ($\Delta H_1, \Delta H_2, \Delta H_3$).
  • Điều kiện biên và ranh giới áp dụng (Boundary Conditions): Mô hình vận hành tối ưu cho các lưu vực sông sườn dốc dốc, không chịu ảnh hưởng của sóng triều truyền ngược từ biển, và mật độ trạm đo mưa mặt đất đạt tối thiểu 1 trạm/$200 - 300\text{ km}^2$.

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

Thiết kế nghiên cứu

Luận án theo đuổi thế giới quan bản thể luận duy thực khách quan (Positivism / Critical Realism), tiếp cận bài toán thủy văn thông qua phương pháp mô hình hóa toán - lý có cấu trúc tiền định (deterministic physical modeling). Thiết kế nghiên cứu được thực hiện đa cấp độ:

  • Cấp độ 1: Mô hình hóa phân tán trên ô lưới sườn dốc (Raster Cell Level - kích thước ô lưới $100\text{ m} \times 100\text{ m}$ đến $250\text{ m} \times 250\text{ m}$).
  • Cấp độ 2: Mô hình hóa bán phân bố dọc theo các đoạn sông (River Reach Level).
  • Cấp độ 3: Mô hình hóa tập trung tại các nút công trình hồ chứa (Point/Node Level).

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm được chuẩn hóa qua 5 bước nghiêm ngặt:

  1. Xây dựng cơ sở dữ liệu GIS và Viễn thám:
    • Dữ liệu DEM phân giải cao xử lý từ Shuttle Radar Topography Mission (SRTM) và bản đồ địa hình tỷ lệ 1:50.000 để trích xuất mạng lưới sông, diện tích lưu vực và độ dốc sườn.
    • Bản đồ thổ nhưỡng thổ nhưỡng Feralit đỏ vàng, đất mùn vàng núi cao, đất phù sa được số hóa và gán các thông số vật lý: chiều sâu tầng đất ẩm $h_{soil}$, độ ẩm ban đầu $\theta_0$, độ dẫn bão hòa $K_s$.
    • Bản đồ thảm phủ thực vật trích xuất từ ảnh vệ tinh Landsat/Sentinel, mã hóa hệ số nhám sườn dốc Manning ($n_{slope}$ từ $0,035$ đến $0,45$) và nhám lòng dẫn ($n_{river}$ từ $0,025$ đến $0,045$).
  2. Kỹ thuật tích hợp mã nguồn: Luận án can thiệp trực tiếp vào mã nguồn C++/Fortran của mô hình MARINE, viết lại hoàn toàn cấu trúc hàm CalcApLat bằng thủ tục tính toán thủy lực Qluoisong, nhúng trực tiếp thư viện giải phương trình vi phân Runge - Kutta và hàm nội suy IDW vào luồng thực thi trung tâm.
  3. Hiệu chuẩn và kiểm định (Calibration & Validation): Áp dụng nguyên tắc chia tách chuỗi dữ liệu độc lập (Split-sample test). Bộ thông số được hiệu chỉnh tự động và bán thủ công trên các trận lũ lớn, sau đó kiểm định trên các trận lũ độc lập khác nhau về biên độ và thời gian xuất hiện.
  4. Hệ thống chỉ tiêu đánh giá độ chính xác:
    • Hệ số hiệu quả Nash - Sutcliffe (NSE): $$NSE = 1 - \frac{\sum_{t=1}^T (Q_{obs}^t - Q_{sim}^t)^2}{\sum_{t=1}^T (Q_{obs}^t - \bar{Q}_{obs})^2}$$
    • Sai số tương đối đỉnh lũ ($PE_p, %$): $$PE_p = \frac{Q_{p,sim} - Q_{p,obs}}{Q_{p,obs}} \times 100%$$
    • Sai số tương đối tổng lượng dòng chảy lũ ($PE_{Vol}, %$): $$PE_{Vol} = \frac{W_{sim} - W_{obs}}{W_{obs}} \times 100%$$
    • Sai số thời gian xuất hiện đỉnh lũ ($\Delta T_p, \text{giờ}$).
| BẢNG SO SÁNH HIỆU QUẢ MÔ PHỎNG: MÔ HÌNH MARINE GỐC VS MÔ HÌNH MARINE CẢI TIẾN             |
| Lưu vực / Trạm Thủy văn  | Trận Lũ Năm  | Mô hình     | NSE         | PEp (%)  | PEVol (%)|
| Sông Cái Nha Trang       | 2009         | MARINE Gốc  | 0.62        | -24.5    | +16.2    |
| (Trạm Đồng Trăng)        |              | CẢI TIẾN    | **0.88**    | **-4.2** | **+3.1** |
|                          | 2013         | MARINE Gốc  | 0.58        | -28.1    | +18.7    |
|                          |              | CẢI TIẾN    | **0.91**    | **+3.5** | **-2.4** |
|                          | 2016         | MARINE Gốc  | 0.65        | -21.0    | +14.3    |
|                          |              | CẢI TIẾN    | **0.86**    | **-5.1** | **+4.0** |
| Sông Dinh Ninh Hòa       | 2010         | MARINE Gốc  | 0.54        | -32.4    | +22.0    |
| (Trạm Ninh Hòa)          |              | CẢI TIẾN    | **0.84**    | **-6.8** | **+5.2** |
|                          | 2016         | MARINE Gốc  | 0.60        | -26.3    | +19.5    |
|                          |              | CẢI TIẾN    | **0.87**    | **+4.1** | **-3.8** |
|                          | 2017         | MARINE Gốc  | 0.59        | -29.0    | +21.1    |
|                          |              | CẢI TIẾN    | **0.85**    | **-5.5** | **+4.6** |
| Sông Cái Phan Rang       | 2016         | MARINE Gốc  | 0.61        | -25.8    | +17.9    |
| (Trạm Tân Mỹ)            |              | CẢI TIẾN    | **0.85**    | **+5.8** | **-4.1** |
|                          | 2018         | MARINE Gốc  | 0.64        | -22.4    | +15.6    |
|                          |              | CẢI TIẾN    | **0.89**    | **-3.9** | **+3.2** |

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

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

  1. Triệt tiêu hoàn toàn hiện tượng "đỉnh lũ giả / đỉnh lũ lệch pha": Mô hình MARINE gốc do tích lũy đại số lưu lượng nên sóng lũ truyền về trạm cửa ra nhanh hơn thực tế từ $2 - 4\text{ giờ}$ và biên độ đỉnh lũ bị đánh giá thấp từ $20% - 32%$. Mô hình cải tiến với sóng động học 1 chiều phi tuyến đã tái tạo chính xác dạng đường quá trình lũ, đưa sai số đỉnh lũ ($PE_p$) về dưới $\pm 7%$ và sai số thời gian đỉnh lũ $\Delta T_p$ về mức $0 - 1\text{ giờ}$.
  2. Định lượng chính xác vai trò cắt/chậm lũ của hồ chứa vùng núi: Tại lưu vực sông Dinh Ninh Hòa, hồ Đá Bàn giữ vai trò chi phối hoàn toàn thủy văn hạ du. Mô hình cải tiến đã mô phỏng thành công dung tích trữ và quá trình xả tràn tự do trong trận lũ 2010 và 2016, giải quyết tình trạng dự báo vượt quá thực tế (overestimation) mà các mô hình không có mô đun hồ thường mắc phải.
  3. Hiện tượng mưa cực đoan phân hóa theo sườn đón gió: Kết quả nội suy IDW kết hợp hồi quy phi tuyến độ cao chỉ ra rằng tại các đỉnh Hòn Giao ($2038\text{ m}$) và Bi Đúp ($2280\text{ m}$), lượng mưa giờ trong các trận bão lịch sử lớn hơn từ $1,8 - 2,5\text{ lần}$ so với lượng mưa đo được tại các trạm đồng bằng ven biển (Nha Trang, Phan Rang). Việc đưa trường mưa hiệu chỉnh theo cao trình vào mô hình là yếu tố quyết định nâng chỉ số NSE từ mức trung bình ($0,60$) lên mức xuất sắc ($>0,85$).
  4. Khả năng ứng dụng dự báo nghiệp vụ thời gian thực: Thử nghiệm dự báo dòng chảy lũ hạn $24\text{h}$ và $48\text{h}$ trên lưu vực sông Cái Nha Trang (các trạm Đồng Trăng, Diên Phú) trong mùa lũ năm 2020 đạt độ tin cậy rất cao, sai số mực nước đỉnh dự báo nằm trong phạm vi cho phép của quy chuẩn dự báo quốc gia.

Implications đa chiều

  • Về mặt học thuật: Chứng minh tính khả thi của việc "lai ghép mềm" giữa mô hình phân bố vật lý 2D bề mặt và mô hình thủy lực 1D lòng dẫn tinh giản, tạo hướng đi mới cho bài toán PUB tại các quốc gia đang phát triển.
  • Về mặt quản lý và chính sách: Cung cấp công cụ tính toán thủy văn độc lập, tin cậy cho Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các tỉnh Khánh Hòa, Ninh Thuận trong việc ban hành lệnh vận hành hồ chứa liên hồ và cảnh báo sơ tán dân cư vùng trũng hạ du.
  • Về mặt kinh tế - xã hội: Giảm thiểu thiệt hại về người và tài sản do lũ lụt bất ngờ tại các đô thị ven biển Nam Trung Bộ (Nha Trang, Ninh Hòa, Phan Rang - Tháp Chàm), đồng thời hỗ trợ vận hành tối ưu hóa phát điện và tích nước phục vụ mùa khô.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions):

  1. Chưa xét đến tác động của thủy triều và nước dâng do bão: Mô hình sóng động học 1 chiều hoạt động trên nguyên lý sóng truyền một chiều từ thượng lưu xuống hạ lưu theo trọng lực, do đó không thể mô phỏng được hiện tượng sóng triều truyền ngược (backwater effects) và ngập lụt vùng cửa sông giáp biển.
  2. Giới hạn trong điều kiện dòng chảy bùn cát và sạt lở đất: Mô hình giả định lòng dẫn cố định, chưa tích hợp mô đun xói lở, bồi lắng lòng dẫn và vận chuyển bùn cát hạt thô - vốn diễn ra rất khốc liệt trong các trận lũ quét vùng núi Nam Trung Bộ.
  3. Độ phân giải thời gian của số liệu khí tượng: Việc phụ thuộc vào số liệu mưa đo giờ thủ công hoặc trạm tự động đôi khi gặp gián đoạn do sự cố đường truyền viễn thông trong mưa bão lớn.

Định hướng nghiên cứu tương lai (5-10 năm tới):

  • Tích hợp mô hình MARINE cải tiến với mô hình thủy lực 2 chiều ngập lụt (2D Hydrodynamic Inundation Model) cho vùng hạ lưu đồng bằng ven biển chịu ảnh hưởng tương tác lũ sông - triều biển.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning / Deep Learning - LSTM, Physics-Informed Neural Networks - PINN) để tự động hóa quá trình tối ưu bộ thông số và hiệu chỉnh sai số dự báo thời gian thực.
  • Đồng hóa dữ liệu viễn thám thời gian thực từ radar thời tiết Doppler và sản phẩm mưa vệ tinh GSMAP/GPM có hiệu chỉnh tức thời bằng mạng lưới trạm mưa tự động đo mặt đất.

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

  • Ảnh hưởng học thuật: Luận án đóng góp một bộ mã nguồn cải tiến có tính mở cao, cung cấp tài liệu tham khảo chuyên sâu cho các chương trình đào tạo sau đại học ngành Thủy văn học, Tài nguyên nước và Biến đổi khí hậu tại Việt Nam.
  • Chuyển giao công nghệ nghiệp vụ: Kết quả nghiên cứu đã được chuyển giao, thử nghiệm trực tiếp tại Đài Khí tượng Thủy văn khu vực Nam Trung Bộ và Trung tâm Dự báo Khí tượng Thủy văn Quốc gia, hỗ trợ đắc lực cho các dự báo viên trong việc ra bản tin cảnh báo lũ sớm.
  • Đóng góp cho khung thể chế quản lý thiên tai: Cung cấp cơ sở khoa học để rà soát, cập nhật Quy trình vận hành liên hồ chứa trên lưu vực sông Cái Nha Trang và sông Cái Phan Rang theo Quyết định của Thủ tướng Chính phủ.

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

|                  CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                   |
| Nhóm đối tượng                 | Giá trị và Lợi ích nhận được             |
| Nghiên cứu sinh & Giới Viện /  | Khung phương pháp luận tích hợp mô hình  |
| Trường đại học chuyên ngành    | mã nguồn mở; giải pháp khắc phục PUB     |
| Cơ quan Dự báo Khí tượng       | Bộ công cụ dự báo lũ phân bố hạn 24h-48h |
| Thủy văn Quốc gia & Khu vực    | vận hành nhanh, không đòi hỏi mặt cắt    |
| Đơn vị Quản lý & Khai thác     | Mô đun diễn toán dòng chảy qua hồ hỗ trợ |
| Hồ chứa Thủy lợi / Thủy điện   | điều tiết tích nước và xả lũ an toàn     |
| Chính quyền địa phương &       | Bản đồ cảnh báo ngập lụt cục bộ phục vụ  |
| Ban Chỉ huy PCTT & TKCN        | công tác sơ tán dân và bảo vệ hạ tầng    |

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?

Đóng góp độc đáo nhất là việc thiết lập thành công mô hình tích hợp cấp độ hợp nhất giữa Lý thuyết dòng chảy tràn 2D sườn dốc Saint-Venant/Green-Ampt và Lý thuyết sóng động học 1D phi tuyến cho mạng lưới sông phân cấp. Luận án đã mở rộng lý thuyết sóng động học cổ điển của Lighthill & Whitham (1955) bằng cách nhúng trực tiếp nguồn gia nhập sườn dốc phân bố liên tục theo không gian từ từng ô lưới vào phương trình vi phân liên tục của từng đoạn sông, giải quyết triệt để bài toán biến dạng sóng lũ mà không cần dữ liệu trắc dọc, trắc ngang lòng dẫn.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu trước đây trên thế giới và Việt Nam?

So với mô hình MARINE gốc tại Pháp (Estupina, 2000; Bessière, 2005) vốn chỉ diễn toán một trục sông chính đơn giản hoặc phải kết nối ngoài với các phần mềm thương mại đắt tiền như HEC-RAS hay MIKE 11, luận án đã tự chủ xây dựng cấu trúc mô hình liền khối (monolithic model architecture). So với các nghiên cứu của Nguyễn Lan Châu (2001), Bùi Đình Lập (2006) tại Việt Nam từng kết nối MARINE với IMECH-1D theo cấp độ nối tiếp (loose coupling), giải pháp tích hợp của luận án loại bỏ hoàn toàn độ trễ truyền dữ liệu trung gian, cho phép mô phỏng tức thời mạng lưới sông nhánh và hồ chứa trong cùng một vòng lặp thời gian.

3. Phát hiện bất ngờ nhất trong quá trình phân tích số liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là tại lưu vực sông Dinh Ninh Hòa và Cái Nha Trang, nếu sử dụng phương pháp đa giác Thiessen để tính mưa bình quân lưu vực, đỉnh lũ tính toán luôn bị lệch pha trễ từ $2 - 4\text{ giờ}$ so với thực tế do không bắt được tâm mưa dịch chuyển trên đỉnh núi cao. Tuy nhiên, khi áp dụng mô hình IDW kết hợp hồi quy phi tuyến độ cao, gradient mưa theo phương thẳng đứng đã tái lập chính xác thời điểm bão hòa đất vùng thượng nguồn, giúp mô hình tái tạo đỉnh lũ trùng khớp với thời gian thực đo với độ sai lệch $\Delta T_p \le 1\text{ giờ}$.

4. Luận án có cung cấp đầy đủ giao thức để tái lặp nghiên cứu (Replication Protocol) không?

Có. Toàn bộ cấu trúc tệp tin đầu vào chuẩn ASCII cho địa hình (DEM), bản đồ thổ nhưỡng, bản đồ lớp phủ thảm thực vật, cấu trúc tệp dữ liệu mạng lưới sông, ma trận kết nối nút - nhánh, bảng quan hệ dung tích hồ chứa VFZ và mã nguồn các hàm toán học cải tiến (Qluoisong, Runge-Kutta 3, IDW-Elevation) đều được hệ thống hóa chi tiết trong thuyết minh và hệ thống phụ lục bảng biểu của luận án.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?

Luận án định hình lộ trình phát triển gồm 3 giai đoạn: (1) Chuẩn hóa gói phần mềm mã nguồn mở MARINE-VN có giao diện đồ họa người dùng (GUI); (2) Tích hợp mô-đun thủy lực ngập lụt 2 chiều giải bằng phương pháp thể tích hữu hạn cho vùng hạ lưu đồng bằng; (3) Xây dựng hệ thống tự động hóa hoàn toàn từ khâu thu nhận dữ liệu viễn thám/mưa radar đến khâu đồng hóa dữ liệu và xuất bản tin dự báo lũ trực tuyến.

Kết luận

  1. Luận án đã giải quyết thành công bài toán khoa học cấp thiết về mô phỏng và dự báo thủy văn phân bố trên các lưu vực sông thiếu số liệu mặt cắt hình học và mưa không gian tại khu vực Nam Trung Bộ.
  2. Cải tiến toàn diện mô hình MARINE thông qua 3 trụ cột kỹ thuật đột phá: Tích hợp mô hình sóng động học 1 chiều phi tuyến cho mạng lưới sông phức tạp; Tích hợp mô đun diễn toán dòng chảy qua hồ chứa bằng thuật toán Runge - Kutta bậc 3; Tích hợp công cụ nội suy mưa không gian IDW kết hợp tương quan phi tuyến độ cao địa hình.
  3. Kiểm chứng thực nghiệm xuất sắc trên 3 lưu vực sông trọng điểm (Sông Cái Nha Trang, Sông Dinh Ninh Hòa, Sông Cái Phan Rang) với chuỗi lũ lịch sử nhiều năm, nâng cao hệ số Nash-Sutcliffe ($NSE$) lên mức $0,84 - 0,91$, kiểm soát sai số đỉnh lũ trong phạm vi hẹp $\pm 7%$.
  4. Ứng dụng thành công vào nghiệp vụ dự báo lũ thời gian thực hạn $24\text{h} - 48\text{h}$ trong mùa lũ năm 2020 trên lưu vực sông Cái Nha Trang, khẳng định tính ổn định và độ tin cậy vượt trội trong điều kiện vận hành thực tế.
  5. Mở ra hướng nghiên cứu mới đầy triển vọng trong việc phát triển các bộ mô hình thủy văn phân bố mã nguồn mở phù hợp với đặc thù tự nhiên và điều kiện kinh tế - xã hội của Việt Nam, đóng góp thiết thực cho sự nghiệp phòng chống thiên tai và phát triển bền vững đất nước.