Tổng quan về luận án
Thiên tai lũ quét (flash flood) tại các lưu vực sông suối miền núi là một trong những hiểm họa tự nhiên khốc liệt, có sức tàn phá lớn nhất đối với tính mạng, tài sản và cơ sở hạ tầng dân sinh. Theo số liệu thống kê của Tổng cục Phòng chống Thiên tai Việt Nam, trung bình mỗi năm cả nước xảy ra từ 10 đến 15 trận lũ quét; riêng các tỉnh miền núi phía Bắc trong hai thập kỷ qua đã ghi nhận trên 300 trận lũ quét khốc liệt. Điển hình tại tỉnh Hà Giang, giai đoạn 2015–2019 ghi nhận 15 trận lũ quét làm 70 người chết, 1 người mất tích, 82 người bị thương và phá hủy gần 19.853 ha hoa màu. Tính chất hình thành dòng chảy xiết đột ngột trong thời gian tập trung nước cực ngắn ($< 3 - 6$ giờ) tại các lưu vực địa hình dốc phân cắt mạnh khiến các mô hình cảnh báo tĩnh truyền thống bộc lộ nhiều điểm nghẽn nghiêm trọng.
flowchart TD
A["Dữ liệu số DEM & Viễn thám"] --> B["Phân chia tiểu lưu vực & Trích xuất hình thái"]
C["Khảo sát thực địa 6 dấu hiệu tràn bờ"] --> D["Thủy lực Manning Q = f(Z)"]
D --> E["Hồi quy thực nghiệm Lưu lượng tràn bờ Qbf"]
F["Mạng lưới đo mưa tự động thời gian thực"] --> G["Mô hình toán thủy văn CTM (Java)"]
G --> H["Cập nhật trạng thái mặt đệm & Dòng chảy ban đầu"]
E & H --> I["Tính toán Ngưỡng mưa kích hoạt FFG"]
J["Mô hình dự báo thời tiết số trị GEM"] --> K["Đánh giá Chỉ số đe dọa lũ quét FFT"]
I & K --> L["Hệ thống cảnh báo sớm trực tuyến cấp Tiểu lưu vực"]
Khoảng trống nghiên cứu (research gap) cốt lõi nằm ở việc thiếu vắng một phương pháp luận khả thi để lượng hóa lưu lượng tràn bờ ($Q_{bf}$ - Bankfull Discharge) và xác định ngưỡng mưa định hướng sinh lũ quét ($FFG$ - Flash Flood Guidance) theo thời gian thực tại các lưu vực nhỏ miền núi không có trạm quan trắc thủy văn liên tục (ungauged basins). Các nghiên cứu trước đây (như dự án VNFFGS hay các bản đồ phân vùng tiềm năng lũ quét FFPI) thường chỉ dựa vào các kịch bản tĩnh, áp dụng công thức kinh nghiệm ngoại suy chưa qua kiểm chứng thực địa chi tiết, hoặc phụ thuộc vào các mô hình dự báo mưa đơn nguồn với sai số lớn. Luận án tiến sĩ với đề tài "Nghiên cứu phương pháp cảnh báo lũ quét cho lưu vực nhỏ miền núi và áp dụng thử nghiệm cho 2 lưu vực Nậm Ly và Nà Nhùng, tỉnh Hà Giang" giải quyết bài toán này thông qua hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:
- Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xác định chính xác lưu lượng tràn bờ $Q_{bf}$ tại các lưu vực sông suối miền núi không có số liệu quan trắc thủy văn liên tục?
- Câu hỏi nghiên cứu 2 (RQ2): Cơ chế toán học nào cho phép cập nhật liên tục độ ẩm đất và trạng thái mặt đệm theo thời gian thực để xác định động ngưỡng mưa $FFG$?
- Câu hỏi nghiên cứu 3 (RQ3): Mô hình tích hợp giữa lưu lượng tràn bờ thực nghiệm, mô hình thủy văn tự phát triển (CTM) và dự báo mưa đa quy mô (GEM) có nâng cao thời gian dự dẫn và độ chính xác của chỉ số đe dọa lũ quét ($FFT$) hay không?
- Giả thuyết nghiên cứu 1 (H1): Mối quan hệ giữa lưu lượng tràn bờ $Q_{bf}$ và các biến số hình thái lưu vực (diện tích $F$, độ dốc lưu vực $S_{lv}$, chiều dài sông chính $L_c$, hệ số tổn thất $CN$) tuân theo quy luật hàm lũy thừa đa biến có ý nghĩa thống kê ($p < 0.05$).
- Giả thuyết nghiên cứu 2 (H2): Ngưỡng mưa kích hoạt lũ quét không phải là một hằng số tĩnh mà là một hàm biến thiên phi tuyến phụ thuộc vào trạng thái bão hòa độ ẩm của lưu vực trước trận mưa.
Khung lý thuyết của nghiên cứu được xây dựng trên sự giao thoa giữa Thuyết địa mạo dòng chảy (Fluvial Geomorphology Theory của Leopold & Wolman), Động lực học dòng thấm phi bão hòa (Horton, Green-Ampt, SCS-CN) và Lý thuyết đường lũ đơn vị tích lũy (Unit Hydrograph Theory). Phạm vi không gian của luận án tập trung tại hai lưu vực thử nghiệm Nậm Ly và Nà Nhùng (tỉnh Hà Giang), theo dõi và kiểm chứng qua 3 đợt lũ quét - ngập lụt điển hình trong mùa mưa lũ năm 2020 (đợt ngày 13-14/6/2020, ngày 05-08/7/2020 và ngày 19-22/7/2020). Nghiên cứu mang lại bước đột phá khi rút ngắn chu kỳ cập nhật dữ liệu cảnh báo xuống từng phút, mở rộng thời gian dự dẫn cảnh báo trước từ 30 đến 120 phút.
Literature Review và Positioning
Tổng quan y văn quốc tế cho thấy các hướng tiếp cận cảnh báo lũ quét đã trải qua nhiều giai đoạn phát triển với những cuộc tranh luận học thuật sâu sắc. Trường phái thứ nhất tiếp cận theo hướng tĩnh dựa trên Chỉ số tiềm năng lũ quét ($FFPI$ - Flash Flood Potential Index) do Jeffrey Zogg và Kevin Deitsch (2013), Greg Smith (2003) và Brewster (2009) khởi xướng tại Hoa Kỳ. Nhóm nghiên cứu này sử dụng công nghệ GIS chồng lớp 4 biến số tĩnh (độ dốc, thổ nhưỡng, lớp phủ thực vật, mật độ rừng) để phân vùng nguy cơ. Tuy nhiên, Forestieri (2016) và Hapuarachchi et al. (2011) đã chỉ trích mạnh mẽ phương pháp này vì bỏ qua hoàn toàn trạng thái động của lưu vực (độ ẩm đất ban đầu và dung tích trữ tức thời), dẫn đến tình trạng sai lệch nghiêm trọng: một trận mưa nhỏ trên nền đất bão hòa có thể gây lũ quét khốc liệt, trong khi mưa lớn trên nền lưu vực khô cằn lại bị giữ lại bởi tổn thất thấm ban đầu.
Trường phái thứ hai tập trung vào ngưỡng mưa kích hoạt, tiêu biểu là Phương pháp Đường tới hạn ($CL$ - Critical Line) của Bộ Xây dựng Nhật Bản (2005) áp dụng tại núi Unzen-Fugen, mô hình SABO tại núi lửa Merapi (Indonesia), và các nghiên cứu mô phỏng Monte Carlo kết hợp mô hình HEC-HMS của Golian et al. (2010) tại lưu vực Madarsoo (Iran). Mặc dù các phương pháp này đã lượng hóa được mối tương quan giữa cường độ mưa giờ và mưa lũy tích hoạt động, chúng vẫn xem ngưỡng mưa là một đường bao cố định, chưa liên kết động với hình thái lòng dẫn và cơ chế tràn bờ cục bộ.
| Tiêu chí so sánh |
Hệ thống ALERT & FFMP (Hoa Kỳ / WMO) |
Phương pháp Đường tới hạn CL (Nhật Bản / Indonesia) |
Mô hình VNFFGS (Việt Nam - 2012-2017) |
Mô hình Tích hợp CTM - $Q_{bf}$ (Luận án) |
| Cơ chế xác định ngưỡng |
Ngưỡng mưa FFG dựa trên lưu vực trung bình từ Radar Doppler |
Đường ranh giới tới hạn ($CL$) giữa mưa giờ và mưa tích lũy |
Ngưỡng FFG tính toán từ mô hình WRF cho tiểu lưu vực 10-30 $\text{km}^2$ |
Ngưỡng FFG động tính từ $Q_{bf}$ thực nghiệm và mô hình CTM |
| Phương pháp xác định $Q_{bf}$ |
Phương trình hồi quy vùng quy mô lớn (USGS / Lumia 2006) |
Gián tiếp qua tần suất dòng chảy bùn cát Sabo |
Công thức kinh nghiệm chưa kiểm định chi tiết miền Bắc |
Điều tra 6 chỉ số thực địa + Thủy lực Manning + Hồi quy đa biến |
| Nguồn dữ liệu mưa |
Đa cảm biến (Radar + Trạm tự động + Vệ tinh) |
Trạm đo mưa tự ghi chuyên dụng 10 phút |
Mô hình số trị WRF đơn lẻ (kém chính xác khi mưa $> 30\text{ mm}$) |
Tích hợp trạm đo tự động liên tục + Mô hình GEM đa quy mô |
| Độ phân giải cảnh báo |
Cấp tiểu lưu vực / Hạt |
Cấp sườn dốc cục bộ ven núi lửa |
Cấp tiểu lưu vực lớn miền Trung - Tây Nguyên |
Chi tiết đến từng tiểu lưu vực, cấp xã và cấp thôn bản |
Trường phái thứ ba khai thác ngưỡng mưa định hướng sinh lũ quét ($FFG$) kết hợp lưu lượng tràn bờ ($Q_{bf}$) theo định nghĩa của Leopold (1964), Williams (1978), Georgakakos (1984) và NOAA (1992). Nghiên cứu của Emmett (2001, 2004) tại dãy núi Rocky (Hoa Kỳ) chứng minh $Q_{bf}$ gắn liền với quá trình hình thành lòng dẫn tự nhiên, có chu kỳ lặp lại khoảng 1.5 năm ($Q_{1.5}$). Tuy nhiên, tại các lưu vực vùng cao có độ dốc lớn và thành phần hạt thô, chu kỳ lặp lại này biến thiên phức tạp từ 1 đến 3 năm. Tại Việt Nam, hệ thống VNFFGS (2012-2017) đã ứng dụng FFG nhưng chủ yếu tập trung cho miền Trung và Tây Nguyên, sử dụng dữ liệu mưa dự báo WRF vốn có sai số lớn đối với các đợt mưa rào cực đoan cục bộ $> 30\text{ mm/h}$.
Luận án định vị chính xác tại điểm giao thoa giữa việc khắc phục sự thiếu hụt số liệu đo đạc mặt cắt thực tế và việc cá thể hóa động lực học dòng chảy cho lưu vực sông miền núi Bắc Bộ. Bằng việc tích hợp trực tiếp 6 dấu hiệu nhận biết mực nước tràn bờ trên thực địa với phương trình hồi quy đa biến và công cụ toán thủy văn CTM tự lập trình, nghiên cứu đã khắc phục triệt để hạn chế cảnh báo khống ở thượng lưu và thiếu chính xác ở hạ lưu.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án mang lại những đóng góp nền tảng cho lý thuyết thủy văn địa mạo và mô hình hóa dòng chảy lưu vực nhỏ:
- Mở rộng lý thuyết cân bằng động lòng dẫn (Dynamic Equilibrium Channel Theory): Công trình của Leopold & Wolman (1957) và Williams (1978) chủ yếu xem lưu lượng tràn bờ $Q_{bf}$ là hàm đơn biến của diện tích lưu vực hứng nước $F$. Nghiên cứu này chứng minh rằng tại các lưu vực núi cao bị chia cắt mạnh, $Q_{bf}$ chịu sự chi phối đồng thời của tổ hợp đa biến không gian bao gồm độ dốc lưu vực ($S_{lv}$), mật độ lưới sông ($n$), chiều dài dòng chính ($L_c$) và đặc trưng tổn thất thấm mặt đệm ($CN$).
- Cơ chế hóa định nghĩa toán học của Ngưỡng mưa định hướng ($FFG$): Luận án chuyển dịch khái niệm định tính của WMO thành bài toán tối ưu tiếp xúc động:
$$FFG = \arg\min_{R_e} \left{ \max_{t} \left[ Q_{t}^{hiện_trạng}(P_{thực_đo}) + Q_{t}^{bổ_sung}(R_e) \right] = Q_{bf} \right}$$
Biểu thức này khẳng định $FFG$ là lượng mưa hiệu quả bổ sung ($R_e$) tối thiểu rơi trên lưu vực trong một khoảng thời gian nhất định sao cho đỉnh của đường quá trình dòng chảy tổng hợp vừa chạm đúng ngưỡng lưu lượng tràn bờ $Q_{bf}$ tại cửa ra.
Khung phân tích độc đáo
Khung phân tích của luận án kết hợp ba trụ cột lý thuyết: (1) Thủy lực lòng dẫn mở (Open Channel Hydraulics) qua phương trình Manning, (2) Động lực học thấm và dòng chảy mặt (Infiltration and Surface Runoff Kinetics) qua mô hình số SCS-CN và đường cong nước rút Horton/Phillip, và (3) Thủy văn toán diễn toán sóng lũ (Hydrologic Flood Routing) qua phương pháp Muskingum và bể chứa tuyến tính.
classDiagram
class BasinComponent {
+Double Area
+Double CurveNumber
+calculateInfiltration(Horton, GreenAmpt)
+calculateSurfaceRunoff(SCS_UnitHydrograph)
+calculateBaseflow(RecessionCurve)
}
class ReachComponent {
+Double Length
+Double Slope
+Double Roughness
+routeMuskingum()
+routeLinearReservoir()
}
class JunctionComponent {
+List upstreamInflows
+computeTotalDischarge()
}
class WarningEngine {
+Double BankfullDischarge
+calculateFFG()
+calculateThreatFFT()
}
BasinComponent --> JunctionComponent : Cung cấp dòng chảy tiểu lưu vực
ReachComponent --> JunctionComponent : Diễn toán lũ dọc sông
JunctionComponent --> WarningEngine : So sánh đỉnh lũ với Qbf
Khung phân tích thiết lập điều kiện biên chặt chẽ: áp dụng cho các lưu vực nhỏ diện tích từ vài $\text{km}^2$ đến dưới $500\text{ km}^2$, mạng lưới sông suối cấp 2-3 có độ dốc đáy lớn ($S_s > 0.02$), lòng dẫn bồi tích cuội sỏi pha cát, loại trừ các lưu vực có sự điều tiết nhân tạo quy mô lớn của các hồ chứa thủy điện đa mục tiêu.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án tuân thủ thế giới quan thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng phương pháp nghiên cứu hỗn hợp kết hợp chặt chẽ giữa đo đạc thực nghiệm ngoài thực địa, viễn thám GIS và mô phỏng số học thời gian thực (real-time numerical modeling).
Thiết kế nghiên cứu đa tầng bao gồm:
- Tầng 1 (Hình thái và Thủy lực thực địa): Khảo sát chi tiết trắc ngang, trắc dọc tại các vị trí xung yếu ven khu dân cư và công trình hạ tầng. Xác định độ cao mực nước tràn bờ ($SBF$) thông qua 6 chỉ báo địa mạo tiêu chuẩn theo Bent (2013), Blanton (2010) và USDA (2014):
- Cao trình bãi sông nơi dòng chảy bắt đầu tràn bờ;
- Cao trình cao nhất tại các điểm bồi tụ trầm tích trong lòng dẫn;
- Điểm gãy thay đổi độ dốc của mái bờ sông;
- Điểm phân tầng thay đổi kích thước cấp phối hạt vật liệu đáy/bờ;
- Cao trình đỉnh của các hốc xói hàm ếch (undercut);
- Ranh giới chuyển tiếp sinh thái thực vật (từ vùng không có thực vật ngập nước sang vùng thảm phủ cố định).
- Tầng 2 (Mô phỏng thủy văn phân bố): Tự xây dựng mô hình toán CTM (Custom Hydrological Model) viết bằng ngôn ngữ Java hướng đối tượng, tích hợp trực tiếp vào nhân hệ thống cảnh báo thay vì sử dụng các phần mềm đóng gói thương mại.
Quy trình nghiên cứu rigorous
Quy trình thu thập và xử lý số liệu đạt chuẩn kiểm định khắt khe:
- Khảo sát mặt cắt và xây dựng quan hệ $Q = f(Z)$: Tại mỗi vị trí khảo sát (như các mặt cắt N4, N5, N6, N7, Q5, Q14 trên lưu vực Nậm Ly và Nà Nhùng), thiết lập đường cong quan hệ mực nước - lưu lượng dựa trên phương trình thủy lực Manning:
$$Q = \frac{1}{n} F R^{2/3} \sqrt{S}$$
Trong đó, $F$ là diện tích ướt ($\text{m}^2$), $R$ là bán kính thủy lực ($\text{m}$), $S$ là độ dốc đáy lòng sông thực đo, $n$ là hệ số nhám lòng dẫn xác định theo bảng phân loại chuẩn của Barnes (1969).
- Khấu trừ tổn thất thấm và phân tách dòng ngầm trong mô hình CTM: Sử dụng phương pháp SCS Curve Number ($CN$) và phương trình cường độ thấm Horton:
$$K_t = K_c + (K_o - K_c) e^{-\alpha t}$$
Hoặc phương trình thấm Philip:
$$K_t = K_c + \frac{A}{t^{1/2}}$$
Trong đó $K_o, K_c$ là cường độ thấm ban đầu và thấm ổn định ($\text{mm/phút}$), $\alpha$ là hệ số triết giảm, $A$ là thông số đặc trưng cho loại đất. Dòng ngầm được phân tách tự động theo phương pháp đường cong nước rút tuyến tính (recession curve).
- Diễn toán lũ lòng dẫn: Áp dụng phương pháp Muskingum và bể chứa tuyến tính để chuyển đổi sóng lũ từ các sườn dốc tiểu lưu vực về các nút liên kết (junctions) và cửa ra (outlets).
Data và phân tích
Toàn bộ dữ liệu không gian được xử lý trên nền tảng DEM độ phân giải cao kết hợp bản đồ địa hình 1:10.000 do Bộ Tài nguyên và Môi trường ban hành, liên kết bản đồ phân loại đất và lớp phủ thực vật để xác định chỉ số $CN$ cho từng tiểu lưu vực thông qua mạng lưới đa giác Thiessen.
Các phương trình hồi quy xác định $Q_{bf}$ được kiểm định bằng phương pháp bình phương nhỏ nhất tổng quát (GLS). Kiểm định ý nghĩa thống kê của các biến độc lập thông qua giá trị $p$-value ($p < 0.05$). Mô hình CTM được hiệu chuẩn và kiểm định độ tin cậy thông qua hệ số hiệu quả Nash-Sutcliffe ($NSE$), sai số thể tích đỉnh lũ ($RE$) so sánh trực tiếp với kết quả mô phỏng từ mô hình chuẩn quốc tế HEC-HMS.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Xác lập phương trình hồi quy thực nghiệm đa biến $Q_{bf}$ cho vùng núi phía Bắc: Kết quả phân tích thống kê đa biến đã lượng hóa mức độ tác động của các yếu tố địa lý - hình thái đến lưu lượng tràn bờ. Thay vì chỉ phụ thuộc vào diện tích $F$, phương trình tối ưu tích hợp diện tích lưu vực, độ dốc lưu vực ($S_{lv}$) và chỉ số số đường cong dòng chảy ($CN$) mang lại hệ số xác định $R^2 > 0.88$ với $p < 0.01$, vượt trội hoàn toàn so với các công thức kinh nghiệm đơn biến truyền thống ($R^2 < 0.62$).
- Quy luật phi tuyến của ngưỡng mưa kích hoạt $FFG$ theo thời gian: Kết quả thử nghiệm trên 2 lưu vực Nậm Ly và Nà Nhùng chỉ ra rằng ngưỡng mưa $FFG$ biến động rất mạnh tùy thuộc vào lượng mưa tích lũy kỳ trước ($API$). Khi lưu vực ở trạng thái khô hạn, giá trị $FFG$ thời đoạn 1 giờ đạt mức $45 - 65\text{ mm/h}$; tuy nhiên, khi lớp đất mặt đã đạt trạng thái bão hòa ẩm sau các đợt mưa kéo dài 3-5 ngày, ngưỡng $FFG$ giảm mạnh xuống chỉ còn $15 - 20\text{ mm/h}$, giải thích hiện tượng lũ quét xảy ra ngay cả trong các cơn mưa có cường độ trung bình.
- Độ tương thích cao giữa mô hình tự phát triển CTM và HEC-HMS: Kết quả so sánh tính toán dòng chảy lũ tại cửa ra tiểu lưu vực Nậm Ly 1 giữa mô hình CTM (nền tảng Java) và phần mềm HEC-HMS cho thấy sự tương đồng cao về cả đỉnh lũ ($Q_{max}$), thời gian xuất hiện đỉnh ($T_p$) và tổng lượng dòng chảy lũ, với sai số đỉnh lũ $< 6.5%$ và chỉ số $NSE > 0.85$.
So sánh đường quá trình dòng chảy tại cửa ra Nậm Ly 1:
Lưu lượng (m³/s)
^
Mưa tích lũy (GEM + Trạm đo)
- Kiểm chứng thực tế thành công qua 3 trận lũ lịch sử năm 2020: Hệ thống đã tái hiện và cảnh báo chính xác diễn biến ngập lụt và lũ quét tại:
- Thị trấn Quảng Nguyên (lưu vực Nậm Ly) trong đợt lũ ngày 13-14/6/2020;
- Thị trấn Quảng Nguyên trong đợt mưa lũ ngày 05-08/7/2020;
- Thôn Cốc Nam, xã Bản Nhùng (lưu vực Nà Nhùng) sau trận lũ quét bùn đá khốc liệt ngày 19-22/7/2020.
Hệ thống cảnh báo thử nghiệm đã phát tín hiệu cảnh báo nguy cơ cao trước khi lũ quét xuất hiện trên thực địa từ 45 đến 90 phút, cung cấp cơ sở tin cậy cho công tác ứng phó khẩn cấp tại địa phương.
Implications đa chiều
- Về mặt lý luận học thuật: Công trình khẳng định tính đúng đắn của việc kết hợp thủy lực địa mạo thực địa với mô hình hóa thủy văn số trị phân bố, mở ra phương pháp luận chuẩn hóa cho bài toán nghiên cứu thủy văn tại các lưu vực không có trạm quan trắc (Prediction in Ungauged Basins - PUB).
- Về mặt phương pháp luận: Khẳng định tính tự chủ công nghệ thông qua việc xây dựng mô hình mã nguồn mở CTM bằng Java, cho phép can thiệp sâu vào cấu trúc thuật toán và dễ dàng nhúng vào các hệ thống máy chủ cảnh báo trực tuyến mà không bị giới hạn bởi bản quyền phần mềm thương mại.
- Về mặt thực tiễn và chính sách: Cung cấp cho Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tỉnh Hà Giang một công cụ trực quan, hỗ trợ ra quyết định phân cấp theo thời gian thực đến tận cấp xã/thôn bản (hiển thị trực quan bản đồ nguy cơ, cấp độ cảnh báo theo màu sắc, xuất bản tin tự động).
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:
- Phạm vi dữ liệu thực nghiệm: Số lượng mặt cắt khảo sát trực tiếp $Q_{bf}$ tập trung chủ yếu trên 2 lưu vực Nậm Ly và Nà Nhùng (Hà Giang). Bộ thông số phương trình hồi quy cần được tái chuẩn hóa khi áp dụng cho các vùng địa chất khác biệt (như vùng karst đá vôi Tây Bắc hoặc đất đỏ bazan Tây Nguyên).
- Độ phân giải của mô hình mưa dự báo GEM: Mô hình GEM mặc dù có nhiều ưu điểm nhưng đối với các ổ mây đối lưu nhiệt cực đoan cục bộ phát triển nhanh dưới 1 giờ, sai số về vị trí tâm mưa vẫn còn tồn tại.
- Chưa tích hợp biến động sạt lở bồi lấp lòng dẫn theo thời gian thực: Mô hình giả định mặt cắt lòng sông ổn định trong suốt trận lũ, chưa xét đến hiện tượng sạt lở bờ gây nghẽn dòng tạm thời (tạo đập bùn đá tự nhiên) rồi vỡ đập đột ngột.
Định hướng nghiên cứu tiếp theo:
- Tích hợp thêm dữ liệu Radar thời tiết phân giải cao (Doppler radar) và đồng hóa dữ liệu vệ tinh đa kênh để nâng cao độ chính xác của trường mưa dự báo tức thời (nowcasting);
- Mở rộng mạng lưới khảo sát thực địa 6 dấu hiệu tràn bờ cho toàn bộ các lưu vực trọng điểm vùng Đông Bắc và Tây Bắc;
- Phát triển mô-đun liên kết thủy động lực học bùn cát 2D để mô phỏng dòng lũ bùn đá hỗn hợp rắn - lỏng.
Tác động và ảnh hưởng
Nghiên cứu có sức lan tỏa sâu rộng trên nhiều bình diện:
- Học thuật: Đóng gó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 Kỹ thuật Tài nguyên nước, Thủy văn học và Quản lý Rủi ro Thiên tai; cung cấp cơ sở dữ liệu thực nghiệm quý giá cho các công bố quốc tế thuộc danh mục ISI/Scopus.
- Quản lý nhà nước: Chuyển giao công nghệ cho các cơ quan khí tượng thủy văn địa phương và Chi cục Thủy lợi, giúp nâng cao năng lực cảnh báo sớm theo Chiến lược quốc gia phòng, chống và giảm nhẹ thiên tai.
- Xã hội: Bảo vệ an toàn tính mạng và sinh kế cho các cộng đồng dân tộc thiểu số tại vùng sâu vùng xa, giảm thiểu thiệt hại kinh tế hàng chục tỷ đồng mỗi năm do lũ quét gây ra.
Đối tượng hưởng lợi
mindmap
root((HỆ SINH THÁI HƯỞNG LỢI))
Nghiên cứu sinh & Giới học thuật
Phương pháp luận PUB
Khung hồi quy Qbf đa biến
Mô hình toán CTM mã nguồn mở
Cơ quan Quản lý & Dự báo KTTV
Công cụ trực quan hóa WebGIS
Thời gian dự dẫn mở rộng 30-120 phút
Bản tin cảnh báo chi tiết cấp Thôn Bản
Chính quyền Địa phương & Cộng đồng
Kế hoạch sơ tán dân chủ động
Bảo vệ hạ tầng dân sinh xung yếu
Giảm thiểu thiệt hại tính mạng và tài sản
- Nghiên cứu sinh và Giới học thuật Thủy văn: Tiếp cận một phương pháp luận hoàn chỉnh về xác định lưu lượng tràn bờ thực nghiệm và công cụ mô phỏng CTM mã nguồn mở để phát triển các hướng nghiên cứu chuyên sâu.
- Cơ quan Dự báo Khí tượng Thủy văn và Quản lý Đê điều: Tiếp nhận công cụ cảnh báo lũ quét trực tuyến có khả năng tích hợp linh hoạt các nguồn dữ liệu mưa tự động, nâng cao chất lượng bản tin dự báo.
- Chính quyền cơ sở và Cộng đồng dân cư vùng cao: Hưởng lợi trực tiếp từ các thông tin cảnh báo sớm có độ tin cậy cao, giúp chủ động sơ tán người và tài sản trước khi lũ tràn về.
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 mở rộng Thuyết địa mạo thủy lực lòng dẫn (Fluvial Hydraulic Geometry Theory) của Leopold & Wolman (1957) và Williams (1978). Thay vì coi lưu lượng tràn bờ $Q_{bf}$ chỉ là một hàm đơn biến tĩnh theo diện tích lưu vực hoặc tần suất lặp lại cố định 1.5 năm, luận án đã chứng minh tại các lưu vực sông suối nhỏ miền núi dốc, $Q_{bf}$ là hàm phụ thuộc phi tuyến vào cấu trúc hình thái đa chiều ($F, S_{lv}, L_c, S_s$) và lớp phủ mặt đệm ($CN$), từ đó thiết lập cầu nối lý thuyết giữa hình thái lòng dẫn và động lực học ngưỡng mưa kích hoạt lũ quét ($FFG$).
2. Điểm mới về mặt phương pháp luận so với các công trình quốc tế và trong nước?
So với các nghiên cứu của Lumia (2006) tại New York (sử dụng 14 biến hồi quy vùng quy mô lớn) hay Golian et al. (2010) tại Iran (mô phỏng Monte Carlo trên HEC-HMS với ngưỡng mưa tĩnh), luận án đã:
- Tích hợp 6 chỉ báo địa mạo thực địa tiêu chuẩn (Bent 2013) để xác định chính xác mực nước tràn bờ $SBF$ tại các mặt cắt thực tế vùng núi Việt Nam;
- Tự phát triển hoàn toàn mô hình toán CTM bằng Java, cho phép liên kết động và cập nhật trạng thái bão hòa độ ẩm đất liên tục theo từng bước thời gian thực đo, khắc phục nhược điểm dự báo mưa sai lệch lớn của hệ thống VNFFGS.
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à sự suy giảm phi tuyến tính cực lớn của ngưỡng mưa $FFG$: trong cùng một tiểu lưu vực, khi độ ẩm đất chuyển từ trạng thái khô sang bão hòa hoàn toàn (sau chuỗi ngày mưa liên tục), ngưỡng mưa kích hoạt lũ quét có thể suy giảm hơn 65% (từ mức $> 60\text{ mm/h}$ xuống chỉ còn xấp xỉ $18\text{ mm/h}$). Điều này giải thích hiện tượng nhiều trận lũ quét kinh hoàng trong lịch sử xảy ra trong những trận mưa không quá lớn về cường độ tức thời nhưng diễn ra trên nền lưu vực đã ngậm no nước.
4. Luận án có cung cấp quy trình nhân bản (replication protocol) cho các khu vực khác không?
Có. Luận án xây dựng quy trình chuẩn 5 bước rõ ràng: (1) Trích xuất đặc trưng lưu vực từ DEM 1:10.000 và bản đồ thổ nhưỡng; (2) Khảo sát thực địa đo đạc trắc ngang/dọc và nhận diện 6 dấu hiệu tràn bờ; (3) Tính toán đường cong $Q=f(Z)$ bằng công thức Manning và thiết lập phương trình hồi quy $Q_{bf}$; (4) Cài đặt mô hình CTM, nạp tham số $CN$, Horton/Muskingum và kết nối trạm mưa tự động; (5) Thiết lập ngưỡng cảnh báo $FFT$ và kết nối giao diện WebGIS.
5. Định hướng chương trình nghiên cứu 10 năm tới được vạch ra như thế nào?
Chương trình nghiên cứu dài hạn bao gồm: (1) Phát triển thuật toán AI/Machine Learning để tự động nhận dạng mực nước tràn bờ từ ảnh viễn thám độ phân giải siêu cao (SAR/UAV); (2) Tích hợp mô hình dự báo mưa thời gian thực bằng Radar phân cực kép; (3) Mở rộng hệ thống cảnh báo tích hợp đa thiên tai (Lũ quét - Sạt lở đất - Lũ bùn đá) phủ kín toàn bộ 14 tỉnh trung du và miền núi phía Bắc.
Kết luận
Luận án tiến sĩ đã hoàn thành xuất sắc các mục tiêu đề ra với những dấu ấn khoa học nổi bật:
- Xây dựng thành công cơ sở khoa học và phương trình hồi quy thực nghiệm xác định lưu lượng tràn bờ $Q_{bf}$ cho các lưu vực nhỏ miền núi thiếu số liệu quan trắc, với độ tin cậy thống kê cao ($R^2 > 0.88, p < 0.01$).
- Phát triển độc lập phần mềm mô hình toán thủy văn CTM (Java) có cấu trúc module linh hoạt, mô phỏng chính xác các quá trình thấm, dòng chảy mặt và dòng ngầm, đạt độ tương thích cao với mô hình chuẩn HEC-HMS ($NSE > 0.85$).
- Lượng hóa thành công phương pháp xác định ngưỡng mưa định hướng sinh lũ quét ($FFG$) động theo thời gian thực, phản ánh chính xác hiện trạng trữ ẩm của lưu vực.
- Đánh giá và tích hợp hiệu quả mô hình mưa dự báo GEM, nâng cao năng lực tính toán chỉ số đe dọa lũ quét ($FFT$).
- Thử nghiệm thành công hệ thống cảnh báo sớm trực tuyến trên 2 lưu vực Nậm Ly và Nà Nhùng (Hà Giang), tái hiện chính xác 3 trận lũ quét lịch sử năm 2020 và mở rộng thời gian dự dẫn an toàn từ 30 đến 120 phút.
- Mở ra 3 hướng nghiên cứu mới: Mô hình hóa lũ bùn đá hai pha trên lưu vực dốc; Ứng dụng công nghệ Internet vạn vật (IoT) trong giám sát mực nước tự động giá thành thấp; và Đồng hóa dữ liệu Radar/Vệ tinh trong cảnh báo thiên tai tức thời.
